Refine
Year of publication
Document Type
- Article (2579) (remove)
Language
- English (2579) (remove)
Keywords
- Holocene (35)
- climate change (31)
- Earthquake source observations (20)
- Climate change (18)
- permafrost (18)
- Himalaya (17)
- Pollen (17)
- Tibetan Plateau (17)
- erosion (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2579) (remove)
Pre-existing weakness zones in the lithosphere such as transform faults/fracture zones and extinct mid-oceanic ridges have been suggested to facilitate subduction initiation in an intra-oceanic environment. Here, we propose that the additional forcing coming from the mantle suction flow is required to trigger the conversion of a fracture zone/transform fault into a converging plate boundary. This suction flow can be induced either from the slab remnants of former converging plate boundaries or/and from slabs of neighbouring active subduction zones. Using 2-D coupled thermo-mechanical models, we show that a sufficiently strong mantle flow is able to convert a fracture zone/transform fault into a subduction zone. However, this process is feasible only if the fracture zone/transform fault is very close to the mid-oceanic ridge. Our numerical model results indicate that time of subduction initiation depends on the velocity, domain size and location of mantle suction flow and age of the oceanic plate.
Accretionary orogens are considered major sites of formation of juvenile continental crust. In the central and southern Andes this is contradicted by two observations: siliciclastic fills of Paleozoic basins in the central Andean segment of the accretionary Terra Australis Orogen consist almost exclusively of shales and mature sandstones; and magmatic rocks connected to the Famatinian (Ordovician) and Late Paleozoic magmatic arcs are predominantly felsic and characterized by significant crustal contamination and strongly unradiogenic Nd isotope compositions. Evidence of juvenile crustal additions is scarce. We present laser ablation (LA)-ICPMS U-Pb ages and LA-MC-ICPMS Hf isotope data of detrital zircons from seven Devonian to Permian turbidite sandstones incorporated into a Late Paleozoic accretionary wedge at the western margin of Gondwana in northern Chile. The combination with Nd whole-rock isotope data permits us to trace the evolution of the South American continental crust through several Proterozoic and Paleozoic orogenic cycles. The analyzed detrital zircon spectra reflect all Proterozoic orogenic cycles representing the step-wise evolution of the accretionary SW Amazonia Orogenic System between 2.0 and 0.9 Ga, followed by the Terra Australis Orogen between 0.9 and 0.25 Ga. The zircon populations are characterized by two prominent maxima reflecting input from Sunsas (Grenville) age magmatic rocks (1.2-0.9 Ga) and from the Ordovician to Silurian Famatinian magmatic arc (0.52-0.42 Ga). Grains of Devonian age are scarce or absent from the analyzed zircon populations. The Hf isotopic compositions of selected dated zircons at the time of their crystallization (epsilon Hf-(T); T = 3.3-0.25 Ga) vary between -18 and +11. All sandstones have a significant juvenile component; between 20 and 50% of the zircons from each sedimentary rock have positive epsilon Hf-(T) and can be considered juvenile. The majority of the juvenile grains have Hf-depleted mantle model ages (Hf T-DM) between 1.55 and 0.8 Ga. the time of the Rondonia-San Ignacio and Sunsas orogenic events on the Amazonia craton. The corresponding whole-rock epsilon Nd-(T) values fot these same rocks are between -8 and -3 indicating a mixture of older evolved and juvenile sources. Nd-depleted mantle model ages (Nd T-DM*) are between 1.5 and 1.2 Ga and coincide broadly with the zircon Hf model ages. Our data indicate that the Paleo- and Mesoproterozoic SW Amazonia Orogenic System, and the subsequent Neoproterozoic and Paleozoic Terra Australis Orogen in the region of the central and southern Andes, developed following two markedly different patterns of accretionary orogenic crustal evolution. The SW Amazonia Orogenic System developed by southwestward growth over approximately 1.1 Ga through a combination of accretion of juvenile material and crustal recycling typical of the extensional or retreating mode of accretionary orogens. In contrast, the central Andean segment of the Terra Australis Orogen evolved from 0.9 to 0.25 Ga in the compressional or advancing mode in a relatively fixed position without the accretion of oceanic crustal units or large scale input of juvenile material to the orogenic crust. Here, recycling mainly of Mesoproterozoic continental crust has been the dominant process of crustal evolution. (C) 2009 Elsevier B.V. All rights reserved.
The formation of the supercontinent Pangaea during the Permo-Triassic gave rise to an extreme monsoonal climate (often termed 'mega-monsoon') that has been documented by numerous palaeo-records. However, considerable debate exists about the role of orbital forcing in causing humid intervals in an otherwise arid climate. To shed new light on the forcing of monsoonal variability in subtropical Pangaea, this study focuses on sediment facies and colour variability of playa and alluvial fan deposits in an outcrop from the late Carnian (ca 225 Ma) in the southern Germanic Basin, south-western Germany. The sediments were deposited against a background of increasingly arid conditions following the humid Carnian Pluvial Event (ca 234 to 232 Ma). The ca 2 center dot 4 Myr long sedimentary succession studied shows a tripartite long-term evolution, starting with a distal mud-flat facies deposited under arid conditions. This phase was followed by a highly variable playa-lake environment that documents more humid conditions and finally a regression of the playa-lake due to a return of arid conditions. The red-green (a*) and lightness (L*) records show that this long-term variability was overprinted by alternating wet/dry cycles driven by orbital precession and ca 405 kyr eccentricity, without significant influence of obliquity. The absence of obliquity in this record indicates that high-latitude forcing played only a minor role in the southern Germanic Basin during the late Carnian. This is different from the subsequent Norian when high-latitude signals became more pronounced, potentially related to the northward drift of the Germanic Basin. The recurring pattern of pluvial events during the late Triassic demonstrates that orbital forcing, in particular eccentricity, stimulated the occurrence and intensity of wet phases. It also highlights the possibility that the Carnian Pluvial Event, although most likely triggered by enhanced volcanic activity, may also have been modified by an orbital stimulus.
Linking deep seismic profiles with regional-scale gravity inversion is a powerful tool to deduce the architecture of rifted margins and their structural evolution. Here we map upper and lower crustal thicknesses of the northern South China Sea (SCS) margin in order to investigate the occurrence of depth-dependent crustal extension from the proximal to the distal margin. By comparing upper and lower crustal stretching factors, we find that the northern margin of the SCS is segmented in three parts: (1) sedimentary basins where upper crust is stretched more than lower crust, (2) distal margin where lower crust is stretched more than upper crust, (3) mostly proximal margin regions where the two layers have similar stretching factors. Our results suggest that sedimentary basins and distal margin prominently feature depth-dependent extension, however accommodated by different processes. While differential thinning within sedimentary basins appears to be governed by lateral pressure variations inducing lower crustal flow, we suggest the distal margin to be affected by a combination of mantle flow-induced lower crustal shearing and sequential fault activity during crustal hyper-extension.
The Younger Dryas event, which began approximately 12,900 years ago, was a period of rapid cooling in the Northern Hemisphere, driven by large-scale reorganizations of patterns of atmospheric and oceanic circulation(1-3). Environmental changes during this period have been documented by both proxy-based reconstructions(3) and model simulations(4), but there is currently no consensus on the exact mechanisms of onset, stabilization or termination of the Younger Dryas(5-8). Here we present high-resolution records from two sediment cores obtained from Lake Krakenes in western Norway and the Nordic seas. Multiple proxies from Lake Krakenes are indicative of rapid alternations between glacial growth and melting during the later Younger Dryas. Meanwhile, reconstructed sea surface temperature and salinity from the Nordic seas show an alternation between sea-ice cover and the influx of warm, salty North Atlantic waters. We suggest that the influx of warm water enabled the westerly wind systems to drift northward, closer to their present-day positions. The winds thus brought relatively warm maritime air to Northern Europe, resulting in rising temperatures and the melting of glaciers. Subsequent input of this fresh meltwater into the ocean spurred the formation of sea ice, which forced the westerly winds back to the south, cooling Northern Europe. We conclude that rapid alternations between these two states immediately preceded the termination of the Younger Dryas and the permanent transition to an interglacial state.
Aiming at the stimulation of intrinsic microbial activity, pulses of pure oxygen or pressurized air were recurrently injected into groundwater polluted with chlorobenzene. To achieve well-controlled conditions and intensive sampling, a large, vertical underground tank was filled with the local unconfined sandy aquifer material. In the course of two individual gas injections, one using pure oxygen and one using pressurized air, the mass transfer of individual gas species between trapped gas phase and groundwater was studied. Field data on the dissolved gas composition in the groundwater were combined with a kinetic model on gas dissolution and transport in porous media. Phase mass transfer of individual gas components caused a temporary enrichment of nitrogen, and to a lower degree of methane, in trapped gas leading to the formation of excess dissolved nitrogen levels downgradient from the dissolving gas phase. By applying a novel gas sampling method for dissolved gases in groundwater it was shown that dissolved nitrogen can be used as a partitioning tracer to indicate complete gas dissolution in porous media.
Thirteen N-butylpyridinium salts, including three monometallic [C4Py](2)[MCl4], nine bimetallic [C4Py](2)[(M1-xMxCl4)-M-a-Cl-b] and one trimetallic compound [C4Py](2)[(M1-y-zMyMz (c) Cl4)-M-a-M-b] (M=Co, Cu, Mn; x=0.25, 0.50 or 0.75 and y=z=0.33), were synthesized and their structure and thermal and electrochemical properties were studied. All compounds are ionic liquids (ILs) with melting points between 69 and 93 degrees C. X-ray diffraction proves that all ILs are isostructural. The conductivity at room temperature is between 10(-4) and 10(-8) S cm(-1). Some Cu-based ILs reach conductivities of 10(-2) S cm(-1), which is, however, probably due to IL dec. This correlates with the optical bandgap measurements indicating the formation of large bandgap semiconductors. At elevated temperatures approaching the melting points, the conductivities reach up to 1.47x10(-1) S cm(-1) at 70 degrees C. The electrochemical stability windows of the ILs are between 2.5 and 3.0 V.
Fifteen N-butylpyridinium salts - five monometallic [C4Py](2)[MBr4] and ten bimetallic [C4Py](2)[(M0.5M0.5Br4)-M-a-Br-b] (M=Co, Cu, Mn, Ni, Zn) - were synthesized, and their structures and thermal and electrochemical properties were studied. All the compounds are ionic liquids (ILs) with melting points between 64 and 101 degrees C. Powder and single-crystal X-ray diffraction show that all ILs are isostructural. The electrochemical stability windows of the ILs are between 2 and 3 V. The conductivities at room temperature are between 10(-5) and 10(-6) S cm(-1). At elevated temperatures, the conductivities reach up to 10(-4) S cm(-1) at 70 degrees C. The structures and properties of the current bromide-based ILs were also compared with those of previous examples using chloride ligands, which illustrated differences and similarities between the two groups of ILs.
The removal, redistribution, and transient storage of sediments in tectonically active mountain belts is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but unambiguous evidence for this potential link has been elusive so far. Here we combine 2D numerical experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks) with field data from intermontane basins in the NW Argentine Andes. Our modeling results suggest that thicker sedimentary fills (>0.7-1.0 km) may suppress basinal faulting processes, while thinner fills (<0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and basin excavation promotes renewed intrabasinal faulting. These results help to better understand the tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean Plateau in northwestern Argentina. For example, the Santa Maria and the Humahuaca basins record intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal deformation patterns in intermontane basins of tectonically active broken forelands. (C) 2018 Elsevier B.V. All rights reserved.
Sedimentary basins in the interior of orogenic plateaus can provide unique insights into the early history of plateau evolution and related geodynamic processes. The northern sectors of the Iranian Plateau of the Arabia-Eurasia collision zone offer the unique possibility to study middle-late Miocene terrestrial clastic and volcaniclastic sediments that allow assessing the nascent stages of collisional plateau formation. In particular, these sedimentary archives allow investigating several debated and poorly understood issues associated with the long-term evolution of the Iranian Plateau, including the regional spatio-temporal characteristics of sedimentation and deformation and the mechanisms of plateau growth. We document that middle-late Miocene crustal shortening and thickening processes led to the growth of a basement-cored range (Takab Range Complex) in the interior of the plateau. This triggered the development of a foreland-basin (Great Pari Basin) to the east between 16.5 and 10.7Ma. By 10.7Ma, a fast progradation of conglomerates over the foreland strata occurred, most likely during a decrease in flexural subsidence triggered by rock uplift along an intraforeland basement-cored range (Mahneshan Range Complex). This was in turn followed by the final incorporation of the foreland deposits into the orogenic system and ensuing compartmentalization of the formerly contiguous foreland into several intermontane basins. Overall, our data suggest that shortening and thickening processes led to the outward and vertical growth of the northern sectors of the Iranian Plateau starting from the middle Miocene. This implies that mantle-flow processes may have had a limited contribution toward building the Iranian Plateau in NW Iran.
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 idea that climatically modulated erosion may impact orogenic processes has challenged geoscientists for decades. Although modeling studies and physical calculations have provided a solid theoretical basis supporting this interaction, to date, field-based work has produced inconclusive results. The central-western Alborz Mountains in the northern sectors of the Arabia-Eurasia collision zone constitute a promising area to explore these potential feedbacks. This region is characterized by asymmetric precipitation superimposed on an orogen with a history of spatiotemporal changes in exhumation rates, deformation patterns, and prolonged, km-scale base-level changes. Our analysis suggests that despite the existence of a strong climatic gradient at least since 17.5 Ma, the early orogenic evolution (from similar to 36 to 9-6 Ma) was characterized by decoupled orographic precipitation and tectonics. In particular, faster exhumation and sedimentation along the more arid southern orogenic flank point to a north-directed accretionary flux and underthrusting of Central Iran. Conversely, from 6 to 3 Ma, erosion rates along the northern orogenic flank became higher than those in the south, where they dropped to minimum values. This change occurred during a similar to 3-Myr-long, km-scale base-level lowering event in the Caspian Sea. We speculate that mass redistribution processes along the northern flank of the Alborz and presumably across all mountain belts adjacent to the South Caspian Basin and more stable areas of the Eurasian plate increased the sediment load in the basin and ultimately led to the underthrusting of the Caspian Basin beneath the Alborz Mountains. This underthrusting in turn triggered a new phase of northward orogenic expansion, transformed the wetter northern flank into a new pro-wedge, and led to the establishment of apparent steady-state conditions along the northern orogenic flank (i.e., rock uplift equal to erosion rates). Conversely, the southern mountain front became the retro-wedge and experienced limited tectonic activity. These observations overall raise the possibility that mass-distribution processes during a pronounced erosion phase driven by base-level changes may have contributed to the inferred regional plate-tectonic reorganization of the northern Arabia-Eurasia collision during the last similar to 5 Ma. (C) 2015 Elsevier B.V. All rights reserved.
The Central Pontides of N Turkey represents a mobile orogenic belt of the southern Eurasian margin that experienced several phases of exhumation associated with the consumption of different branches of the Neo-Tethys Ocean and the amalgamation of continental domains. Our new low-temperature thermochronology data help to constrain the timing of these episodes, providing new insights into associated geodynamic processes. In particular, our data suggest that exhumation occurred at (1) similar to 110 to 90Ma, most likely during tectonic accretion and exhumation of metamorphic rocks from the subduction zone; (2) from similar to 60 to 40Ma, during the collision of the Kirehir and Anatolide-Tauride microcontinental domains with the Eurasian margin; (3) from similar to 0 to 25Ma, either during the early stages of the Arabia-Eurasia collision (soft collision) when the Arabian passive margin reached the trench, implying 70 to 530km of subduction of the Arabian passive margin, or during a phase of trench advance predating hard collision at similar to 20Ma; and (4) similar to 11Ma to the present, during transpression associated with the westward motion of Anatolia. Our findings document the punctuated nature of fault-related exhumation, with episodes of fast cooling followed by periods of slow cooling or subsidence, the role of inverted normal faults in controlling the Paleogene exhumation pattern, and of the North Anatolian Fault in dictating the most recent pattern of exhumation.
The southern foreland basin of the Alborz Mountains of northern Iran is characterized by an approximately 7.3-km-thick sequence of Miocene sedimentary rocks, constituting three basin-wde coarsening-upward units spanning a period of 10(6)years. We assess available magnetostratigraphy, paleoclimatic reconstructions, stratal architecture, records of depositional environments, and sediment-provenance data to characterize the relationships between tectonically-generated accommodation space (A) and sediment supply (S). Our analysis allows an inversion of the stratigraphy for particular forcing mechanisms, documenting causal relationships, and providing a basis to decipher the relative contributions of tectonics and climate (inferred changes in precipitation) in controlling sediment supply to the foreland basin. Specifically, A/S>1, typical of each basal unit (17.5-16.0, 13.8-13.1 and 10.3-9.6Ma), is associated with sharp facies retrogradation and reflects substantial tectonic subsidence. Within these time intervals, arid climatic conditions, changes in sediment provenance, and accelerated exhumation in the orogen suggest that sediment supply was most likely driven by high uplift rates. Conversely, A/S<1 (13.8 and 13.8-11Ma, units 1, and 2) reflects facies progradation during a sharp decline in tectonic subsidence caused by localized intra-basinal uplift. During these time intervals, climate continued to be arid and exhumation active, suggesting that sediment supply was again controlled by tectonics. A/S<1, at 11-10.3Ma and 9-6-7.6Ma (and possibly 6.2; top of units 2 and 3), is also associated with two episodes of extensive progradation, but during wetter phases. The first episode appears to have been linked to a pulse in sediment supply driven by an increase in precipitation. The second episode reflects a balance between a climatically-induced increase in sediment supply and a reduction of subsidence through the incorporation of the proximal foreland into the orogenic wedge. This in turn caused an expansion of the catchment and a consequent further increase in sediment supply.
A poorly understood lag time of 15-20 m.y. exists between the initial Arabia-Eurasia continental collision in late Eocene to early Oligocene time and the acceleration of tectonic and sedimentary processes across the collision zone in the early to late Miocene. The late Eocene to Miocene-Pliocene clastic and shallow-marine sedimentary rocks of the Kond, Eyvanekey, and Semnan Basins in the Alborz Mountains (northern Iran) offer the possibility to track the evolution of this orogen in the framework of collision processes. A transition from volcaniclastic submarine deposits to shallow-marine evaporites and terrestrial sediments occurred shortly after 36 Ma in association with reversals in sediment provenance, strata tilting, and erosional unroofing. These events followed the termination of subduction arc magmatism and marked a changeover from an extensional to a contractional regime in response to initiation of continental collision with the subduction of stretched Arabian lithosphere. This early stage of collision produced topographic relief associated with shallow foreland basins, suggesting that shortening and tectonic loading occurred at low rates. Starting from the early Miocene (17.5 Ma), flexural subsidence in response to foreland basin initiation occurred. Fast sediment accumulation rates and erosional unroofing trends point to acceleration of shortening by the early Miocene. We suggest that the lag time between the initiation of continental collision (36 Ma) and the acceleration of regional deformation (20-17.5 Ma) reflects a two-stage collision process, involving the "soft" collision of stretched lithosphere at first and "hard" collision following the arrival of unstretched Arabian continental litho sphere in the subduction zone.
With this paper, we assess the present-day conductive thermal field of the Glueckstadt Graben in NW Germany that is characterized by large salt walls and diapirs structuring the graben fill. We use a finite element method to calculate the 3D steady-state conductive thermal field based on a lithosphere-scale 3D structural model that resolves the first-order structural characteristics of the graben and its underlying lithosphere. Model predictions are validated against measured temperatures in six deep wells. Our investigations show that the interaction of thickness distributions and thermal rock properties of the different geological layers is of major importance for the distribution of temperatures in the deep subsurface of the Glueckstadt Graben. However, the local temperatures may result from the superposed effects of different controlling factors. Especially, the upper sedimentary part of the model exhibits huge lateral temperature variations, which correlate spatially with the shape of the thermally highly conductive Permian salt layer. Variations in thickness and geometry of the salt cause two major effects, which provoke considerable lateral temperature variations for a given depth. (1) The "chimney effect" causes more efficient heat transport within salt diapirs. As a consequence positive thermal anomalies develop in the upper part and above salt structures, where the latter are covered by much less conductive sediments. In contrast, negative thermal anomalies are noticeable underneath salt structures. (2) The "thermal blanketing effect" is caused by thermally low conductive sediments that provoke the local storage of heat where these insulating sediments are present. The latter effect leads to both local and regional thermal anomalies. Locally, this translates to higher temperatures where salt margin synclines are filled with thick insulating clastic sediments. For the regional anomalies the cumulative insulating effects of the entire sediment fill results in a long-wavelength variation of temperatures in response to heat refraction effects caused by the contrast between insulating sediments and highly conductive crystalline crust. Finally, the longest wavelength of temperature variations is caused by the depth position of the isothermal lithosphere-asthenosphere boundary defining the regional variations of the overall geothermal gradient. We find that a conductive thermal model predicts observed temperatures reasonably well for five of the six available wells, whereas the steady-state conductive approach appears not to be valid for the sixth well.
Sedimentological, provenance, and detrital thermochronological results for basin fill at the modern deformation front of the northern Andes (6 degrees N latitude) provide a long-term, Eocene to Pliocene record of foreland-basin sedimentation along the Eastern Cordillera !Janos basin boundary in Colombia. Lithofacies assemblages and paleocurrent orientations in the upward-coarsening, 5-km-thick succession of the Nunchia syncline reveal a systematic shift from craton-derived, shallow-marine distal foreland (back-bulge) accumulation in the Mirador Formation, to orogen-sourced, deltaic, and coastalinfluenced sedimentation of the distal to medial foreland (foredeep) in the Carbonera and Leon Formations, to anastomosing fluvial and distributive braided fluvial megafan systems of the proximal foreland (foredeep to wedge-top) basin in the lower and upper Guayabo Formation. These changes in depositional processes and sediment dispersal are supported by up-section variations in detrital zircon U-Pb and (U-Th)/He ages that record exhumation of evolving, compartmentalized sediment source areas in the Eastern Cordillera. The data are interpreted in terms of a progressive eastward advance in foldand-thrust deformation, with late Eocene Oligocene deformation in the axial zone of the Eastern Cordillera along the western edge of Floresta basin (Soapaga thrust), early Miocene reactivation (inversion) of the eastern margin of the Mesozoic rift system (Pajarito and Guaicaramo thrusts), and middle late Miocene propagation of a footwall shortcut fault (Vopal thrust) that created the Nunchia syncline in a wedge-top (piggyback) setting of the eastern foothills along the transition from the Eastern Cordillera to Harms foreland basin. Collectively, the data presented here for the frontal Eastern Cordillera define a general in-sequence pattern of eastwardadvancing fold-and-thrust deformation during Cenozoic east-west shortening in the Colombian Andes.
Groundwater recharge (GWR) is one of the most challenging water fluxes to estimate, as it relies on observed data that are often limited in many developing countries.
This study developed an innovative water budget method using satellite products for estimating the spatially distributed GWR at monthly and annual scales in tropical wet sedimentary regions despite cloudy conditions.
The distinctive features proposed in this study include the capacity to address 1) evapotranspiration estimations in tropical wet regions frequently overlaid by substantial cloud cover; and 2) seasonal root-zone water storage estimations in sedimentary regions prone to monthly variations.
The method also utilises satellite-based information of the precipitation and surface runoff. The GWR was estimated and validated for the hydrologically contrasting years 2016 and 2017 over a tropical wet sedimentary region located in North-eastern Brazil, which has substantial potential for groundwater abstraction.
This study showed that applying a cloud-cleaning procedure based on monthly compositions of biophysical data enables the production of a reasonable proxy for evapotranspiration able to estimate groundwater by the water budget method.
The resulting GWR rates were 219 (2016) and 302 (2017) mm yr(-1), showing good correlations (CC = 0.68 to 0.83) and slight underestimations (PBIAS =-13 to-9%) when compared with the referenced estimates obtained by the water table fluctuation method for 23 monitoring wells. Sensitivity analysis shows that water storage changes account for +19% to-22% of our monthly evaluation.
The satellite-based approach consistently demonstrated that the consideration of cloud-cleaned evapotranspiration and root-zone soil water storage changes are essential for a proper estimation of spatially distributed GWR in tropical wet sedimentary regions because of their weather seasonality and cloudy conditions.
The lack of process-based classification procedures may lead to unrealistic hyetograph design due to complex oscillation of rainfall depths when assimilated at high temporal resolutions. Four consecutive years of sub-hourly rainfall data were assimilated in three study areas (Guaraira, GEB, Sao Joao do Cariri, CEB, and Aiuaba, AEB) under distinct climates (very hot semi-arid and tropical wet). This study aimed to define rainfall events (for Minimum Inter-event Time, MIT, and Minimum Rainfall Depth, MRD, equal to 30 min and 1.016 mm, respectively), classify their hyetograph types (rectangular, R, unimodal with left-skewed, UL, right-skewed, UR, and centred peaks, UC, bimodal, B, and shapeless, SL), and compare their key rainfall properties (frequency, duration, depth, rate and peak). A rain pulse aggregation process allowed for reshaping SL-events for six different time spans varying from 2 to 30 min. The results revealed that the coastal area held predominantly R-events (64% events and 49% rainfall depth), in western semi-arid prevailed UL-events (57% events and 63% rainfall depth), whereas in eastern semi-arid mostly were R-events (61% events and 30% rainfall depth) similar to coastal area. It is concluded that each cloud formation type had important effects on hyetograph properties, differentiating them even within the same climate.
The Eastern Mediterranean is the most seismically active region in Europe due to the complex interactions of the Arabian, African, and Eurasian tectonic plates. Deformation is achieved by faulting in the brittle crust, distributed flow in the viscoelastic lower-crust and mantle, and Hellenic subduction, but the long-term partitioning of these mechanisms is still unknown. We exploit an extensive suite of geodetic observations to build a kinematic model connecting strike-slip deformation, extension, subduction, and shear localization across Anatolia and the Aegean Sea by mapping the distribution of slip and strain accumulation on major active geological structures. We find that tectonic escape is facilitated by a plate-boundary-like, translithospheric shear zone extending from the Gulf of Evia to the Turkish-Iranian Plateau that underlies the surface trace of the North Anatolian Fault. Additional deformation in Anatolia is taken up by a series of smaller-scale conjugate shear zones that reach the upper mantle, the largest of which is located beneath the East Anatolian Fault. Rapid north-south extension in the western part of the system, driven primarily by Hellenic Trench retreat, is accommodated by rotation and broadening of the North Anatolian mantle shear zone from the Sea of Marmara across the north Aegean Sea, and by a system of distributed transform faults and rifts including the rapidly extending Gulf of Corinth in central Greece and the active grabens of western Turkey. Africa-Eurasia convergence along the Hellenic Arc occurs at a median rate of 49.8mm yr(-1) in a largely trench-normal direction except near eastern Crete where variably oriented slip on the megathrust coincides with mixed-mode and strike-slip deformation in the overlying accretionary wedge near the Ptolemy-Pliny-Strabo trenches. Our kinematic model illustrates the competing roles the North Anatolian mantle shear zone, Hellenic Trench, overlying mantle wedge, and active crustal faults play in accommodating tectonic indentation, slab rollback and associated Aegean extension. Viscoelastic flow in the lower crust and upper mantle dominate the surface velocity field across much of Anatolia and a clear transition to megathrust-related slab pull occurs in western Turkey, the Aegean Sea and Greece. Crustal scale faults and the Hellenic wedge contribute only a minor amount to the large-scale, regional pattern of Eastern Mediterranean interseismic surface deformation.
In Germany, the irrigation sector accounts for only 1% of water use. In recent years, however, this sector has attracted more attention due to the occurrence of severe drought periods. Irrigation scheduling systems could support adaptation strategies but little is known about current providers, performance and users. In this study we aimed to depict the current situation of the existence and functioning of irrigation scheduling systems available in Germany. Six methods were identified and assessed based on direct interviews with end-users and a comparative analysis. The results showed a positive feedback from the users. However, the recommendations were rarely implemented, while only the seasonal irrigation requirement was considered to support actual water abstraction. These results were corroborated by the comparative analysis. Five of the six irrigation scheduling systems estimated the seasonal irrigation amount consistently, while wider differences were found by looking at the irrigation season and at the number of irrigations. Overall, it is found that irrigation support systems are valuable tools for supporting adaptation strategies to fast changes in agro-environmental conditions. However, specific assessments based on real measurements should be considered in order to improve the performance of the systems and provide more consistent support to end-users. (c) 2019 John Wiley & Sons, Ltd.
An effective strategy for combining variance- and distribution-based global sensitivity analysis
(2020)
We present a new strategy for performing global sensitivity analysis capable to estimate main and interaction effects from a generic sampling design. The new strategy is based on a meaningful combination of varianceand distribution-based approaches. The strategy is tested on four analytic functions and on a hydrological model. Results show that the analysis is consistent with the state-of-the-art Saltelli/Jansen formula but to better quantify the interaction effect between the input factors when the output distribution is skewed. Moreover, the estimation of the sensitivity indices is much more robust requiring a smaller number of simulations runs. Specific settings and alternative methods that can be integrated in the new strategy are also discussed. Overall, the strategy is considered as a new simple and effective tool for performing global sensitivity analysis that can be easily integrated in any environmental modelling framework.
Soil moisture dynamics are affected by complex interactions among several factors. Understanding the relative importance of these factors is still an important challenge in the study of water fluxes and solute transport in unsaturated media. In this study, the spatio-temporal variability of surface soil moisture was investigated in a 10 ha flat cropped field located in northern Italy. Soil moisture was measured on a regular 50 x 50 m grid on seven dates during the growing season. For each measurement campaign, the spatial variability of the soil moisture was compared with the spatial variability of the soil texture and crop properties. In particular, to better understand the role of the vegetation, the spatio-temporal variability of two different parameters - leaf area index and crop height - was monitored on eight dates at different crop development stages. Statistical and geostatistical analysis was then applied to explore the interactions between these variables. In agreement with other studies, the results show that the soil moisture variability changes according to the average value within the field, with the standard deviation reaching a maximum value under intermediate mean soil moisture conditions and the coefficient of variation decreasing exponentially with increasing mean soil moisture. The controls of soil moisture variability change according to the average soil moisture within the field. Under wet conditions, the spatial distribution of the soil moisture reflects the variability of the soil texture. Under dry conditions, the spatial distribution of the soil moisture is affected mostly by the spatial variability of the vegetation. The interaction between these two factors is more important under intermediate soil moisture conditions. These results confirm the importance of considering the average soil moisture conditions within a field when investigating the controls affecting the spatial variability of soil moisture. This study highlights the importance of considering the spatio-temporal variability of the vegetation in investigating soil moisture dynamics, especially under intermediate and dry soil moisture conditions. The results of this study have important implications in different hydrological applications, such as for sampling design, ranking stability application, indirect measurements of soil properties and model parameterisation.
Cosmic-Ray neutron sensing (CRS) is a unique approach to measure soil moisture at field scale filling the gap of current methodologies. However, CRS signal is affected by all the hydrogen pools on the land surface and understanding their relative importance plays an important role for the application of the method e.g., validation of remote sensing products and data assimilation. In this study, a soil moisture scaling approach is proposed to estimate directly the correct CRS soil moisture based on the soil moisture profile measured at least in one position within the field. The approach has the advantage to avoid the need to introduce one correction for each hydrogen contribution and to estimate indirectly all the related time-varying hydrogen pools. Based on the data collected in three crop seasons, the scaling approach shows its ability to identify and to quantify the seasonal biomass water equivalent. Additionally, the analysis conducted at sub-daily time resolution is able to quantify the daily vertical redistribution of the water biomass and the rainfall interception, showing promising applications of the CRS method also for these types of measurements. Overall, the study underlines how not only soil moisture but all the specific hydrological processes in the soil-plant-atmosphere continuum should be considered for a proper evaluation of the CRS signal. For this scope, the scaling approach reveals to be a simple and pragmatic analysis that can be easily extended to other experimental sites. (C) 2015 Elsevier B.V. All rights reserved.
Cosmic-ray neutron sensing (CRNS) is a promising proximal soil sensing technique to estimate soil moisture at intermediate scale and high temporal resolution. However, the signal shows complex and non-unique response to all hydrogen pools near the land surface, providing some challenges for soil moisture estimation in practical applications. Aims of the study were 1) to assess the uncertainty of CRNS as a stand-alone approach to estimate volumetric soil moisture in cropped field 2) to identify the causes of this uncertainty 3) and possible improvements. Two experimental sites in Germany were equipped with a CRNS probe and point-scale soil moisture network. Additional monitoring activities were conducted during the crop growing season to characterize the soil-plant systems. This data is used to identify and quantify the different sources of uncertainty (factors). An uncertainty analysis, based on Monte Carlo approach, is applied to propagate these uncertainties to CRNS soil moisture estimations. In addition, a sensitivity analysis based on the Sobol’ method is performed to identify the most important factors explaining this uncertainty. Results show that CRNS soil moisture compares well to the soil moisture network when these point-scale values are weighted to account for the spatial sensitivity of the signal and other sources of hydrogen (lattice water and organic carbon) are added to the water content. However, the performance decreases when CRNS is considered as a stand-alone method to retrieve the actual (non-weighted) volumetric soil moisture. The support volume (penetration depth and radius) shows also a considerable uncertainty, especially in relatively dry soil moisture conditions. Four of the seven factors analyzed (the vertical soil moisture profile, bulk density, incoming neutron correction and the calibrated parameter N0) were found to play an important role. Among the possible improvements identified, a simple correction factor based on vertical point-scale soil moisture profiles shows to be a promising approach to account for the sensitivity of the CRNS signal to the upper soil layers.
The present study proposes a General Probabilistic Framework (GPF) for uncertainty and global sensitivity analysis of deterministic models in which, in addition to scalar inputs, non-scalar and correlated inputs can be considered as well. The analysis is conducted with the variance-based approach of Sobol/Saltelli where first and total sensitivity indices are estimated. The results of the framework can be used in a loop for model improvement, parameter estimation or model simplification. The framework is applied to SWAP, a 113 hydrological model for the transport of water, solutes and heat in unsaturated and saturated soils. The sources of uncertainty are grouped in five main classes: model structure (soil discretization), input (weather data), time-varying (crop) parameters, scalar parameters (soil properties) and observations (measured soil moisture). For each source of uncertainty, different realizations are created based on direct monitoring activities. Uncertainty of evapotranspiration, soil moisture in the root zone and bottom fluxes below the root zone are considered in the analysis. The results show that the sources of uncertainty are different for each output considered and it is necessary to consider multiple output variables for a proper assessment of the model. Improvements on the performance of the model can be achieved reducing the uncertainty in the observations, in the soil parameters and in the weather data. Overall, the study shows the capability of the GPF to quantify the relative contribution of the different sources of uncertainty and to identify the priorities required to improve the performance of the model. The proposed framework can be extended to a wide variety of modelling applications, also when direct measurements of model output are not available.
Soil properties show high heterogeneity at different spatial scales and their correct characterization remains a crucial challenge over large areas. The aim of the study is to quantify the impact of different types of uncertainties that arise from the unresolved soil spatial variability on simulated hydrological states and fluxes. Three perturbation methods are presented for the characterization of uncertainties in soil properties. The methods are applied on the soil map of the upper Neckar catchment (Germany), as an example. The uncertainties are propagated through the distributed mesoscale hydrological model (mHM) to assess the impact on the simulated states and fluxes. The model outputs are analysed by aggregating the results at different spatial and temporal scales. These results show that the impact of the different uncertainties introduced in the original soil map is equivalent when the simulated model outputs are analysed at the model grid resolution (i.e. 500 m). However, several differences are identified by aggregating states and fluxes at different spatial scales (by subcatchments of different sizes or coarsening the grid resolution). Streamflow is only sensitive to the perturbation of long spatial structures while distributed states and fluxes (e.g. soil moisture and groundwater recharge) are only sensitive to the local noise introduced to the original soil properties. A clear identification of the temporal and spatial scale for which finer-resolution soil information is (or is not) relevant is unlikely to be universal. However, the comparison of the impacts on the different hydrological components can be used to prioritize the model improvements in specific applications, either by collecting new measurements or by calibration and data assimilation approaches. In conclusion, the study underlines the importance of a correct characterization of uncertainty in soil properties. With that, soil maps with additional information regarding the unresolved soil spatial variability would provide strong support to hydrological modelling applications.
In this paper, we examine the influence of the 27 October 2012, M-w 7.8 earthquake on landslide occurrence in the southern half of Haida Gwaii (formerly Queen Charlotte Islands), British Columbia, Canada. Our 1350 km(2) study area is undisturbed, primarily forested terrain that has not experienced road building or timber harvesting. Our inventory of landslide polygons is based on optical airborne and spaceborne images acquired between 2007 and 2018, from which we extracted and mapped 446 individual landslides (an average of 33 landslides per 100 km(2)). The landslide rate in years without major earthquakes averages 19.4 per year, or 1.4/100 km(2)/year, and the annual average area covered by non-seismically triggered landslides is 35 ha/year. The number of landslides identified in imagery closely following the 2012 earthquake, and probably triggered by it, is 244 or an average of about 18 landslides per 100 km(2). These landslides cover a total area of 461 ha. In the following years-2013-2016 and 2016-2018-the number of landslides fell, respectively, to 26 and 13.5 landslides per year. In non-earthquake years, most landslides happen on south-facing slopes, facing the prevailing winds. In contrast, during or immediately after the earthquake, up to 32% of the landslides occurred on north and northwest-facing slopes. Although we could not find imagery from the day after the earthquake, overview reconnaissance flights 10 and 16 days later showed that most of the landslides were recent, suggesting they were co-seismic.
A hydrochemical approach to quantify the role of return flow in a surface flow-dominated catchment
(2017)
Stormflow generation in headwater catchments dominated by subsurface flow has been studied extensively, yet catchments dominated by surface flow have received less attention. We addressed this by testing whether stormflow chemistry is controlled by either (a) the event-water signature of overland flow, or (b) the pre-event water signature of return flow. We used a high-resolution hydrochemical data set of stormflow and end-members of multiple storms in an end-member mixing analysis to determine the number of end-members needed to explain stormflow, characterize and identify potential end-members, calculate their contributions to stormflow, and develop a conceptual model of stormflow. The arrangement and relative positioning of end-members in stormflow mixing space suggest that saturation excess overland flow (26-48%) and return flow from two different subsurface storage pools (17-53%) are both similarly important for stormflow. These results suggest that pipes and fractures are important flow paths to rapidly release stored water and highlight the value of within-event resolution hydrochemical data to assess the full range and dynamics of flow paths.
End member mixing analysis (EMMA) is a commonly applied method to identify and quantify the dominant runoff producing sources of water. It employs tracers to determine the dimensionality of the hydrologic system. Many EMMA studies have been conducted using two to six tracers, with some of the main tracers being Ca, Na, Cl(-), water isotopes, and alkalinity. Few studies use larger tracer sets including minor trace elements such as Li, Rb, Sr, and Ba. None of the studies has addressed the question of the tracer set size and composition, despite the fact that these determine which and how many end members (EM) will be identified. We examine how tracer set size and composition affects the conceptual model that results from an EMMA. We developed an automatic procedure that conducts EMMA while iteratively changing tracer set size and composition. We used a set of 14 tracers and 9 EMs. The validity of the resulting conceptual models was investigated under the aspects of dimensionality, EM combinations, and contributions to stream water. From the 16,369 possibilities, 23 delivered plausible results. The resulting conceptual models are highly sensitive to the tracer set size and composition. The moderate reproducibility of EM contributions indicates a still missing EM. It also emphasizes that the major elements are not always the most useful tracers and that larger tracer sets have an enhanced capacity to avoid false conclusions about catchment functioning. The presented approach produces results that may not be apparent from the traditional approach and it is a first step to add the idea of statistical significance to the EMMA approach.
The spatial distribution of soil types is controlled by a set of environmental factors such as climate, organisms, parent material and topography as well as time and space. A change of these factors will lead to a change in the spatial distribution of soil types. In this study, we use a digital soil mapping approach to improve our knowledge about major soil type distributing factors in the steppe regions of Inner Mongolia (China) which currently undergo tremendous environmental change, e.g. climate and land use change. We use Random Forests in an effort to map Reference Soil Groups according to the World Reference Base for Soil Resources (WRB) in the Xilin River catchment. We benefit from the superior prediction capabilities of RF and additional interpretive results in order to identify the major environmental factors that control spatial patterns of soil types. The nine WRB soil groups that were identified and spatially predicted for the study area are Arenosol, Calcisol, Cambisol, Chernozem, Cryosol, Gleysol, Kastanozem, Phaeozem and Regosol.
Model and prediction performances of the RF model are high with an Out-of-Bag error of 51.6% for the model and a misclassification error for the predicted map of 28.9%. The main controlling factors of soil type distribution are land use, a set of topographic variables, geology and climate. However, land use and climate are of major importance and topography and geology are of minor importance. The visualizations of the predictions, the variable importance measures as result of RF and the comparisons of these with the spatial distribution of the environmental factors delivered additional, quantitative information of these controlling factors and revealed that intensively grazed areas are subjected to soil degradation. However, most of the area is still governed by natural soil forming processes which are driven by climate, topography and geology. Most importantly though, our study revealed that a shift towards warmer temperatures and lower precipitation regimes will lead to a change of the spatial distribution of RSGs towards steppe soils that store less carbon, i.e. a decrease of spatial extent of Phaeozems and an increase of spatial extent of Chernozems and Kastanozems.
Combinations of runoff characteristics are commonly used to represent distinct conceptual models of stormflow generation. In this study, three runoff characteristics: hydrograph response, time source of runoff water, and flow path are used to classify catchments. Published data from the scientific literature are used to provide evidence from small, forested catchments. Each catchment was assigned to one of the eight conceptual models, depending on the combination of quick/slow response, old/new water, and overland/subsurface flow. A standard procedure was developed to objectively diagnose the predominant conceptual model of stormflow generation for each catchment and assess its temporal and spatial support. The literature survey yielded 42 catchments, of which 30 catchments provide a complete set of qualitative runoff characteristics resulting in one of the eight conceptual models. The majority of these catchments classify as subsurface flow path dominated. No catchments were found for conceptual models representing combinations of quick response-new water-subsurface flow (SSF), slow-new-SSF, slow-old-overland flow (OF) nor new-slow-OF. Of the 30 qualitatively classified catchments, 24 provide a complete set of quantitative measures. In summary, the field support is strong for 19 subsurface-dominated catchments and is weak for 5 surface flow path dominated catchments (six catchments had insufficient quantitative data). Two alternative explanations exist for the imbalance of field support between the two flow path classes: (1) the selection of research catchments in past field studies was mainly to explain quick hydrograph response in subsurface dominated catchments; (2) catchments with prevailing subsurface flow paths are more common in nature. We conclude that the selection of research catchments needs to cover a wider variety of environmental conditions which should lead to a broader, and more widely applicable, spectrum of resulting conceptual models and process mechanisms. This is a prerequisite in studies where catchment organization and similarity approaches are used to develop catchment classification systems in order to regionalize stormflow.
Draft Genome Sequence of Nocardioides alcanivorans NGK65(T), a Hexadecane-Degrading Bacterium
(2022)
The Gram-positive bacterium Nocardioides alcanivorans NGK65(T) was isolated from plastic-polluted soil and cultivated on medium with polyethylene as the single carbon source. Nanopore sequencing revealed the presence of candidate enzymes for the biodegradation of polyethylene. Here, we report the draft genome of this newly described member of the terrestrial plastisphere.
This review presents a compositional database of primary anatectic granitoid magmas, entirely based on melt inclusions (MI) in high-grade metamorphic rocks. Although MI are well known to igneous petrologists and have been extensively studied in intrusive and extrusive rocks, MI in crustal rocks that have undergone anatexis (migmatites and granulites) are a novel subject of research. They are generally trapped along the heating path by peritectic phases produced by incongruent melting reactions. Primary MI in high-grade metamorphic rocks are small, commonly 5-10 pm in diameter, and their most common mineral host is peritectic garnet. In most cases inclusions have crystallized into a cryptocrystalline aggregate and contain a granitoid phase assemblage (nanogranitoid inclusions) with quartz, K-feldspar, plagioclase, and one or two mica depending on the particular circumstances. After their experimental remelting under high-confining pressure, nanogranitoid MI can be analyzed combining several techniques (EMP, LA-ICP-MS, NanoSIMS, Raman). The trapped melt is granitic and metaluminous to peraluminous, and sometimes granodioritic, tonalitic, and trondhjemitic in composition, in agreement with the different P-T-a(H2o) conditions of melting and protolith composition, and overlap the composition of experimental glasses produced at similar conditions. Being trapped along the up-temperature trajectory as opposed to classic MI in igneous rocks formed during down-temperature magma crystallization fundamental information provided by nanogranitoid MI is the pristine composition of the natural primary anatectic melt for the specific rock under investigation. So far similar to 600 nanogranitoid MI, coming from several occurrences from different geologic and geodynamic settings and ages, have been characterized. Although the compiled MI database should be expanded to other potential sources of crustal magmas, MI data collected so far can be already used as natural "starting-point" compositions to track the processes involved in formation and evolution of granitoid magmas.
We present a wavelet coherence method that is capable of displaying local coherence information between two seismic stations in the sense of a spectrogram. We have analyzed the vertical components of a 20-min-long time series from four stations that were situated in the seismic near field of Stromboli volcano. Typical volcanic seismic signals recorded in the near field of Stromboli volcano consist of continuous volcanic tremor superimposed on frequent Strombolian explosion signals. The tremor exhibits a banded and frequency-stable structure, whereas the broadband explosion signals span two or three frequency decades. We demonstrate that signals related to explosion earthquakes are strongly correlated within the network over 1.5 frequency decades. Using synthetic data, we show how coherent signal portions can be extracted out of noisy data using a coherence-filtering method. A time delay analysis using coherence information results in a coarse source location estimation that lies within the crater region. With the exception of randomly fluctuating coherence peaks, low correlations have been observed in the characteristic bands that are assumed to be generated by continuous tremor. In the low-frequency band that is related to the ocean microseisms (period approximate to 4-8 sec), we observe mostly high correlation that breaks down during the appearance of explosion earthquake signals. Based on further analysis using the inverse wavelet transformation, we propose a model that describes the breakdown phenomenon as a superposition of two independent events
We report the results of our investigations on the catchment area, surface sediments, and hydrology of the monsoonal Lonar Lake, central India. Our results indicate that the lake is currently stratified with an anoxic bottom layer, and there is a spatial heterogeneity in the sensitivity of sediment parameters to different environmental processes. In the shallow (0-5 m) near shore oxic-suboxic environments the lithogenic and terrestrial organic content is high and spatially variable, and the organics show degradation in the oxic part. Due to aerial exposure resulting from lake level changes of at least 3m, the evaporitic carbonates are not completely preserved. In the deep water (>5 m) anoxic environment the lithogenics are uniformly distributed and the delta C-13 is an indicator not only for aquatic vs. terrestrial plants but also of lake pH and salinity. The isotopic composition of the evaporites is dependent not only on the isotopic composition of source water (monsoon rainfall and stream inflow) and evaporation, but is also influenced by proximity to the isotopically depleted stream inflow. We conclude that in the deep water environment lithogenic content, and isotopic composition of organic matter can be used for palaeoenvironmental reconstruction.
The imagination of clearly separated core-shell structures is already outdated by the fact, that the nanoparticle core-shell structures remain in terms of efficiency behind their respective bulk material due to intermixing between core and shell dopant ions. In order to optimize the photoluminescence of core-shell UCNP the intermixing should be as small as possible and therefore, key parameters of this process need to be identified. In the present work the Ln(III) ion migration in the host lattices NaYF4 and NaGdF4 was monitored. These investigations have been performed by laser spectroscopy with help of lanthanide resonance energy transfer (LRET) between Eu(III) as donor and Pr(III) or Nd(III) as acceptor. The LRET is evaluated based on the Forster theory. The findings corroborate the literature and point out the migration of ions in the host lattices. Based on the introduced LRET model, the acceptor concentration in the surrounding of one donor depends clearly on the design of the applied core-shell-shell nanoparticles. In general, thinner intermediate insulating shells lead to higher acceptor concentration, stronger quenching of the Eu(III) donor and subsequently stronger sensitization of the Pr(III) or the Nd(III) acceptors. The choice of the host lattice as well as of the synthesis temperature are parameters to be considered for the intermixing process.
Seismic wide-angle data were collected along a 40-km-long profile centered at the geothermal research well GrSk 3/90 in the Northeast German Basin. Tomographic inversion of travel time data provided a velocity and a vertical velocity gradient model, indicative of Cenozoic to Pre-Permian sediments. Wide-angle reflections are modeled and interpreted as top Zechstein and top Pre-Permian. Changes in velocity gradients are interpreted as the transition from mechanical to chemical compaction at 2-3 km depth, and localized salt structures are imaged, suggesting a previously unknown salt pillow in the southern part of the seismic profile. The Zechstein salt shows decreased velocities in the adjacent salt pillows compared to the salt lows, which is confirmed by sonic log data. This decrease in velocity could be explained by the mobilization of less dense salt, which moved and formed the salt pillows, whereas the denser salt remained in place at the salt lows. We interpret a narrow subvertical low-velocity zone under the salt pillow at GrSk 3/ 90 as a fault in the deep Permian to Pre-Permian. This WNW-ESE trending fault influenced the location of the salt tectonics and led to the formation of a fault-bounded graben in the Rotliegend sandstones with optimal mechanical conditions for geothermal production. Thermal modeling showed that salt pillows are related to chimney effects, a decrease in temperature, and increasing velocity. The assumed variations in salt lithology, density, and strain must thus be even higher to compensate for the temperature effect.
Temperature is a key factor controlling plant growth and vitality in the temperate climates of the mid-latitudes like in vast parts of the European continent. Beyond the effect of average conditions, the timings and magnitudes of temperature extremes play a particularly crucial role, which needs to be better understood in the context of projected future rises in the frequency and/or intensity of such events. In this work, we employ event coincidence analysis (ECA) to quantify the likelihood of simultaneous occurrences of extremes in daytime land surface temperature anomalies (LSTAD) and the normalized difference vegetation index (NDVI). We perform this analysis for entire Europe based upon remote sensing data, differentiating between three periods corresponding to different stages of plant development during the growing season. In addition, we analyze the typical elevation and land cover type of the regions showing significantly large event coincidences rates to identify the most severely affected vegetation types. Our results reveal distinct spatio-temporal impact patterns in terms of extraordinarily large co-occurrence rates between several combinations of temperature and NDVI extremes. Croplands are among the most frequently affected land cover types, while elevation is found to have only a minor effect on the spatial distribution of corresponding extreme weather impacts. These findings provide important insights into the vulnerability of European terrestrial ecosystems to extreme temperature events and demonstrate how event-based statistics like ECA can provide a valuable perspective on environmental nexuses.
The Seismic Hazard Inferred from Tectonics based on the Global Strain Rate Map (SHIFT_GSRM) earthquake forecast was designed to provide high-resolution estimates of global shallow seismicity to be used in seismic hazard assessment. This model combines geodetic strain rates with global earthquake parameters to characterize long-term rates of seismic moment and earthquake activity. Although SHIFT_GSRM properly computes seismicity rates in seismically active continental regions, it underestimates earthquake rates in subduction zones by an average factor of approximately 3. We present a complementary method to SHIFT_GSRM to more accurately forecast earthquake rates in 37 subduction segments, based on the conservation of moment principle and the use of regional interface seismicity parameters, such as subduction dip angles, corner magnitudes, and coupled seismogenic thicknesses. In seven progressive steps, we find that SHIFT_GSRM earthquake-rate underpredictions are mainly due to the utilization of a global probability function of seismic moment release that poorly captures the great variability among subduction megathrust interfaces. Retrospective test results show that the forecast is consistent with the observations during the 1 January 1977 to 31 December 2014 period. Moreover, successful pseudoprospective evaluations for the 1 January 2015 to 31 December 2018 period demonstrate the power of the regionalized earthquake model to properly estimate subduction-zone seismicity.
The first step in the estimation of probabilistic seismic hazard in a region commonly consists of the definition and characterization of the relevant seismic sources. Because in low-seismicity regions seismicity is often rather diffuse and faults are difficult to identify, large areal source zones are mostly used. The corresponding hypothesis is that seismicity is uniformly distributed inside each areal seismic source zone. In this study, the impact of this hypothesis on the probabilistic hazard estimation is quantified through the generation of synthetic spatial seismicity distributions. Fractal seismicity distributions are generated inside a given source zone and probabilistic hazard is computed for a set of sites located inside this zone. In our study, the impact of the spatial seismicity distribution is defined as the deviation from the hazard value obtained for a spatially uniform seismicity distribution. From the generation of a large number of synthetic distributions, the correlation between the fractal dimension D and the impact is derived. The results show that the assumption of spatially uniform seismicity tends to bias the hazard to higher values. The correlation can be used to determine the systematic biases and uncertainties for hazard estimations in real cases, where the fractal dimension has been determined. We apply the technique in Germany (Cologne area) and in France (Alps).
Ground-motion prediction equations (GMPE) are essential in probabilistic seismic hazard studies for estimating the ground motions generated by the seismic sources. In low-seismicity regions, only weak motions are available during the lifetime of accelerometric networks, and the equations selected for the probabilistic studies are usually models established from foreign data. Although most GMPEs have been developed for magnitudes 5 and above, the minimum magnitude often used in probabilistic studies in low-seismicity regions is smaller. Disaggregations have shown that, at return periods of engineering interest, magnitudes less than 5 may be contributing to the hazard. This paper presents the testing of several GMPEs selected in current international and national probabilistic projects against weak motions recorded in France (191 recordings with source-site distances up to 300 km, 3:8 <= M-w <= 4:5). The method is based on the log-likelihood value proposed by Scherbaum et al. (2009). The best-fitting models (approximately 2:5 <= LLH <= 3:5) over the whole frequency range are the Cauzzi and Faccioli (2008), Akkar and Bommer (2010), and Abrahamson and Silva (2008) models. No significant regional variation of ground motions is highlighted, and the magnitude scaling could be the predominant factor in the control of ground-motion amplitudes. Furthermore, we take advantage of a rich Japanese dataset to run tests on randomly selected low-magnitude subsets, and confirm that a dataset of similar to 190 observations, the same size as the French dataset, is large enough to obtain stable LLH estimates. Additionally we perform the tests against larger magnitudes (5-7) from the Japanese dataset. The ranking of models is partially modified, indicating a magnitude scaling effect for some of the models, and showing that extrapolating testing results obtained from low-magnitude ranges to higher magnitude ranges is not straightforward.
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
The results of this study clearly identify four key parameters controlling the estimation of probabilistic seismic hazard assessment (PSHA) in France in the framework of the Cornell-McGuire method. Results in terms of peak ground acceleration demonstrate the equally high impact, at all return periods, of the choice of truncation of the predicted ground-motion distribution (at + 2sigma) and of the choice between two different magnitude-intensity correlations. The choice of minimum magnitude (3.5/4.5) on hazard estimates can have an important impact at small return periods (<1000 years), whereas the maximum magnitude (6.5/7.0), on the other hand, is not a key parameter even at large return periods (10,000 years). This hierarchy of impacts is maintained at lower frequencies down to 5 Hz. Below 5 Hz, the choice of the maximum magnitude has a much greater impact, whereas the impact due to the choice of the minimum magnitude disappears. Moreover, variability due to catalog uncertainties is also quantified; these uncertainties that underly all hazard results can engender as high a variability as the controlling parameters. Parameter impacts, calculated at the centers of each source zone, show a linear trend with the seismicity models of the zone, demonstrating the lack of contributions coming from neighboring zones. Indeed, the region of influence that contributes to the PSHA estimate at a given site decreases with increasing return periods. The resulting overall variability in hazard estimates due to input uncertainties is quantified through a logic tree, obtained coefficients of variation vary between 10% and 20%. Until better physical models are obtained, the uncertainty on hazard estimates may be reduced by working on an appropriate magnitude-intensity correlation
New age data have been obtained to time constrain the recent Quaternary volcanism of El Hierro (Canary Islands) and to estimate its recurrence rate. We have carried out Ar-40/Ar-39 geochronology on samples spanning the entire volcanostratigraphic sequence of the island and C-14 geochronology on the most recent eruption on the northeast rift of the island: 2280 +/- 30 yr BP. We combine the new absolute data with a revision of published ages onshore, some of which were identified through geomorphological criteria (relative data). We present a revised and updated chronology of volcanism for the last 33 ka that we use to estimate the maximum eruptive recurrence of the island. The number of events per year determined is 9.7 x 10(-4) for the emerged part of the island, which means that, as a minimum, one eruption has occurred approximately every 1000 years. This highlights the need of more geochronological data to better constrain the eruptive recurrence of El Hierro. (C) 2015 Elsevier Ltd. All rights reserved.
Beta diversity is a conceptual link between diversity at local and regional scales. Various additional methodologies of quantifying this and related phenomena have been applied. Among them, measures of pairwise (dis)similarity of sites are particularly popular. Undersampling, i.e. not recording all taxa present at a site, is a common situation in ecological data. Bias in many metrics related to beta diversity must be expected, but only few studies have explicitly investigated the properties of various measures under undersampling conditions. On the basis of an empirical data set, representing near-complete local inventories of the Lepidoptera from an isolated Pacific island, as well as simulated communities with varying properties, we mimicked different levels of undersampling. We used 14 different approaches to quantify beta diversity, among them dataset-wide multiplicative partitioning (i.e. true beta diversity') and pairwise site x site dissimilarities. We compared their values from incomplete samples to true results from the full data. We used these comparisons to quantify undersampling bias and we calculated correlations of the dissimilarity measures of undersampled data with complete data of sites. Almost all tested metrics showed bias and low correlations under moderate to severe undersampling conditions (as well as deteriorating precision, i.e. large chance effects on results). Measures that used only species incidence were very sensitive to undersampling, while abundance-based metrics with high dependency on the distribution of the most common taxa were particularly robust. Simulated data showed sensitivity of results to the abundance distribution, confirming that data sets of high evenness and/or the application of metrics that are strongly affected by rare species are particularly sensitive to undersampling. The class of beta measure to be used should depend on the research question being asked as different metrics can lead to quite different conclusions even without undersampling effects. For each class of metric, there is a trade-off between robustness to undersampling and sensitivity to rare species. In consequence, using incidence-based metrics carries a particular risk of false conclusions when undersampled data are involved. Developing bias corrections for such metrics would be desirable.
The seismicity pattern along the San Andreas fault near Parkfield and Cholame, California, varies distinctly over a length of only fifty kilometres. Within the brittle crust, the presence of frictionally weak minerals, fault-weakening high fluid pressures and chemical weakening are considered possible causes of an anomalously weak fault northwest of Parkfield(1-4). Non-volcanic tremor from lower-crustal and upper-mantle depths(5-7) is most pronounced about thirty kilometres southeast of Parkfield and is thought to be associated with high pore-fluid pressures at depth(8). Here we present geophysical evidence of fluids migrating into the creeping section of the San Andreas fault that seem to originate in the region of the uppermost mantle that also stimulates tremor, and evidence that along-strike variations in tremor activity and amplitude are related to strength variations in the lower crust and upper mantle. Interconnected fluids can explain a deep zone of anomalously low electrical resistivity that has been imaged by magnetotelluric data southwest of the Parkfield-Cholame segment. Near Cholame, where fluids seem to be trapped below a high-resistivity cap, tremor concentrates adjacent to the inferred fluids within a mechanically strong zone of high resistivity. By contrast, sub-vertical zones of low resistivity breach the entire crust near the drill hole of the San Andreas Fault Observatory at Depth, northwest of Parkfield, and imply pathways for deep fluids into the eastern fault block, coincident with a mechanically weak crust and the lower tremor amplitudes in the lower crust. Fluid influx to the fault system is consistent with hypotheses of fault-weakening high fluid pressures in the brittle crust.
Eight d-metal-containing N-butylpyridinium ionic liquids (ILs) with the nominal composition (C4Py)2[Ni0.5M0.5Cl4] or (C4Py)2[Zn0.5M0.5Cl4] (M = Cu, Co, Mn, Ni, Zn; C4Py = N-butylpyridinium) were synthesized, characterized, and investigated for their optical properties. Single crystal and powder X-ray analysis shows that the compounds are isostructural to existing examples based on other d-metal ions. Inductively coupled plasma optical emission spectroscopy measurements confirm that the metal/metal ratio is around 50 : 50. UV-Vis spectroscopy shows that the optical absorption can be tuned by selection of the constituent metals. Moreover, the compounds can act as an optical sensor for the detection of gases such as ammonia as demonstrated via a simple prototype setup.
In the Sri Lankan Highlands erosion and chemical weathering rates are among the lowest for global mountain denudation. In this tropical humid setting, highly weathered deep saprolite profiles have developed from high-grade metamorphic charnockite during spheroidal weathering of the bedrock. The spheroidal weathering produces rounded corestones and spalled rindlets at the rock-saprolite interface. We used detailed textural, mineralogical and chemical analyses to reconstruct the sequence of weathering reactions and their causes. The first mineral attacked by weathering was found to be pyroxene initiated by in situ Fe oxidation. Volumetric calculations suggest that this oxidation leads to the generation of porosity due to the formation of micro-fractures allowing for fluid transport and subsequent dissolution of biotite and plagioclase. The rapid ensuing plagioclase weathering leads to formation of high secondary porosity in the corestone over a distance of only a few cm and eventually to the final disaggregation of bedrock to saprolite. The first secondary phases are oxides or amorphous precipitates from which secondary minerals (mainly gibbsite, kaolinite and goethite) form. As oxidation is the first weathering reaction, the supply of O-2 is a rate-limiting factor for chemical weathering. Hence, the supply of O-2 and its consumption at depth connects processes at the weathering front with those at the Earth's surface in a feedback mechanism. The strength of the feedback depends on the relative weight of advective versus diffusive transport of O-2 through the weathering profile. The feedback will be stronger with dominating diffusive transport. The low weathering rate is explained by the nature of this feedback that is ultimately dependent on the transport of O-2 through the whole regolith, and on lithological factors such as low bedrock porosity and the amount of Fe-bearing primary minerals. Tectonic quiescence in this region and low pre-development erosion rate (attributed to a dense vegetation cover) minimize the rejuvenation of the thick and cohesive regolith column, finally leading to low denudation rates. (C) 2015 Elsevier B.V. All rights reserved.
The Misho complex in Northwest Iran is a prominent topographic massif bounded by well known active faults. Our new structural analysis of this area indicates that faulting has important role in the exhumation of this complex. The conjugate orientation of the North and South Misho Faults caused uplift in the Misho and exhumation of the Precambrian crystalline basement. Our structural and stratigraphic data shows that rapid uplift could have been initiation since the 21-22 Ma and exhumation rate was about 0.16 to 0.24 km/Ma. To refine this age, we performed U/Pb analysis of detrital zircon from the Upper Red Formation using LA-ICP-MS. We conducted AFT analysis on 6 basement samples from the hanging wall and 1 sample from the Upper Red Formation in the footwall NMF. Uplift in the hanging wall of NMF led to resting of sample 916 marl. This geochronologic and thermochronologic data shows that exhumation in the MC is diachronously along strike and affected by faults. The phase of exhumation is documented in the study area and entire Iranian plateau is related to the final closure of the Neo-Tethys and northward motion of the Arabian Plate.
There are far-reaching conceptual similarities between bi-static surface georadar and post-stack, "zero-offset" seismic reflection data, which is expressed in largely identical processing flows. One important difference is, however, that standard deconvolution algorithms routinely used to enhance the vertical resolution of seismic data are notoriously problematic or even detrimental to the overall signal quality when applied to surface georadar data. We have explored various options for alleviating this problem and have tested them on a geologically well-constrained surface georadar dataset. Standard stochastic and direct deterministic deconvolution approaches proved to be largely unsatisfactory. While least-squares-type deterministic deconvolution showed some promise, the inherent uncertainties involved in estimating the source wavelet introduced some artificial "ringiness". In contrast, we found spectral balancing approaches to be effective, practical and robust means for enhancing the vertical resolution of surface georadar data, particularly, but not exclusively, in the uppermost part of the georadar section, which is notoriously plagued by the interference of the direct air- and groundwaves. For the data considered in this study, it can be argued that band- limited spectral blueing may provide somewhat better results than standard band-limited spectral whitening, particularly in the uppermost part of the section affected by the interference of the air- and groundwaves. Interestingly, this finding is consistent with the fact that the amplitude spectrum resulting from least-squares-type deterministic deconvolution is characterized by a systematic enhancement of higher frequencies at the expense of lower frequencies and hence is blue rather than white. It is also consistent with increasing evidence that spectral "blueness" is a seemingly universal, albeit enigmatic, property of the distribution of reflection coefficients in the Earth. Our results therefore indicate that spectral balancing techniques in general and spectral blueing in particular represent simple, yet effective means of enhancing the vertical resolution of surface georadar data and, in many cases, could turn out to be a preferable alternative to standard deconvolution approaches.
Carbonate minerals are common in both marine and lacustrine records, and are frequently used for paleoenvironmental reconstructions. The sedimentary sequence of the endorheic Dead Sea and its precursors contain aragonite laminae that provide a detailed sedimentary archive of climatic, hydrologic, limnologic and environmental conditions since the Pleistocene. However, the interpretation of these archives requires a detailed understanding of the constraints and mechanisms affecting CaCO3 precipitation, which are still debated. The implications of aragonite precipitation in the Dead Sea and in its late Pleistocene predecessor (Lake Lisan) were investigated in this study by mixing natural and synthetic brines with a synthetic bicarbonate solution that mimics flash-floods composition, with and without the addition of extracellular polymeric substances (EPS). Aragonite precipitation was monitored, and precipitation rates and carbonate yields were calculated and are discussed with respect to modern aquatic environments. The experimental insights on aragonite precipitation are then integrated with microfacies analyses in order to reconstruct and constrain prevailing limnogeological processes and their hydroclimatic drivers under low (interglacial) and high (glacial) lake level stands. Aragonite precipitation took place within days to several weeks after the mixing of the brines with a synthetic bicarbonate solution. Incubation time was proportional to bicarbonate concentration, and precipitation rates were partially influenced by ionic strength. Additionally, extracellular polymeric substances inhibited aragonite precipitation for several months. As for the lake's water budget, our calculations suggest that the precipitation of a typical aragonite lamina (0.5 mm thick) during high lake stand requires unreasonable freshwater inflow from either surface or subsurface sources. This discrepancy can be resolved by considering one or a combination of the following scenarios; (1) discontinuous aragonite deposition over parts of the lake floor; (2) supply of additional carbonate flux (or fluxes) to the lake from aeolian dust and the remobilization and dissolution of dust deposits at the watershed; (3) carbonate production via oxidation of organic carbon by sulfate-reducing bacteria. Altogether, it is suggested that aragonite laminae thickness cannot be directly interpreted for quantitatively reconstructing the hydrological balance for the entire lake, they may still prove valuable for identifying inherent hydroclimatic periodicities at a single site.
In the comment on "Varves of the Dead Sea sedimentary record." Quaternary Science Reviews 215 (Ben Dor et al., 2019): 173-184. by R. Bookman, two recently published papers are suggested to prove that the interpretation of the laminated sedimentary sequence of the Dead Sea, deposited mostly during MIS2 and Holocene pluvials, as annual deposits (i.e., varves) is wrong. In the following response, we delineate several lines of evidence which coalesce to demonstrate that based on the vast majority of evidence, including some of the evidence provided in the comment itself, the interpretation of these sediments as varves is the more likely scientific conclusion. We further discuss the evidence brought up in the comment and its irrelevance and lack of robustness for addressing the question under discussion.
The sedimentary record of the Dead Sea provides an exceptional high-resolution archive of past climate changes in the drought-sensitive eastern Mediterranean-Levant, a key region for the development of humankind at the boundary of global climate belts. Moreover, it is the only deep hypersaline lake known to have deposited long sequences of finely laminated, annually deposited sediments (i.e. varves) of varied compositions, including aragonite, gypsum, halite and clastic sediments. Vast efforts have been made over the years to decipher the environmental information stored in these evaporitic-clastic sequences spanning from the Pleistocene Lake Amora to the Holocene Dead Sea. A general characterisation of sediment facies has been derived from exposed sediment sections, as well as from shallow- and deep-water sediment cores. During high lake stands and episodes of positive water budget, mostly during glacial times, alternating aragonite and detritus laminae (‘aad’ facies) were accumulated, whereas during low lake stands and droughts, prevailing during interglacials, laminated detritus (‘ld’ facies) and laminated halite (‘lh’ facies) dominate the sequence. In this paper, we (i) review the three types of laminated sediments of the Dead Sea sedimentary record (‘aad’, ‘ld’ and ‘lh’ facies), (ii) discuss their modes of formation, deposition and accumulation, and their interpretation as varves, and (iii) illustrate how Dead Sea varves are utilized for palaeoclimate reconstructions and for establishing floating chronologies.
Assessment of climate change impact on discharge of the lakhmass catchment (Northwest Tunisia)
(2022)
The Mediterranean region is increasingly recognized as a climate change hotspot but is highly underrepresented in hydrological climate change studies. This study aims to investigate the climate change effects on the hydrology of Lakhmass catchment in Tunisia. Lakhmass catchment is a part of the Medium Valley of Medjerda in northwestern Tunisia that drains an area of 126 km(2). First, the Hydrologiska Byrans Vattenbalansavdelning light (HBV-light) model was calibrated and validated successfully at a daily time step to simulate discharge during the 1981-1986 period. The Nash Sutcliffe Efficiency and Percent bias (NSE, PBIAS) were (0.80, +2.0%) and (0.53, -9.5%) for calibration (September 1982-August 1984) and validation (September 1984-August 1986) periods, respectively. Second, HBV-light model was considered as a predictive tool to simulate discharge in a baseline period (1981-2009) and future projections using data (precipitation and temperature) from thirteen combinations of General Circulation Models (GCMs) and Regional Climatic Models (RCMs). We used two trajectories of Representative Concentration Pathways, RCP4.5 and RCP8.5, suggested by the Intergovernmental Panel on Climate Change (IPCC). Each RCP is divided into three projection periods: near-term (2010-2039), mid-term (2040-2069) and long-term (2070-2099). For both scenarios, a decrease in precipitation and discharge will be expected with an increase in air temperature and a reduction in precipitation with almost 5% for every +1 degrees C of global warming. By long-term (2070-2099) projection period, results suggested an increase in temperature with about 2.7 degrees C and 4 degrees C, and a decrease in precipitation of approximately 7.5% and 15% under RCP4.5 and RCP8.5, respectively. This will likely result in a reduction of discharge of 12.5% and 36.6% under RCP4.5 and RCP8.5, respectively. This situation calls for early climate change adaptation measures under a participatory approach, including multiple stakeholders and water users.
Subduction zone magmas are more oxidised on eruption than those at mid-ocean ridges. This is attributed either to oxidising components, derived from subducted lithosphere (slab) and added to the mantle wedge, or to oxidation processes occurring during magma ascent via differentiation. Here we provide direct evidence for contributions of oxidising slab agents to melts trapped in the sub-arc mantle. Measurements of sulfur (S) valence state in sub-arc mantle peridotites identify sulfate, both as crystalline anhydrite (CaSO4) and dissolved SO42- in spinel-hosted glass (formerly melt) inclusions. Copper-rich sulfide precipitates in the inclusions and increased Fe3+/Sigma Fe in spinel record a S6+-Fe2+ redox coupling during melt percolation through the sub-arc mantle. Sulfate-rich glass inclusions exhibit high U/Th, Pb/Ce, Sr/Nd and delta S-34 (+ 7 to + 11%), indicating the involvement of dehydration products of serpentinised slab rocks in their parental melt sources. These observations provide a link between liberated slab components and oxidised arc magmas.
The rhodolithic slope deposits of a Burdigalian carbonate platform in Sardinia near Sedini were analyzed to reconstruct facies and palaeobathymetry. There is a distinct red-algal growth zonation along the platform slope. The clinoform rollover area consists of coralline-algal bindstones, which downslope change into a zone where rhodoliths are locally fused by progressive encrustation. Mid-slope rhodoliths are moderately branched, and downslope rhodoliths have fruticose protuberances, resulting in branching rhodolith growth patterns. There is a sharp change from the rhodolitic rudstones to the basinal, bivalve-dominated rudstones at the clinoform bottomsets. Red-algal genera identified include Sporolithon, Lithophyllum, Spongites, Hydrolithon, Mesophyllum, Lithoporella, Neogoniolithon, and other mastophoroids and melobesioids. Genera and subfamilies show a zonation along the clinoforms, allowing palaeobathymetric estimates. The clinoform rollovers formed at a water depth of around 40 m and the bottomsets around 60 m. Results from geometrical reconstruction show that coral reefs in the inner platform formed at water depths of around 20 m. Therefore, the Sedini carbonate platform is an example of a reef-bearing platform in which the edge or the platform-interior reefs do not build up to sea level.
The aim of this paper is to estimate the Hurst parameter of Fractional Gaussian Noise (FGN) using Bayesian inference. We propose an estimation technique that takes into account the full correlation structure of this process. Instead of using the integrated time series and then applying an estimator for its Hurst exponent, we propose to use the noise signal directly. As an application we analyze the time series of the Nile River, where we find a posterior distribution which is compatible with previous findings. In addition, our technique provides natural error bars for the Hurst exponent.
The displacement histories of the San Jacinto and southernmost San Andreas fault zones are constrained by offset data with ages in the range of 5 Ma to 5 ka. Apparent discrepancies between long- and short-term average displacement rates can be reconciled with a time-variable rate model. In this model, the displacement rate on the San Andreas decelerated from similar to35 mm/yr at 1.5 Ma to as low as 9 +/- 4 mm/yr by 90 ka. Over this same time period, the rate on the San Jacinto fault zone accelerated from an initial value of zero to a rate of 26 +/- 4 mm/yr. The data also imply that the rate of the San Andreas fault accelerated since ca. 90 ka, from similar to9 mm/yr to the modern rate of 27 +/- 4 mm/yr, whereas the San Jacinto decelerated from 26 +/- 4 mm/yr to the modern rate of 8 +/- 4 mm/yr. The time scale of these changes is significantly longer than the earthquake cycle, but shorter than time scales characteristic of lithospheric-scalle dynamics. The emergence of the San Jacinto fault zone ca. 1.5 Ma coincided with the development of a major restraining bend in the San Andreas fault zone, suggesting that the formation of new subparallell faults could be driven by conditions that inhibit displacement on preexisting faults
Analysis of past and present stimulation projects reveals that the temporal evolution and growth of maximum observed moment magnitudes may be linked directly to the injected fluid volume and hydraulic energy. Overall evolution of seismic moment seems independent of the tectonic stress regime and is most likely governed by reservoir specific parameters, such as the preexisting structural inventory. Data suggest that magnitudes can grow either in a stable way, indicating the constant propagation of self-arrested ruptures, or unbound, for which the maximum magnitude is only limited by the size of tectonic faults and fault connectivity. Transition between the two states may occur at any time during injection or not at all. Monitoring and traffic light systems used during stimulations need to account for the possibility of unstable rupture propagation from the very beginning of injection by observing the entire seismicity evolution in near-real time and at high resolution for an immediate reaction in injection strategy.
Plain Language Summary Predicting and controlling the size of earthquakes caused by fluid injection is currently the major concern of many projects associated with geothermal energy production. Here, we analyze the magnitude and seismic moment evolution with injection parameters for prominent geothermal and scientific projects to date. Evolution of seismicity seems to be largely independent of the tectonic stress background and seemingly depends on reservoir specific characteristics. We find that the maximum observed magnitudes relate linearly to the injected volume or hydraulic energy. A linear relation suggests stable growth of induced ruptures, as predicted by current models, or rupture growth may no longer depend on the stimulated volume but on tectonics. A system may change between the two states during the course of fluid injection. Close-by and high-resolution monitoring of seismic and hydraulic parameters in near-real time may help identify these fundamental changes in ample time to change injection strategy and manage maximum magnitudes.
Preparatory mechanisms accompanying or leading to nucleation of larger earthquakes have been observed at both laboratory and field scales, but conditions favoring the occurrence of observable preparatory processes are still largely unknown. In particular, it remains a matter of debate why some earthquakes occur spontaneously without noticeable precursors as opposed to events that are preceded by an extended failure process. In this study, we have generated new high-resolution seismicity catalogs framing the occurrence of 20 M-L > 2.5 earthquakes at The Geysers geothermal field in California. To this end, a seismicity catalog of the 11 days framing each large event was created. We selected 20 sequences sampling different hypocentral depths and hydraulic conditions within the field. Seismic activity and magnitude frequency distributions displayed by the different earthquake sequences are correlated with their location within the reservoir. Sequences located in the northwestern part of the reservoir show overall increased seismic activity and low b values, while the southeastern part is dominated by decreased seismic activity and higher b values. Periods of high injection coincide with high b values and vice versa. These observations potentially reflect varying differential and mean stresses and damage of the reservoir rocks across the field. About 50% of analyzed sequences exhibit no change in seismicity rate in response to the large main event. However, we find complex waveforms at the onset of the main earthquake, suggesting that small ruptures spontaneously grow into or trigger larger events.
Tectonic and metamorphic data for the Central Alps (Switzerland/Italy) are used to discuss this classic example of a Barrovian metamorphic terrain, notably the evolution of its thermal structure in space and time. Available P-T-t data indicate variable contributions of advective and conductive heat transport during collision and subsequent cooling and exhumation. Some areas experienced a prolonged period of partial melting while other areas, at the same time, show but moderate heating. The Barrow-type metamorphic field gradient observed in the final orogen is the result of two distinct tectonic processes, with their related advective and conductive heat transport processes. The two tectonic processes are (1) accretion of material within a subduction channel related to decompression and emplacement of high-pressure units in the middle crust and (2) wedging and related nappe formation in the continental lower plate. The second process postdates the first one. Wedging and underthrusting of continental lower plate material produces heat input into lower crustal levels, and this process is responsible for predominantly conductive heat transport in the overlying units. The interacting processes lead to different maximum temperatures at different times, producing the final Barrovian metamorphic field gradient. The south experienced rapid cooling, whereas the north shows moderate cooling rates. This discrepancy principally reflects differences in the temperature distribution in the deeper crust prior to cooling. Differences in the local thermal gradient that prevailed before the cooling also determined the relationships between cooling rate and exhumation rate in the different areas. Citation: Berger, A., S. M. Schmid, M. Engi, R. Bousquet, and M. Wiederkehr (2011), Mechanisms of mass and heat transport during Barrovian metamorphism: A discussion based on field evidence from the Central Alps (Switzerland/northern Italy), Tectonics, 30, TC1007, doi:10.1029/2009TC002622.
Tectonic and climatic control on evolution of rift lakes in the Central Kenya Rift, East Africa
(2009)
The long-term histories of the neighboring Nakuru-Elmenteita and Naivasha lake basins in the Central Kenya Rift illustrate the relative importance of tectonic versus climatic effects on rift-lake evolution and the formation of disparate sedimentary environments. Although modem climate conditions in the Central Kenya Rift are very similar for these basins, hydrology and hydrochemistry of present-day lakes Nakuru, Elmenteita and Naivasha contrast dramatically due to tectonically controlled differences in basin geometries, catchment size, and fluvial processes. In this study, we use eighteen C-14 and Ar-40/Ar-39 dated fluvio-lacustrine sedimentary sections to unravel the spatiotemporal evolution of the lake basins in response to tectonic and climatic influences. We reconstruct paleoclimatic and ecological trends recorded in these basins based on fossil diatom assemblages and geologic field mapping. Our study shows a tendency towards increasing alkalinity and shrinkage of water bodies in both lake basins during the last million years. Ongoing volcano-tectonic segmentation of the lake basins, as well as reorganization of upstream drainage networks have led to contrasting hydrologic regimes with adjacent alkaline and freshwater conditions. During extreme wet periods in the past, such as during the early Holocene climate optimum, lake levels were high and all basins evolved toward freshwater systems. During drier periods some of these lakes revert back to alkaline conditions, while others maintain freshwater characteristics. Our results have important implications for the use and interpretation of lake sediment as climate archives in tectonically active regions and emphasize the need to deconvolve lacustrine records with respect to tectonics versus climatic forcing mechanisms.
Three diatomite beds exposed in the Ol Njorowa Gorge south of Lake Naivasha, Central Kenya Rift, document three major lake-level highstands between 175 and 60 kyr BP. Diatom transfer-function estimates of hydrological and hydrochemical parameters suggest that a deep and large freshwater lake existed during the highstands at 135 and 80 kyr BP. In contrast, a shallower but more expanded freshwater lake existed at 110 kyr BP. The best analog for the most extreme highstand at 135 kyr BP is the highstand during the Early Holocene humid period from 10 to 6 kyr BP. The environmental conditions as reconstructed from diatom assemblages suggest long-lasting episodes of increased humidity during the high lake periods. This contrasts to the modern situation with a relatively shallow Lake Naivasha characterized by rapid water level fluctuations within a few decades. The most likely cause for the variable hydrological conditions since 175 kyr BP is orbitally driven insolation changes on the equator and increased lateral moisture transport from the ocean.
We modeled the two most extreme highstands of Lake Naivasha during the last 175 k.y. to estimate potential precipitation/ evaporation changes in this basin. In a first step, the bathymetry of the paleolakes at f135 and 9 k.y. BP was reconstructed from sediment cores and surface outcrops. Second, we modeled the paleohydrologic budget during the highstands using a simplified coupled energy mass-balance model. Our results show that the hydrologic and hence the climate conditions at f135 and 9 k.y. BP were similar, but significantly different from today. The main difference is a f15% higher value in precipitation compared to the present. An adaptation and migration of vegetation in the cause of climate changes would result in a f30% increase in precipitation. The most likely cause for such a wetter climate at f135 and 9 k.y. BP is a more intense intertropical convergence and increased precipitation in East Africa.
Three spinel lherzolite xenoliths from Mt. Quincan (Queensland, northeastern Australia) were studied with special attention to their enclosed fluid inclusions. The xenoliths are deformed, have porphyroclastic textures and overall show very similar petrographic features. The only significant difference is manifested in the abundance of fluid inclusions in the samples, mostly in orthopyroxene porphyroclasts. Xenolith JMTQ11 is fluid inclusion-free, whereas xenolith JMTQ20 shows a high abundance of fluid inclusions (fluid inclusion-rich). Xenolith JMTQ45 represents a transitional state between the previous two, as it contains only a small amount of fluid inclusions (fluid inclusion-bearing). Previous studies revealed that these xenoliths and the entrapped fluid inclusions represent a former addition of a MORB-type fluid to the pre-existing lithosphere, resulting from asthenosphere upwelling. There is a progressive enrichment in LREE, Nb, Sr and Ti from the fluid inclusion-free xenolith through the fluid inclusion-bearing one to the fluid inclusion-rich lherzolite. This suggests an increase in the extent of the interaction between the fluid-rich melt and the lherzolite wallrock. In addition, the same interaction is considered to be responsible for the formation of pargasitic amphibole as well. The presence of fluid inclusions indicates fluid migration at mantle depth, and their association with exsolution lamellae in orthopyroxene suggests fluid entrapment following the continental rifting (thermal relaxation) during cooling. A series of analyses, including microthermometry coupled with Raman spectroscopy, FTIR hyperspectral imaging, and Focused Ion Beam-Scanning Electron Microscopy (FIB-SEM) was carried out on the fluid inclusions. Based on the results, the entrapped high-density fluid is composed of 7589 mol% CO2, 918 mol% H2O, 0.11.7 mol% N-2 and <= 0.5 mol% H2S with dissolved trace elements (melt component). Our findings suggest that the metasomatic fluid phase could have been either a fluid/fluid-rich silicate melt released from the deeper asthenosphere, or a coexisting incipient fluid-rich silicate melt. Further cooling, possibly due to thermal relaxation and the upward migration of the fluid phase, caused the investigated lherzolites to reach pargasite stability conditions. We conclude that pargasite, even if only present in very limited modal proportions, can be a common phase at spinel lherzolite stability in the lithospheric upper mantle in continental rift back-arc settings. Studies of fluid inclusions indicate that significant CO2 release from the asthenosphere in a continental rifting environment is resulting from asthenosphere upwelling and its addition to the lithospheric mantle together with fluid-rich melt lherzolite interaction that leaves a CO2-rich fluid behind.
Scarce research about small European rivers from non-human impacted areas to determine their natural background state has been undertaken. During the annual hydrological cycle of 2008-9 the patterns of particulate supply (SPM, POC, PON, Al, Cd, Co, Cr, Cu, Fe, Ni, Pb, V, Zn) from the rivers Sor, Mera Landro, Lourido and Landoi to the Northern Galician Rias (SW Bay of Biscay) were tackled. No differences in the composition of the SPM were detected for the studied rivers regarding Al, Fe and POC but the relative percentage of particulate trace elements (PTE) discriminate the rivers. So, Cr, Co and Ni in the Lourido, and Landoi rivers, and Cu in the Mera River, are controlled by watershed minerals of Ortegal Geological Complex while for the rest rivers PTE are by granitic and Ollo de Sapo bedrock watershed. Therefore, the imprint of PTE in the parental rocks of the river basins is reflected on the coastal sediments of the Rias. The main process controlling the dynamics and variations of chemical elements in the particulate form is the river discharge. This fact exemplifies that these rivers presents a natural behavior not being highly influenced by anthropogenic activities.
During the last 5 Ma the Earth's ocean-atmosphere system passed through several major transitions, many of which are discussed as possible triggers for human evolution. A classic in this context is the possible influence of the closure of the Panama Strait, the intensification of Northern Hemisphere Glaciation, a stepwise increase in aridity in Africa, and the first appearance of the genus Homo about 2.5 - 2.7 Ma ago. Apart from the fact that the correlation between these events does not necessarily imply causality, many attempts to establish a relationship between climate and evolution fail due to the challenge of precisely localizing an a priori unknown number of changes potentially underlying complex climate records. The kernel-based Bayesian inference approach applied here allows inferring the location, generic shape, and temporal scale of multiple transitions in established records of Plio-Pleistocene African climate. By defining a transparent probabilistic analysis strategy, we are able to identify conjoint changes occurring across the investigated terrigenous dust records from Ocean Drilling Programme (ODP) sites in the Atlantic Ocean (ODP 659), Arabian (ODP 721/722) and Mediterranean Sea (ODP 967). The study indicates a two-step transition in the African climate proxy records at (2.35-2.10) Ma and (1.70 - 1.50) Ma, that may be associated with the reorganization of the Hadley-Walker Circulation. .
The effects of climate and topography on soil physico-chemical and microbial parameters were studied along an extensive latitudinal climate gradient in the Coastal Cordillera of Chile (26 degrees-38 degrees S). The study sites encompass arid (Pan de Azucar), semiarid (Santa Gracia), mediterranean (La Campana) and humid (Nahuelbuta) climates and vegetation, ranging from arid desert, dominated by biological soil crusts (biocrusts), semiarid shrubland and mediterranean sclerophyllous forest, where biocrusts are present but do have a seasonal pattern to temperate-mixed forest, where biocrusts only occur as an early pioneering development stage after disturbance. All soils originate from granitic parent materials and show very strong differences in pedogenesis intensity and soil depth. Most of the investigated physical, chemical and microbiological soil properties showed distinct trends along the climate gradient. Further, abrupt changes between the arid northernmost study site and the other semi-arid to humid sites can be shown, which indicate non-linearity and thresholds along the climate gradient. Clay and total organic carbon contents (TOC) as well as Ah horizons and solum depths increased from arid to humid climates, whereas bulk density (BD), pH values and base saturation (BS) decreased. These properties demonstrate the accumulation of organic matter, clay formation and element leaching as key-pedogenic processes with increasing humidity. However, the soils in the northern arid climate do not follow this overall latitudinal trend, because texture and BD are largely controlled by aeolian input of dust and sea salts spray followed by the formation of secondary evaporate minerals. Total soil DNA concentrations and TOC increased from arid to humid sites, while areal coverage by biocrusts exhibited an opposite trend. Relative bacterial and archaeal abundances were lower in the arid site, but for the other sites the local variability exceeds the variability along the climate gradient. Differences in soil properties between topographic positions were most pronounced at the study sites with the mediterranean and humid climate, whereas microbial abundances were independent on topography across all study sites. In general, the regional climate is the strongest controlling factor for pedogenesis and microbial parameters in soils developed from the same parent material. Topographic position along individual slopes of limited length augmented this effect only under humid conditions, where water erosion likely relocated particles and elements downward. The change from alkaline to neutral soil pH between the arid and the semi-arid site coincided with qualitative differences in soil formation as well as microbial habitats. This also reflects non-linear relationships of pedogenic and microbial processes in soils depending on climate with a sharp threshold between arid and semi-arid conditions. Therefore, the soils on the transition between arid and semi-arid conditions are especially sensitive and may be well used as indicators of long and medium-term climate changes. Concluding, the unique latitudinal precipitation gradient in the Coastal Cordillera of Chile is predestined to investigate the effects of the main soil forming factor - climate - on pedogenic processes.
Terrigenous sediment supply, marine transport, and depositional processes along tectonically active margins are key to decoding turbidite successions as potential archives of climatic and seismic forcings. Sequence stratigraphic models predict coarse-grained sediment delivery to deep-marine sites mainly during sea-level fall and lowstand. Marine siliciclastic deposition during transgressions and highstands has been attributed to sustained connectivity between terrigenous sources and marine sinks facilitated by narrow shelves. To decipher the controls on Holocene highstand turbidite deposition, we analyzed 12 sediment cores from spatially discrete, coeval turbidite systems along the Chile margin (29 degrees-40 degrees S) with changing climatic and geomorphic characteristics but uniform changes in sea level. Sediment cores from intraslope basins in north-central Chile (29 degrees-33 degrees S) offshore a narrow to absent shelf record a shut-off of turbidite deposition during the Holocene due to postglacial aridification. In contrast, core sites in south-central Chile (36 degrees-40 degrees S) offshore a wide shelf record frequent turbidite deposition during highstand conditions. Two core sites are linked to the Biobio river-canyon system and receive sediment directly from the river mouth. However, intraslope basins are not connected via canyons to fluvial systems but yield even higher turbidite frequencies. High sediment supply combined with a wide shelf and an undercurrent moving sediment toward the shelf edge appear to control Holocene turbidite sedimentation and distribution. Shelf undercurrents may play an important role in lateral sediment transport and supply to the deep sea and need to be accounted for in sediment-mass balances.
The details of how narrow, orogen-parallel ocean basins are filled with sediment by large axial submarine channels is important to understand because these depositional systems commonly form in through-like basins in various tectonic settings. The Magallanes foreland basin is an excellent location to study an orogen-parallel deep-marine system. Conglomerate lenses of the Upper Cretaceous Cerro Toro Formation have been previously interpreted to represent the fill of a single submarine channel (48 km wide, >100 km long) that funneled coarse detritus southward along the basin axis. This interpretation was based on lithologic correlations. New U/Pb dating of zircons from volcanic ashes and sandstones, coupled with strontium isotope stratigraphy, refine the controls on depositional ages and provenance. Results demonstrate that north-south oriented conglomerate lenses are contemporaneous within error limits (ca. 8482 Ma) supporting that they represent parts of an axial channel belt. Channel deposits 20 km west of the axial location are 8782 Ma in age. These channels are partly contemporaneous with the ones within the axial channel belt, making it likely that they represent feeders to the axial channel system. The northern Cerro Toro Formation spans a Turonian to Campanian interval (ca. 9082 Ma) whereas the formation top, 70 km to the south, is as young as ca. 76 Ma. KolmogorovSmirnoff statistical analysis on detrital zircon age distributions shows that the northern uppermost Cerro Toro Formation yields a statistically different age distribution than other samples from the same formation but shows no difference relative to the overlying Tres Pasos Formation. These results suggest the partly coeval deposition of both formations. Integration of previously acquired geochronologic and stratigraphic data with new data show a pronounced southward younging pattern in all four marine formations in the Magallanes Basin. Highly diachronous infilling may be an important depositional pattern for narrow, orogen-parallel ocean basins.
Newly acquired high-resolution bathymetric data (with 5 m and 2 m grid sizes) from the continental shelf off Concepcion (Chile), in combination with seismic reflection profiles, reveal a distinctly different evolution for the Biobio submarine canyon compared to that of one of its tributaries. Both canyons are incised into the shelf of the active margin. Whereas the inner shelf appears to be mantled with unconsolidated sediment, the outer shelf shows the influence of strong bottom currents that form drifts of loose sediment and transport -material into the Biobio submarine canyon and onto the continental slope.
The main stem of the Biobio Canyon is connected to the mouth of the Biobio River and currently provides a conduit for terrestrial sediment from the continental shelf to the deep seafloor. In contrast, the head of its tributary closest to the coast is located similar to 24 km offshore of the present-day coastline at 120 m water depth, and it is subject to passive sedimentation. However, canyon activity within the study area is interpreted to be controlled not only by the direct input of fluvial sediments into the canyon head facilitated by the river-mouth to canyon-head connection, but also by input from southward-directed bottom currents and possibly longshore drift. In addition, about 24 km offshore of the present-day coastline, the main stem of the Biobio Canyon has steep canyon walls next to sites of active tectonic deformation that are prone to wall failure. Mass-failure events may also foster turbidity currents and contribute to canyon feeding. In contrast, the tributary has less steep canyon walls with limited evidence of canyon-wall failure and is located down-system of bottom currents from the Biobio Canyon. It consequently receives neither fluvial nor longshore sediments. Therefore, the canyon's connectivity to fluvial or longshore sediment delivery pathways is affected by the distance of the canyon head from the coastline and the orientation of the canyon axis relative to the direction of bottom currents.
The ability of a submarine canyon to act as an active conduit for large quantities of terrestrial sediment toward the deep sea during sea-level highstands may be controlled by several different conditions simultaneously. These include bottom current direction, structural deformation of the seafloor affecting canyon location and orientation as well as canyon-wall failure, shelf gradient and associated distance from the canyon head to the coast, and fluvial networks. The complex interplay between these factors may vary even within an individual canyon system, resulting in distinct levels of canyon activity on a regional scale.
Understanding how Earth-surface processes respond to past climatic perturbations is crucial for making informed predictions about future impacts of climate change on sediment "uxes. Sedimentary records provide the archives for inferring these processes, but their interpretation is compromised by our incomplete understanding of how sediment-routing systems respond to millennial-scale climate cycles. We analyzed seven sediment cores recovered from marine turbidite depositional sites along the Chile continental margin. The sites span a pronounced arid-to-humid gradient with variable relief and related sediment connectivity of terrestrial and marine environments. These sites allowed us to study event related depositional processes in different climatic and geomorphic settings from the Last Glacial Maximum to the present day. The three sites reveal a steep decline of turbidite deposition during deglaciation. High rates of sea-level rise postdate the decline in turbidite deposition. Comparison with paleoclimate proxies documents that the spatio-temporal sedimentary pattern rather mirrors the deglacial humidity decrease and concomitant warming with no resolvable lag times. Our results let us infer that declining deglacial humidity decreased "uvial sediment supply. This signal propagated rapidly through the highly connected systems into the marine sink in north-central Chile. In contrast, in south-central Chile, connectivity between the Andean erosional zone and the "uvial transfer zone probably decreased abruptly by sediment trapping in piedmont lakes related to deglaciation, resulting in a sudden decrease of sediment supply to the ocean. Additionally, reduced moisture supply may have contributed to the rapid decline of turbidite deposition. These different causes result in similar depositional patterns in the marine sinks. We conclude that turbiditic strata may constitute reliable recorders of climate change across a wide range of climatic zones and geomorphic conditions. However, the underlying causes for similar signal manifestations in the sinks may differ, ranging from maintained high system connectivity to abrupt connectivity loss. (C) 2017 Elsevier B.V. All rights reserved.
Deposits of submarine debris flows can build up substantial topography on the sea floor. The resulting sea floor morphology can strongly influence the pathways of and deposition from subsequent turbidity currents. Map views of sea floor morphology are available for parts of the modern sea floor and from high-resolution seismic-reflection data. However, these data sets usually lack lithological information. In contrast, outcrops provide cross-sectional and lateral stratigraphic details of deep-water strata with superb lithological control but provide little information on sea floor morphology. Here, a methodology is presented that extracts fundamental lithological information from sediment core and well logs with a novel calibration between core, well-logs and seismic attributes within a large submarine axial channel belt in the Tertiary Molasse foreland basin, Austria. This channel belt was the course of multiple debris-flow and turbidity current events, and the fill consists of interbedded layers deposited by both of these processes. Using the core-well-seismic calibration, three-dimensional lithofacies proportion volumes were created. These volumes enable the interpretation of the three-dimensional distribution of the important lithofacies and thus the investigation of sea floor morphology produced by debris-flow events and its impact on succeeding turbidite deposition. These results show that the distribution of debris-flow deposits follows a relatively regular pattern of levees and lobes. When subsequent high-density turbidity currents encountered this mounded debris-flow topography, they slowed and deposited a portion of their sandy high-density loads just upstream of morphological highs. Understanding the depositional patterns of debris flows is key to understanding and predicting the location and character of associated sandstone accumulations. This detailed model of the filling style and the resulting stratigraphic architecture of a debris-flow dominated deep-marine depositional system can be used as an analogue for similar modern and ancient systems.
The Acheulean technological tradition, characterized by a large (>10 cm) flake-based component, represents a significant technological advance over the Oldowan. Although stone tool assemblages attributed to the Acheulean have been reported from as early as circa 1.6-1.75 Ma, the characteristics of these earliest occurrences and comparisons with later assemblages have not been reported in detail. Here, we provide a newly established chronometric calibration for the Acheulean assemblages of the Konso Formation, southern Ethiopia, which span the time period similar to 1.75 to <1.0 Ma. The earliest Konso Acheulean is chronologically indistinguishable from the assemblage recently published as the world's earliest with an age of similar to 1.75 Ma at Kokiselei, west of Lake Turkana, Kenya. This Konso assemblage is characterized by a combination of large picks and crude bifaces/unifaces made predominantly on large flake blanks. An increase in the number of flake scars was observed within the Konso Formation handaxe assemblages through time, but this was less so with picks. The Konso evidence suggests that both picks and handaxes were essential components of the Acheulean from its initial stages and that the two probably differed in function. The temporal refinement seen, especially in the handaxe forms at Konso, implies enhanced function through time, perhaps in processing carcasses with long and stable cutting edges. The documentation of the earliest Acheulean at similar to 1.75 Ma in both northern Kenya and southern Ethiopia suggests that behavioral novelties were being established in a regional scale at that time, paralleling the emergence of Homo erectus-like hominid morphology.
Constructing a hidden Markov Model based earthquake detector: application to induced seismicity
(2012)
The triggering or detection of seismic events out of a continuous seismic data stream is one of the key issues of an automatic or semi-automatic seismic monitoring system. In the case of dense networks, either local or global, most of the implemented trigger algorithms are based on a large number of active stations. However, in the case of only few available stations or small events, for example, like in monitoring volcanoes or hydrothermal power plants, common triggers often show high false alarms. In such cases detection algorithms are of interest, which show reasonable performance when operating even on a single station. In this context, we apply Hidden Markov Models (HMM) which are algorithms borrowed from speech recognition. However, many pitfalls need to be avoided to apply speech recognition technology directly to earthquake detection. We show the fit of the model parameters in an innovative way. State clustering is introduced to refine the intrinsically assumed time dependency of the HMMs and we explain the effect coda has on the recognition results. The methodology is then used for the detection of anthropogenicly induced earthquakes for which we demonstrate for a period of 3.9 months of continuous data that the single station HMM earthquake detector can achieve similar detection rates as a common trigger in combination with coincidence sums over two stations. To show the general applicability of state clustering we apply the proposed method also to earthquake classification at Mt. Merapi volcano, Indonesia.
The marine-controlled source electromagnetic method (mCSEM) is employed for studying the electrical characteristics and fluid contents of sedimentary reservoirs. However, the success rate of the method can be improved significantly by finding the sources of electromagnetic noise and addressing the challenge posed by them at larger offsets where the reservoir signal is often weak. I have studied the mCSEM data and reporting an electromagnetic noise. The strength of the noise is observed 1600 times stronger than the seafloor mCSEM signal at 0.1 Hz. Moreover, the noise and the transmitted mCSEM signals are found coherent in interstation recordings. These readings suggest the severity of the noise. The source investigation presuming the observed noise as an infragravity wave failed to match the response. Then, the role of microseisms is investigated. Microseism causes oscillation of the seafloor and produces electromagnetic disturbances by the dynamics of water. I have used various conditions for a proper discrimination of the noise as microseisms. This mechanism is clearly illustrated with the help of a conceptual diagram. The role of the directionality is part of the study, which is argued for having a significant role in the generation of microseisms. In this paper, a new algorithm is presented and is used for calculating the coherency. The algorithm helps in mapping the coherency value simultaneously in time and frequency domains.
A new solid-state material, N-butyl pyridinium diiodido argentate(I), is synthesized using a simple and effective one-pot approach. In the solid state, the compound exhibits 1D ([AgI2](-))(n) chains that are stabilized by the N-butyl pyridinium cation. The 1D structure is further manifested by the formation of long, needle-like crystals, as revealed from electron microscopy. As the general composition is derived from metal halide-based ionic liquids, the compound has a low melting point of 100-101 degrees C, as confirmed by differential scanning calorimetry. Most importantly, the compound has a conductivity of 10(-6) S cm(-1) at room temperature. At higher temperatures the conductivity increases and reaches to 10(-4 )S cm(-1) at 70 degrees C. In contrast to AgI, however, the current material has a highly anisotropic 1D arrangement of the ionic domains. This provides direct and tuneable access to fast and anisotropic ionic conduction. The material is thus a significant step forward beyond current ion conductors and a highly promising prototype for the rational design of highly conductive ionic solid-state conductors for battery or solar cell applications.
Ground‐motion prediction equations (GMPEs) are calibrated to predict the intensity of ground shaking at any given location, based on earthquake magnitude, source‐to‐site distance, local soil amplifications, and other parameters. GMPEs are generally assumed to be independent of time; however, evidence is increasing that large earthquakes modify the shallow soil conditions and those of the fault zone for months or years. These changes may affect the intensity of shaking and result in time‐dependent effects that can potentially be resolved by analyzing between‐event residuals (residuals between observed and predicted ground motion for individual earthquakes averaged over all stations). Here, we analyze a data set of about 65,000 recordings for about 1400 earthquakes in the moment magnitude range 2.5–6.5 that occurred in central Italy from 2008 to 2017 to capture the temporal variability of the ground shaking at high frequency. We first compute between‐event residuals for each earthquake in the Fourier domain with respect to a GMPE developed ad hoc for the analyzed data set. The between‐events show large changes after the occurrence of mainshocks such as the 2009 Mw 6.3 L'Aquila, the 2016 Mw 6.2 Amatrice, and Mw 6.5 Norcia earthquakes. Within the time span of a few months after the mainshocks, the between‐event contribution to the ground shaking varies by a factor 7. In particular, we find a large drop in the between‐events in the aftermath of the L'Aquila earthquake, followed by a slow positive trend that leads to a recovery interrupted by a new drop at the beginning of 2014. We also quantify the frequency‐dependent correlation between the Brune stress drop Δσ and the between‐events. We find that the temporal changes of Δσ resemble those of the between‐event residuals; in particular, during the period when the between‐events show the positive trend, the average logarithm of Δσ increases with an annual rate of 0.19 (i.e., the amplification factor for Δσ is 1.56 per year). Breakpoint analysis located a change in the linear trend coefficients of Δσ versus time in February 2014, although no large earthquakes occurred at that time. Finally, the temporal variability of Δσ mirrors the relative seismic‐velocity variations observed in previous studies for the same area and period, suggesting that both crack healing along the main fault system and healing of microcracks distributed at shallow depths throughout the surrounding region might be necessary to explain the wider observations of postearthquake recovery.
We present the results of a consistency check performed over the flatfile extracted from the engineering strong motion (ESM) database. The flatfile includes 23,014 recordings from 2179 earthquakes in the magnitude range from 3.5 to 7.8 that occurred since the 1970s in Europe and Middle East, as presented in the companion article by Lanzano et al. (Bull Earthq Eng, 2018a). The consistency check is developed by analyzing different residual distributions obtained from ad-hoc ground motion prediction equations for the absolute spectral acceleration (SA), displacement and Fourier amplitude spectra (FAS). Only recordings from earthquakes shallower than 40 km are considered in the analysis. The between-event, between-station and event-and-station corrected residuals are computed by applying a mixed-effect regression. We identified those earthquakes, stations, and recordings showing the largest deviations from the GMPE median predictions, and also evaluated the statistical uncertainty on the median model to get insights on the applicable magnitude–distance ranges and the usable period (or frequency) range. We observed that robust median predictions are obtained up to 8 s for SA and up to 20 Hz for FAS, although median predictions for Mw ≥ 7 show significantly larger uncertainties with ‘bumps’ starting above 5 s for SA and below 0.3 Hz for FAS. The between-station variance dominates over the other residual variances, and the dependence of the between-station residuals on logarithm of Vs30 is well-described by a piece-wise linear function with period-dependent slopes and hinge velocity around 580 m/s. Finally, we compared the between-event residuals obtained by considering two different sources of moment magnitude. The results show that, at long periods, the between-event terms from the two regressions have a weak correlation and the overall between-event variability is dissimilar, highlighting the importance of magnitude source in the regression results.
In this work, we analyse continuous measurements of microseisms to assess the reliability of the fundamental resonance frequency estimated by means of the horizontal-to-vertical (H/V) spectral ratio within the 0.1-1 Hz frequency range, using short-period sensors (natural period of 1 s). We apply the H/V technique to recordings of stations installed in two alluvial basins with different sedimentary cover thicknesses-the Lower Rhine Embayment (Germany) and the Gubbio Plain (Central Italy). The spectral ratios are estimated over the time-frequency domain, and we discuss the reliability of the results considering both the variability of the microseism activity and the amplitude of the instrumental noise. We show that microseisms measured by short period sensors allow the retrieval of fundamental resonance frequencies greater than about 0.1-0.2 Hz, with this lower frequency bound depending on the relative amplitude of the microseism signal and the self-noise of the instruments. In particular, we show an example where the considered short-period sensor is connected to instruments characterized by an instrumental noise level which allows detecting only fundamental frequencies greater than about 0.4 Hz. Since the frequency at which the peak of the H/V spectral ratio is biased depends upon the seismic signal-to-instrument noise ratio, the power spectral amplitude of instrumental self- noise should be always considered when interpreting the frequency of the peak as the fundamental resonance frequency of the investigated site.
We derive a set of regional ground-motion prediction equations (GMPEs) in the Fourier amplitude spectra (FAS-GMPE) and in the spectral acceleration (SA-GMPE) domains for the purpose of interpreting the between-event residuals in terms of source parameter variability. We analyze a dataset of about 65,000 recordings generated by 1400 earthquakes (moment magnitude 2: 5 <= M-w <= 6: 5, hypocentral distance R-hypo <= 150 km) that occurred in central Italy between January 2008 and October 2017. In a companion article (Bindi, Spallarossa, et al., 2018), the nonparametric acceleration source spectra were interpreted in terms of omega-square models modified to account for deviations from a high-frequency flat plateau through a parameter named k(source). Here, the GMPEs are derived considering the moment (M-w), the local (M-L), and the energy (M-e) magnitude scales, and the between-event residuals are computed as random effects. We show that the between-event residuals for the FAS-GMPE implementing M-w are correlated with stress drop, with correlation coefficients increasing with increasing frequency up to about 10 Hz. Contrariwise, the correlation is weak for the FAS-GMPEs implementing M-L and M-e, in particular between 2 and 5 Hz, where most of the corner frequencies lie. At higher frequencies, all models show a strong correlation with k(source). The correlation with the source parameters reflects in a different behavior of the standard deviation tau of the between-event residuals with frequency. Although tau is smaller for the FAS-GMPE using M-w below 1.5 Hz, at higher frequencies, the model implementing either M-L or M-e shows smaller values, with a reduction of about 30% at 3 Hz (i.e., from 0.3 for M-w to 0.1 for M-L). We conclude that considering magnitude scales informative for the stress-drop variability allows to reduce the between-event variability with a significant impact on the hazard assessment, in particular for studies in which the ergodic assumption on site is removed.
Rapidly changing climate in the Northern Hemisphere and associated socio-economic impacts require reliable understanding of lake systems as important freshwater resources and sensitive sentinels of environmental change. To better understand time-series data in lake sediment cores, it is necessary to gain information on within-lake spatial variabilities of environmental indicator data. Therefore, we retrieved a set of 38 samples from the sediment surface along spatial habitat gradients in the boreal, deep, and yet pristine Lake Bolshoe Toko in southern Yakutia, Russia. Our methods comprise laboratory analyses of the sediments for multiple proxy parameters, including diatom and chironomid taxonomy, oxygen isotopes from diatom silica, grain-size distributions, elemental compositions (XRF), organic carbon content, and mineralogy (XRD). We analysed the lake water for cations, anions, and isotopes. Our results show that the diatom assemblages are strongly influenced by water depth and dominated by planktonic species, i.e. Pliocaenicus bolshetokoensis. Species richness and diversity are higher in the northern part of the lake basin, associated with the availability of benthic, i.e. periphytic, niches in shallower waters. delta O-18(diatom) values are higher in the deeper south-western part of the lake, probably related to water temperature differences. The highest amount of the chironomid taxa underrepresented in the training set used for palaeoclimate inference was found close to the Utuk River and at southern littoral and profundal sites. Abiotic sediment components are not symmetrically distributed in the lake basin, but vary along restricted areas of differential environmental forcing. Grain size and organic matter are mainly controlled by both river input and water depth. Mineral (XRD) data distributions are influenced by the methamorphic lithology of the Stanovoy mountain range, while elements (XRF) are intermingled due to catchment and diagenetic differences. We conclude that the lake represents a valuable archive for multiproxy environmental reconstruction based on diatoms (including oxygen isotopes), chironomids, and sediment-geochemical parameters. Our analyses suggest multiple coring locations preferably at intermediate depth in the northern basin and the deep part in the central basin, to account for representative bioindicator distributions and higher temporal resolution, respectively.
Although the climate development over the Holocene in the Northern Hemisphere is well known, palaeolimnological climate reconstructions reveal spatiotemporal variability in northern Eurasia. Here we present a multi-proxy study from north-eastern Siberia combining sediment geochemistry, and diatom and pollen data from lake-sediment cores covering the last 38,000 cal. years. Our results show major changes in pyrite content and fragilarioid diatom species distributions, indicating prolonged seasonal lake-ice cover between similar to 13,500 and similar to 8900 cal. years BP and possibly during the 8200 cal. years BP cold event. A pollen-based climate reconstruction generated a mean July temperature of 17.8 degrees C during the Holocene Thermal Maximum (HTM) between similar to 8900 and similar to 4500 cal. years BP. Naviculoid diatoms appear in the late Holocene indicating a shortening of the seasonal ice cover that continues today. Our results reveal a strong correlation between the applied terrestrial and aquatic indicators and natural seasonal climate dynamics in the Holocene. Planktonic diatoms show a strong response to changes in the lake ecosystem due to recent climate warming in the Anthropocene. We assess other palaeolimnological studies to infer the spatiotemporal pattern of the HTM and affirm that the timing of its onset, a difference of up to 3000 years from north to south, can be well explained by climatic teleconnections. The westerlies brought cold air to this part of Siberia until the Laurentide ice sheet vanished 7000 years ago. The apparent delayed ending of the HTM in the central Siberian record can be ascribed to the exceedance of ecological thresholds trailing behind increases in winter temperatures and decreases in contrast in insolation between seasons during the mid to late Holocene as well as lacking differentiation between summer and winter trends in paleolimnological reconstructions. (C) 2015 Elsevier Ltd. All rights reserved.
Natural thermoluminescence (TL) in rocks reflects a dynamic equilibrium between radiation-induced TL growth and decay via thermal and athermal pathways. When rocks exhume through Earth's crust and cool from high to low temperature, this equilibrium level increases as the temperature dependent thermal decay decreases. This phenomenon can be exploited to extract thermal histories of rocks. The main advantage of TL is that a single TL glow curve has a wide range of thermal stabilities (lifetime <ka to Ba), and hence can provide multiple constraints on thermal histories. Here we constrain the distribution of kinetic parameters of TL in feldspar using a glow curve deconvolution method and fitting infinitesimal sub-peaks using a general order kinetic model. Each peak corresponds to a different but closely located energy level E. Forward modeling is applied for different time-temperature histories to estimate the sensitivity and limitation of each signal for different cooling rates. The results show that it is possible to constrain thermal histories between ∼30 °C and ∼80 °C. The results also illustrate that shallower traps, i.e. with lower activation energies, can be exploited to constrain lower cooling histories >100 °C/Ma, whereas deeper traps, i.e. with higher activation energies, provide constraints on thermal histories for higher cooling rates (>300 °C/Ma). Finally, we show how the path of rock exhumation (i.e., depth vs. time) can be constrained using an inverse approach. The newly developed methodology is applied to rapidly cooled samples from the Namche Barwa massif, eastern Himalaya to suggest a trend in exhumation rate with time that follows an inverse correlation with global temperature and glaciers equilibrium altitude line (ELA).
The concentration of sulfur that can be dissolved in a silicate liquid is of fundamental importance because it is closely associated with several major Earth-related processes. Considerable effort has been made to understand the interplay between the effects of silicate melt composition and its capac-ity to retain sulfur, but the dependence on pressure and temperature is mostly based on experiments performed at pressures and temperatures below 6 GPa and 2073 K. Here we present a study of the effects of pressure and temperature on sulfur content at sulfide saturation of a peridotitic liquid. We performed 14 multi-anvil experiments using a peridotitic starting composition, and we produced 25 new measurements at conditions ranging from 7 to 23 GPa and 2173 to 2623 K. We analyzed the recovered samples using both electron microprobe and laser ablation ICP-MS. We compiled our data together with previously published data that were obtained at lower P-T conditions and with various silicate melt compositions. We present a new model based on this combined data set that encompasses the entire range of upper mantle pressure-temperature conditions, along with the effect of a wide range of silicate melt compositions. Our findings are consistent with earlier work based on extrapolation from lower-pressure and lower-temperature experiments and show a decrease of sulfur content at sulfide saturation (SCSS) with increasing pressure and an increase of SCSS with increasing temperature. We have extrapolated our results to pressure-temperature conditions of the Earth's primitive magma ocean, and show that FeS will exsolve from the molten silicate and can effectively be extracted to the core by a process that has been termed the "Hadean Matte." We also discuss briefly the implications of our results for the lunar magma ocean.
High pressure and high temperature experiments performed with laser-heated diamond anvil cells (LH-DAC) are being extensively used in geosciences to study matter at conditions prevailing in planetary interiors. Due to the size of the apparatus itself, the samples that are produced are extremely small, on the order of few tens of micrometers. There are several ways to analyze the samples and extract physical, chemical or structural information, using either in situ or ex situ methods. In this paper, we compare two nanoprobe techniques, namely nano-XRF and NanoSIMS, that can be used to analyze recovered samples synthetized in a LH-DAC. With these techniques, it is possible to extract the spatial distribution of chemical elements in the samples. We show the results for several standards and discuss the importance of proper calibration for the acquisition of quantifiable results. We used these two nanoprobe techniques to retrieve elemental ratios of dilute species (few tens of ppm) in quenched experimental molten samples relevant for the formation of the iron-rich core of the Earth. We finally discuss the applications of such probes to constrain the partitioning of trace elements between metal and silicate phases, with a focus on moderately siderophile elements, tungsten and molybdenum.
Earthquake rupture length and width estimates are in demand in many seismological applications. Earthquake magnitude estimates are often available, whereas the geometrical extensions of the rupture fault mostly are lacking. Therefore, scaling relations are needed to derive length and width from magnitude. Most frequently used are the relationships of Wells and Coppersmith (1994) derived on the basis of a large dataset including all slip types with the exception of thrust faulting events in subduction environments. However, there are many applications dealing with earthquakes in subduction zones because of their high seismic and tsunamigenic potential. There are no well-established scaling relations for moment magnitude and length/width for subduction events. Within this study, we compiled a large database of source parameter estimates of 283 earthquakes. All focal mechanisms are represented, but special focus is set on (large) subduction zone events, in particular. Scaling relations were fitted with linear least-square as well as orthogonal regression and analyzed regarding the difference between continental and subduction zone/oceanic relationships. Additionally, the effect of technical progress in earthquake parameter estimation on scaling relations was tested as well as the influence of different fault mechanisms. For a given moment magnitude we found shorter but wider rupture areas of thrust events compared to Wells and Coppersmith (1994). The thrust event relationships for pure continental and pure subduction zone rupture areas were found to be almost identical. The scaling relations differ significantly for slip types. The exclusion of events prior to 1964 when the worldwide standard seismic network was established resulted in a remarkable effect on strike-slip scaling relations: the data do not show any saturation of rupture width of strike- slip earthquakes. Generally, rupture area seems to scale with mean slip independent of magnitude. The aspect ratio L/W, however, depends on moment and differs for each slip type.
Tsunami early warning (TEW) is a challenging task as a decision has to be made within few minutes on the basis of incomplete and error-prone data. Deterministic warning systems have difficulties in integrating and quantifying the intrinsic uncertainties. In contrast, probabilistic approaches provide a framework that handles uncertainties in a natural way. Recently, we have proposed a method using Bayesian networks (BNs) that takes into account the uncertainties of seismic source parameter estimates in TEW. In this follow-up study, the method is applied to 10 recent large earthquakes offshore Sumatra and tested for its performance. We have evaluated both the general model performance given the best knowledge we have today about the source parameters of the 10 events and the corresponding response on seismic source information evaluated in real-time. We find that the resulting site-specific warning level probabilities represent well the available tsunami wave measurements and observations. Difficulties occur in the real-time tsunami assessment if the moment magnitude estimate is severely over- or underestimated. In general, the probabilistic analysis reveals a considerably large range of uncertainties in the near-field TEW. By quantifying the uncertainties the BN analysis provides important additional information to a decision maker in a warning centre to deal with the complexity in TEW and to reason under uncertainty.
The northward indentation of the Pamir salient into the Tarim basin at the western syntaxis of the India-Asia collision zone is the focus of controversial models linking lithospheric to surface and atmospheric processes. Here we report on tectonic events recorded in the most complete and best-dated sedimentary sequences from the western Tarim basin flanking the eastern Pamir (the Aertashi section), based on sedimentologic, provenance, and magnetostratigraphic analyses. Increased tectonic subsidence and a shift from marine to continental fluvio-deltaic deposition at 41Ma indicate that far-field deformation from the south started to affect the Tarim region. A sediment accumulation hiatus from 24.3 to 21.6Ma followed by deposition of proximal conglomerates is linked to fault propagation into the Tarim basin. From 21.6 to 15.0Ma, increasing accumulation rates of fining upward clastics is interpreted as the expression of a major dextral transtensional system linking the Kunlun to the Tian Shan ahead of the northward Pamir indentation. At 15.0Ma, the appearance of North Pamir-sourced conglomerates followed at 11Ma by Central Pamir-sourced volcanics coincides with a shift to E-W compression, clockwise vertical-axis rotations and the onset of growth strata associated with the activation of the local east vergent Qimugen thrust wedge. Together, this enables us to interpret that Pamir indentation into Tarim had started by 24.3Ma, reached the study location by 15.0Ma and had passed it by 11Ma, providing kinematic constraints on proposed tectonic models involving intracontinental subduction and delamination.
The Pamirs represent the indented westward continuation of the northern margin of the Tibetan Plateau, dividing the Tarim and Tajik basins. Their evolution may be a key factor influencing aridification of the Asian interior, yet the tectonics of the Pamir Salient are poorly understood. We present a provenance study of the Aertashi section, a Paleogene to late Neogene clastic succession deposited in the Tarim basin to the north of the NW margin of Tibet (the West Kunlun) and to the east of the Pamirs. Our detrital zircon U-Pb ages coupled with zircon fission track, bulk rock Sm-Nd, and petrography data document changes in contributing source terranes during the Oligocene to Miocene, which can be correlated to regional tectonics. We propose a model for the evolution of the Pamir and West Kunlun (WKL), in which the WKL formed topography since at least similar to 200 Ma. By similar to 25 Ma, movement along the Pamir-bounding faults such as the Kashgar-Yecheng Transfer System had commenced, marking the onset of Pamir indentation into the Tarim-Tajik basin. This is coincident with basinward expansion of the northern WKL margin, which changed the palaeodrainage pattern within the Kunlun, progressively cutting off the more southerly WKL sources from the Tarim basin. An abrupt change in the provenance and facies of sediments at Aertashi has a maximum age of 14 Ma; this change records when the Pamir indenter had propagated sufficiently far north that the North Pamir was now located proximal to the Aertashi region.
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