Refine
Has Fulltext
- no (214) (remove)
Year of publication
- 2014 (214) (remove)
Document Type
- Article (196)
- Preprint (8)
- Doctoral Thesis (6)
- Review (4)
Is part of the Bibliography
- yes (214)
Keywords
- Earthquake source observations (7)
- Holocene (6)
- Tibetan Plateau (4)
- Erosion (3)
- Himalaya (3)
- Seismicity and tectonics (3)
- Aleatory variability (2)
- Andes (2)
- Baseline shift (2)
- Body waves (2)
Institute
- Institut für Geowissenschaften (214) (remove)
The complementary advantages of high-rate Global Positioning System (GPS) and accelerometer observations for measuring seismic ground motion have been recognised in previous research. Here we propose an approach of tight integration of GPS and accelerometer measurements. The baseline shifts of the accelerometer are introduced as unknown parameters and estimated by a random walk process in the Precise Point Positioning (PPP) solution. To demonstrate the performance of the new strategy, we carried out several experiments using collocated GPS and accelerometer. The experimental results show that the baseline shifts of the accelerometer are automatically corrected, and high precision coseismic information of strong ground motion can be obtained in real-time. Additionally, the convergence and precision of the PPP is improved by the combined solution.
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.
The basin-scale spatial variability in lipid biomarker proxies in lacustrine sediments, which are established tools for studying continental environmental change, has rarely been examined. It is often implicitly assumed that a lake sediment core provides an average integral of catchment sources. Here we evaluated the distribution of lipid biomarkers in a modern ecosystem and compared it with the sedimentary record. We analyzed lipid biomarkers in terrestrial and aquatic organisms and in lake surface sediments from 17 locations within the saline-alkaline Lonar crater lake in central India. Terrestrial vegetation and lake surface sediments were characterized by relatively high average chain length (ACL) index values (29.6-32.8) of leaf wax n-alkanes, consistent with suggestions that plants in drier and warmer climates produce longer chain alkyl lipids than plants in cooler and humid areas. A heterogeneous spatial distribution of ACL values in lake surface sediments was found: at locations away from the shore, the values were highest (31 or more), possibly indicating different sources and/or transport of terrestrial biomarkers. In floating, benthic microbial mats and surface sediment, n-heptadecane, carotenoids, diploptene, phytol and tetrahymanol occurred in large amounts. Interestingly, these biomarkers of a unique bacterial community were found in substantially higher concentrations in nearshore sediment samples. We suggest that human influence and subsequent nutrient supply resulted in increased primary productivity, leading to an unusually high concentration of tetrahymanol in the nearshore sediments. In summary, the data showed that substantial heterogeneity existed within the lake, but leaf wax n-alkanes in a core from the center of the lake represented an integral of catchment conditions. However, lake level fluctuation may potentially affect aquatic lipid biomarker distributions in lacustrine sediments, in addition to source changes.
The complementary advantages of GPS and seismic measurements are well recognized in seismotectonic monitoring studies. Therefore, integrated processing of the two data streams has been proposed recently in an attempt to obtain accurate and reliable information of surface displacements associated with earthquakes. A hitherto still critical issue in the integrated processing is real-time detection and precise estimation of the transient baseline error in the seismic records. Here, we report on a new approach by introducing the seismic acceleration corrected by baseline errors into the state equation system. The correction is performed and regularly updated in short epochs (with increments which may be as short as seconds), so that station position, velocity, and acceleration can be constrained very tightly and baseline error can be estimated as a random-walk process. With the adapted state equation system, our study highlights the use of a new approach developed for integrated processing of GPS and seismic data by means of sequential least-squares adjustment. The efficiency of our approach is demonstrated and validated using simulated, experimental, and real datasets. The latter were collected at collocated GPS and seismic stations around the 4 April 2010, E1 Mayor-Cucapah earthquake (Mw, 7.2). The results have shown that baseline errors of the strong-motion sensors are corrected precisely and high-precision seismic displacements are real-timely obtained by the new approach.
We focus on the relation between seismic and total postseismic afterslip following the Maule M-w 8.8 earthquake on 2010 February 27 in central Chile. First, we calculate the cumulative slip released by aftershock seismicity. We do this by summing up the aftershock regions and slip estimated from scaling relations. Comparing the cumulative seismic slip with afterslip modelswe showthat seismic slip of individual aftershocks exceeds locally the inverted afterslip model from geodetic constraints. As the afterslip model implicitly contains the displacements from the aftershocks, this reflects the tendency of afterslip models to smear out the actual slip pattern. However, it also suggests that locally slip for a number of the larger aftershocks exceeds the aseismic slip in spite of the fact that the total equivalent moment of the afterslip exceeds the cumulative moment of aftershocks by a large factor. This effect, seen weakly for the Maule 2010 and also for the Tohoku 2011 earthquake, can be explained by taking into account the uncertainties of the seismicity and afterslip models. In spite of uncertainties, the hypocentral region of the Nias 2005 earthquake is suggested to release a large fraction of moment almost purely seismically. Therefore, these aftershocks are not driven solely by the afterslip but instead their slip areas have probably been stressed by interseismic loading and the mainshock rupture. In a second step, we divide the megathrust of the Maule 2010 rupture into discrete cells and count the number of aftershocks that occur within 50 km of the centre of each cell as a function of time. We then compare this number to a time-dependent afterslip model by defining the 'afterslip to aftershock ratio' (ASAR) for each cell as the slope of the best fitting line when the afterslip at time t is plotted against aftershock count. Although we find a linear relation between afterslip and aftershocks for most cells, there is significant variability in ASAR in both the downdip and along-strike directions of the megathrust. We compare the spatial distribution of ASAR with the spatial distribution of seismic coupling, coseismic slip and Bouguer gravity anomaly, and in each case we find no significant correlation.
The Mohorovicic discontinuity, Moho for short, which marks the boundary between crust and mantle, is the main first-order structure within the lithosphere. Geodynamics and tectonic evolution determine its depth level and properties. Here, we present a map of the Moho in central Europe across the Teisseyre-Tornquist Zone, a region for which a number of previous studies are available. Our results are based on homogeneous and consistent processing of P-and S-receiver functions for the largest passive seismological data set in this region yet, consisting of more than 40 000 receiver functions from almost 500 station. Besides, we also provide new results for the crustal vP/vS ratio for the whole area.
Our results are in good agreement with previous, more localized receiver function studies, as well as with the interpretation of seismic profiles, while at the same time resolving a higher level of detail than previous maps covering the area, for example regarding the Eifel Plume region, Rhine Graben and northern Alps. The close correspondence with the seismic data regarding crustal structure also increases confidence in use of the data in crustal corrections and the imaging of deeper structure, for which no independent seismic information is available.
In addition to the pronounced, stepwise transition from crustal thicknesses of 30 km in Phanerozoic Europe to more than 45 beneath the East European Craton, we can distinguish other terrane boundaries based on Moho depth as well as average crustal v(P)/v(S) ratio and Moho phase amplitudes. The terranes with distinct crustal properties span a wide range of ages, from Palaeoproterozoic in Lithuania to Cenozoic in the Alps, reflecting the complex tectonic history of Europe. Crustal thickness and properties in the study area are also markedly influenced by tectonic overprinting, for example the formation of the Central European Basin System, and the European Cenozoic Rift System. In the areas affected by Cenozoic rifting and volcanism, thinning of the crust corresponds to lithospheric updoming reported in recent surface wave and S-receiver function studies, as expected for thermally induced deformation. The same correlation applies for crustal thickening, not only across the Trans-European Suture Zone, but also within the southern part of the Bohemian Massif.
We examine the use of ambient noise cross-correlation tomography for shallow site characterization using a modified two-step approach. Initially, we extract Rayleigh wave traveltimes from correlation traces of vertical component seismic recordings from a local network installed in Mygdonia basin, northern Greece. The obtained Rayleigh wave traveltimes show significant spatial variability, as well as distance and frequency dependence due to the 3-D structure of the area, dispersion, and anelastic attenuation effects. The traveltime data sets are inverted through a surface wave tomography approach to determine group velocity maps for each frequency. The proposed tomographic inversion involves the use of approximate Fresnel volumes and interfrequency smoothing constraints to stabilize the results. In the last step, we determine a final 3-D velocity model using a node-based Monte Carlo 1-D dispersion curve inversion. The reliability of the final 3-D velocity model is examined by spatial and depth resolution analysis, as well as by inversion for different model parameterizations. The obtained results are in very good agreement with previous findings from seismic and other geophysical methods. The new 3-D VS model provides additional structural constraints for the shallow sediments and bedrock structure of the northern Mygdonia basin up to the depth of similar to 200-250 m. Present work results suggest that the migration of ambient tomography techniques from large scales (tens or hundreds of km) to local scales (few hundred meters) is possible but cannot be used as a black box technique for 3-D modeling and detailed geotechnical site characterization.
In this study, we suggest a novel approach for the retrieval of regional moment tensors for earthquakes with small to moderate magnitudes. The first modification is the combined inversion of broadband and short-period waveform data. The broadband waveforms are inverted in a frequency range suitable for surface waves, whereas for the short-period data a frequency range suitable for body waves is applied. The second modification is the use of first-motion body-wave polarities to select the most probable solution out of all solutions from inversion. To combine three different criteria for selecting the most probable solution (i.e., residual from inversion, double-couple content of solution, number of nonmatching first-motion body-wave polarities), the L2 norm is applied to the normalized parameters. We chose five earthquakes within the Alborz mountains, Iran, as a case study (3.1 <= M-w <= 4.1). In this area, several factors exacerbate the difficulty of performing inversion for moment tensors, for example, a heterogeneous station network and large azimuthal gaps. We have demonstrated that our approach supplies reliable moment tensors when inversion from broadband data alone fails. In one case, we successfully retrieved a stable solution from short-period waveform data alone. Thus, our approach enables successful determination of seismic moment tensors wherever a sparse network of broadband stations has thus far prevented it.
The Ceres earthquake of 29 September 1969 is the largest known earthquake in southern Africa. Digitized analog recordings from Worldwide Standardized Seismographic Network stations (Powell and Fries, 1964) are used to retrieve the point source moment tensor and the most likely centroid depth of the event using full waveform modeling. A scalar seismic moment of 2.2-2.4 x 10(18) N center dot m corresponding to a moment magnitude of 6.2-6.3 is found. The analysis confirms the pure strike-slip mechanism previously determined from onset polarities by Green and Bloch (1971). Overall good agreement with the fault orientation previously estimated from local aftershock recordings is found. The centroid depth can be constrained to be less than 15 km. In a second analysis step, we use a higher order moment tensor based inversion scheme for simple extended rupture models to constrain the lateral fault dimensions. We find rupture propagated unilaterally for 4.7 s from east-southwest to west-northwest for about 17 km ( average rupture velocity of about 3: 1 km/s).
We have redetermined focal depths of moderate and major earthquakes with reported lower-crust and upper-mantle depths that have occurred in Tien-Shan, since the availability of broad-band array data. Records of earthquakes at global arrays have been used for identification and modelling of depth phases in order to make accurate estimation of focal depths. Our results show that half of the purportedly deep earthquakes are indeed originating from depths attributable to middle-crust and lower-crust regions. Also one exceptional event in the northern foreland of Tien-Shan in Junggar Basin is located in the upper mantle at the depth of 64 km. Such unusually deep earthquakes for intraplate continental tectonic domain are all located at the margin of Tien-Shan with its adjacent stable blocks and at least some of them have occurred where the brittle behaviour of continental rocks is not highly expected. The reverse mechanisms of all these earthquakes and their proximity to formerly subducting and later colliding and underplating stable blocks and their interactions with overlying Tien-Shan are clues to explain this extremity.
We track a bilateral rupture propagation lasting similar to 160 s, with its dominant branch rupturing northeastwards at about 3 kms(-1). The area of maximum energy emission is offset from the maximum coseismic slip but matches the zone where most plate interface aftershocks occur. Along dip, energy is preferentially released from two disconnected interface belts, and a distinct jump from the shallower belt to the deeper one is visible after about 20 s from the onset. However, both belts keep on being active until the end of the rupture. These belts approximately match the position of the interface aftershocks, which are split into two clusters of events at different depths, thus suggesting the existence of a repeated transition from stick-slip to creeping frictional regime.
Soils in various places of the Panama Canal Watershed feature a low saturated hydraulic conductivity (K-s) at shallow depth, which promotes overland-flow generation and associated flashy catchment responses. In undisturbed forests of these areas, overland flow is concentrated in flow lines that extend the channel network and provide hydrological connectivity between hillslopes and streams. To understand the dynamics of overland-flow connectivity, as well as the impact of connectivity on catchment response, we studied an undisturbed headwater catchment by monitoring overland-flow occurrence in all flow lines and discharge, suspended sediment, and total phosphorus at the catchment outlet. We find that connectivity is strongly influenced by seasonal variation in antecedent wetness and can develop even under light rainfall conditions. Connectivity increased rapidly as rainfall frequency increased, eventually leading to full connectivity and surficial drainage of entire hillslopes. Connectivity was nonlinearly related to catchment response. However, additional information on factors such as overland-flow volume would be required to constrain relationships between connectivity, stormflow, and the export of suspended sediment and phosphorus. The effort to monitor those factors would be substantial, so we advocate applying the established links between rain event characteristics, drainage network expansion by flow lines, and catchment response for predictive modeling and catchment classification in forests of the Panama Canal Watershed and in similar regions elsewhere.
Sedimentation in the floodplains of the Mekong Delta, Vietnam. Part I: suspended sediment dynamics
(2014)
Suspended sediment is the primary source for a sustainable agro-ecosystem in the Mekong Delta by providing nutrient input for the subsequent cropping season. In addition, the suspended sediment concentration (SSC) plays an important role in the erosion and deposition processes in the Delta; that is, it influences the morphologic development and may counteract the deltaic subsidence and sea level rise. Despite this importance, little is known about the dynamics of suspended sediment in the floodplains of the Mekong Delta. In particular, quantitative analyses are lacking mainly because of data scarcity with respect to the inundation processes in the floodplains. In 2008, therefore, a comprehensive in situ system to monitor the dynamics of suspended sediment in a study area located in the Plain of Reeds was established, aiming at the characterization and quantification of suspended sediment dynamics in the deeply inundated parts of the Vietnamese part of the Mekong Delta. The monitoring system was equipped with seven water quality-monitoring stations. They have a robust design and autonomous power supply suitable for operation on inundated floodplains, enabling the collection of reliable data over a long period of time with a high temporal resolution. The data analysis shows that the general seasonal dynamics of suspended sediment transport in the Delta is controlled by two main mechanisms: the flood wave of the Mekong River and the tidal backwater influences from the coast. In the channel network, SSC decreases exponentially with distance from the Mekong River. The anthropogenic influence on SSC could also be identified for two periods: at the start of the floodplain inundation and at the end of the flood period, when subsequent paddy rice crops are prepared. Based on the results, we recommend an operation scheme for the sluice gates, which intends to distribute the sediment and thus the nutrients equally over the floodplain.
Whilst sophisticated multiphase fluid flow models are routinely employed to understand behaviour of oil and gas reservoirs, high-resolution data describing the three-dimensional (3D) distribution of rock characteristics is rarely available to populate models. We present a new approach to developing a quantitative understanding of the effect of individual controls on the distribution of petrophysical properties and their impact on fluid flow. This involves simulating flow through high-detail permeability architectures generated by forward modelling of the coupled depositional-diagenetic evolution of isolated platforms using CARB3D(+). This workflow is exemplified by an investigation of interactions between subsidence and climate, and their expression in spatial variations in reservoir quality in an isolated carbonate platform of similar size and subsidence history to the Triassic Latemar Platform.
Dissolutional lowering during subaerial exposure controls platform-top graininess via platform top hydrodynamics during the subsequent transgression. Dissolved carbonate is reprecipitated as cements by percolating meteoric waters. However, associated subsurface meteoric dissolution generates significant secondary porosity under a more humid climate. Slower subsidence enhances diagenetic overprinting during repeated exposure events. Single-phase streamline simulations show how early diagenesis develops more permeable fairways within the finer-grained condensed units that can act as thief zones for flow from the grainier but less diagenetically altered cyclic units.
Numerical simulation of fluid-flow processes in a 3D high-resolution carbonate reservoir analogue
(2014)
A high-resolution three-dimensional (3D) outcrop model of a Jurassic carbonate ramp was used in order to perform a series of detailed and systematic flow simulations. The aim of this study was to test the impact of small- and large-scale geological features on reservoir performance and oil recovery. The digital outcrop model contains a wide range of sedimentological, diagenetic and structural features, including discontinuity surfaces, shoal bodies, mud mounds, oyster bioherms and fractures. Flow simulations are performed for numerical well testing and secondary oil recovery. Numerical well testing enables synthetic but systematic pressure responses to be generated for different geological features observed in the outcrops. This allows us to assess and rank the relative impact of specific geological features on reservoir performance. The outcome documents that, owing to the realistic representation of matrix heterogeneity, most diagenetic and structural features cannot be linked to a unique pressure signature. Instead, reservoir performance is controlled by subseismic faults and oyster bioherms acting as thief zones. Numerical simulations of secondary recovery processes reveal strong channelling of fluid flow into high-permeability layers as the primary control for oil recovery. However, appropriate reservoir-engineering solutions, such as optimizing well placement and injection fluid, can reduce channelling and increase oil recovery.
Carbonate reservoirs pose significant challenges for reservoir modelling and flow prediction due to heterogeneities in rock properties, limits to seismic resolution and limited constraints on subsurface data. Hence, a systematic and streamlined approach is needed to construct geological models and to quickly evaluate key sensitivities in the flow models. This paper discusses results from a reservoir analogue study of a Middle Jurassic carbonate ramp in the High Atlas Mountains of Morocco that has stratigraphic and structural similarities to selected Middle East reservoirs. For this purpose, high-resolution geological models were constructed from the integration of sedimentological, diagenetic and structural studies in the area. The models are approximately 1200 x 1250 m in size, and only faults (no fractures) with offsets greater than 1 m are included. Novel methods have been applied to test the response of flow simulations to the presence or absence of specific geological features, including proxies for hardgrounds, stylolites, patch reefs, and mollusc banks, as a way to guide the level of detail that is suitable for modelling objectives. Our general conclusion from the study is that the continuity of any geological feature with extreme permeability (high or low) has the most significant impact on flow.
In situ UV laser spot Ar-40/Ar-39 analyses of distinct phengite types in eclogite-facies rocks from the Sesia-Lanzo Zone (Western Alps, Italy) were combined with SIMS boron isotope analyses as well as boron (B) and lithium (Li) concentration data to link geochronological information with constraints on fluid-rock interaction. In weakly deformed samples, apparent Ar-40/Ar-39 ages of phengite cores span a range of similar to 20 Ma, but inverse isochrons define two distinct main high-pressure (HP) phengite core crystallization periods of 88-82 and 77-74 Ma, respectively. The younger cores have on average lower B contents (similar to 36 mu g/g) than the older ones (similar to 43-48 mu g/g), suggesting that loss of B and resetting of the Ar isotopic system were related. Phengite cores have variable delta B-11 values (-18 parts per thousand to -10 parts per thousand), indicating the lack of km scale B homogenization during HP crystallization.
Overprinted phengite rims in the weakly deformed samples generally yield younger apparent Ar-40/Ar-39 ages than the respective cores. They also show variable effects of heterogeneous excess 40 Ar incorporation and Ar loss. One acceptable inverse isochron age of 77.1 +/- 1.1 Ma for rims surrounding older cores (82.6 +/- 0.6 Ma) overlaps with the second period of core crystallization. Compared to the phengite cores, all rims have lower B and Li abundances but similar delta B-11 values (-15 parts per thousand to -9 parts per thousand), reflecting internal redistribution of B and Li and internal fluid buffering of the B isotopic composition during rim growth. The combined observation of younger Ar-40/Ar-39 ages and boron loss, yielding comparable values of both parameters only in cores and rims of different samples, is best explained by a selective metasomatic overprint. In low permeability samples, this overprint caused recrystallization of phengite rims, whereas higher permeability in other samples led to complete recrystallization of phengite grains.
Strongly deformed samples from a several km long, blueschist-facies shear zone contain mylonitic phengite that forms a tightly clustered group of relatively young apparent Ar-40/Ar-39 ages (64.7-68.8 Ma), yielding an inverse isochron age of 65.0 +/- 3.0 Ma. Almost complete B and Li removal in mylonitic phengite is due to leaching into a fluid. The B isotopic composition is significantly heavier than in phengites from the weakly deformed samples, indicating an external control by a high-delta B-11 fluid (delta B-11 = + 7 +/- 4 parts per thousand). We interpret this result as reflecting phengite recrystallization related to deformation and associated fluid flow in the shear zone. This event also caused partial resetting of the Ar isotope system and further B loss in more permeable rocks of the adjacent unit. We conclude that geochemical evidence for pervasive or limited fluid flow is crucial for the interpretation of Ar-40/Ar-39 data in partially metasomatized rocks.
The Pamir is the western continuation of Tibet and the site of some of the highest mountains on Earth, yet comparatively little is known about its crustal and tectonic evolution and erosional history. Both Tibet and the Pamir are characterized by similar terranes and sutures that can be correlated along strike, although the details of such correlations remain controversial. The erosional history of the Pamir with respect to Tibet is significantly different as well: Most of Tibet has been characterized by internal drainage and low erosion rates since the early Cenozoic; in contrast, the Pamir is externally drained and topographically more rugged, and it has a strongly asymmetric drainage pattern. Here, we report 700 new U-Pb and Lu-Hf isotope determinations and >300 Ar-40/Ar-39 ages from detrital minerals derived from rivers in China draining the northeastern Pamir and >1000 apatite fission-track (AFT) ages from 12 rivers in Tajikistan and China draining the northeastern, central, and southern Pamir. U-Pb ages from rivers draining the northeastern Pamir are Mesozoic to Proterozoic and show affinity with the Songpan-Ganzi terrane of northern Tibet, whereas rivers draining the central and southern Pamir are mainly Mesozoic and show some affinity with the Qiangtang terrane of central Tibet. The epsilon(Hf) values are juvenile, between 15 and -5, for the northeastern Pamir and juvenile to moderately evolved, between 10 and -40, for the central and southern Pamir. Detrital mica Ar-40/Ar-39 ages for the northeastern Pamir (eastern drainages) are generally older than ages from the central and southern Pamir (western drainages), indicating younger or lower-magnitude exhumation of the northeastern Pamir compared to the central and southern Pamir. AFT data show strong Miocene-Pliocene signals at the orogen scale, indicating rapid erosion at the regional scale. Despite localized exhumation of the Mustagh-Ata and Kongur-Shan domes, average erosion rates for the northeastern Pamir are up to one order of magnitude lower than erosion rates recorded by the central and southern Pamir. Deeper exhumation of the central and southern Pamir is associated with tectonic exhumation of central Pamir domes. Deeper exhumation coincides with western and asymmetric drainages and with higher precipitation today, suggesting an orographic effect on exhumation. A younging-southward trend of cooling ages may reflect tectonic processes. Overall, cooling ages derived from the Pamir are younger than ages recorded in Tibet, indicating younger and higher magnitudes of erosion in the Pamir.
Modification of the landscape by glacial erosion reflects the dynamic interplay of climate through temperature, precipitation, and prevailing wind direction, and tectonics through rock uplift and exhumation rate, lithology, and range and fault geometry. We investigate these relationships in the northeast Pamir Mountains using mapping and dating of moraines and terraces to determine the glacial history. We analyze modem glacial morphology to determine glacier area, spacing, headwall relief, debris cover, and equilibrium line altitude (ELA) using the area x altitude balance ratio (AABR), toe-to-headwall altitude ratio (THAR) and toe-to-summit altitude method (TSAM) for 156 glaciers and compare this to lithologic, tectonic, and climatic data We observe a pronounced asymmetry in glacial ELA, area, debris cover, and headwall relief that we interpret to reflect both structural and climatic control: glaciers on the downwind (eastern) side of the range are larger, more debris covered, have steeper headwalls, and tend to erode headward, truncating the smaller glaciers of the upwind, fault-controlled side of the range. We explain this by the transfer of moisture deep into the range as wind-blown or avalanched snow and by limitations imposed on glacial area on the upwind side of the range by the geometry of the Kongur extensional system (KES). The correspondence between rapid exhumation along the KES and maxima in glacier debris cover and headwall relief and minimums in all measures of ELA suggest that taller glacier headwalls develop in a response to more rapid exhumation rates. However, we find that glaciers in the Muji valley did not extend beyond the range front until at least 43 ka, in contrast to extensive glaciation since 300 ka in the south around the high peaks, a pattern which does not clearly reflect uplift rate. Instead, the difference in glacial history and the presence of large peaks (Muztagh Ata and Kongur Shan) with flanking glaciers likely reflects lithologic control (i.e., the location of crustal gneiss domes) and the formation of peaks that rise above the ELA and escape the glacial buzzsaw. (C) 2014 Elsevier B.V. All rights reserved.
Cenozoic magnetostratigraphy and magnetic properties of the southern Issyk-Kul basin, Kyrgyzstan
(2014)
We present paleomagnetic data from the northern flank of the Tianshan range, southeast of Lake Issyk-Kul (Kyrgyzstan). 613 cores were collected in two parallel sections with a total thickness of 960 m (Chon Kyzylsuu, CK) and 990 m Jeti Oguz, JO), as well as 48 cores at six sites in a nearby anticline. Rock magnetic analyses identify both magnetite and hematite in the fluvial-lacustrine sediments. The concentration of both minerals, the magnetite:hematite ratio, and the average magnetite grain size increase upward in both sections. Anisotropy of anhysteretic remanent magnetization defines a tectonic fabric with sub-horizontal maximum axes that parallel the strike direction together with intermediate and minimum axes that streak out about a great circle orthogonal to the maximum axes suggestive of a tectonic fabric emplaced during folding. Stepwise thermal demagnetization isolates interpretable magnetization components in 284 samples that define 26 polarity chrons in CK and 19 in JO. A positive fold test, dual polarities and systematic changes in rock-magnetic parameters with depth suggest that the high temperature magnetization component was acquired coevally with deposition. An age model based on a visual magnetostratigraphic correlation of both sections with the geomagnetic polarity time scale defines absolute ages from 26.0 to 13.3 Ma, with a fairly constant sedimentation rate of 9-10 cm/ka. A correlation based on a numerical algorithm arrives at a slightly different conclusion, with deposition ages from 25.2 to 11.0 Ma and sedimentation rates from 5 to 8 cm/ka. In comparison with sedimentation rates found at other magnetostratigraphic sections in the Tianshan realm, we infer that the sedimentary record in this part of the Issyk-Kul Basin precedes the more rapid phase of uplift of the Kyrgyz Tianshan. The onset of deposition and concomitant erosion of the adjacent Terskey Range is in good agreement with independent assessments of the exhumation history of this mountain range, with erosion increasing at 25-20 Ma and accelerating after 11-13 Ma. (C) 2014 Elsevier B.V. All rights reserved.
Differential exhumation in the Puna Plateau and Eastern Cordillera of NW Argentina is controlled by inherited paleostructures and resulting paleotopography related to the Cretaceous Salta Rift paleomargins. The Ceno zoic deformation front related to the development of the Andean retro-arc orogenic system is generally associated with >4 km of exhumation, which is recorded by Cenozoic apatite fi ssion-track (AFT) and (U-Th-[Sm])/He ages (He ages) in the Eastern Cordillera of NW Argentina. New AFT ages from the top of the Nevado de Cachi document Oligocene (ca. 28 Ma) cooling, which, combined with existing data, indicates exhumation of this range between ca. 28 Ma and ca. 14 Ma. However, some of the highest ranges in the Eastern Cordillera preserve Cretaceous ages indicative of limited Cenozoic exhumation. Samples collected from an similar to 3-km-elevation transect along the northern part of the Sierra de Quilmes paleorift fl ank (Laguna Brava) show AFT ages between ca. 80 and ca. 50 Ma and He ages between ca. 45 and ca. 10 Ma. Another set of samples from an similar to 1-km-elevation transect farther to the southwest (La Quebrada) shows Cretaceous AFT ages between ca. 116 Ma and ca. 76 Ma, and mainly Cretaceous He ages, in agreement with AFT data. Analysis of existing AFT and He ages from the area once occupied by the Salta Rift reveals a pattern characterized by Cretaceous ages along paleorift highs and Cenozoic ages within paleorift hanging-wall basins and later foreland basin depocenters. This pattern is interrupted by the Sierras Pampeanas at similar to 28 degrees S, which record mid-Cenozoic ages. Our data are consistent with a complex inherited pattern of pre-Andean paleostructures, likely associated with paleotopography, which was beveled by the Cenozoic regional foreland basin and reactivated during the late Neogene (ca. <10 Ma), strongly controlling the magnitude of Cenozoic uplift and exhumation and thus cooling age distribution. This, combined with variable lithologic erodibility, resulted in an irregular distribution of themochronological ages.
New low-temperature thermochronological data from 80 samples in eastern Kyrgyzstan are combined with previously published data from 61 samples to constrain exhumation in a number of mountain ranges in the Central Kyrgyz Tien Shan. All sampled ranges are found to have a broadly consistent Cenozoic exhumation history, characterized by initially low cooling rates (<1 degrees C/Myr) followed by a series of increases in exhumation that occurred diachronously across the region in the late Cenozoic that are interpreted to record the onset of deformation in different mountain ranges. Combined with geological estimates for the onset of proximal deformation, our data suggest that the Central Kyrgyz Tien Shan started deforming in the late Oligocene-early Miocene, leading to the development of several, widely spaced mountain ranges separated by large intermontane basins. Subsequently, more ranges have been constructed in response to significant shortening increases across the Central Kyrgyz Tien Shan, notably in the late Miocene. The order of range construction is interpreted to reflect variations in the susceptibility of inherited structures to reactivation. Reactivated structures are also shown to have significance along strike variations in fault vergence and displacement, which have influenced the development and growth of individual mountain ranges. Moreover, the timing of deformation allows the former extent of many intermontane basins that have since been partitioned to be inferred; this can be linked to the highly time-transgressive onset of late Cenozoic coarse clastic sedimentation.
The literature on the costing of mitigation measures for reducing impacts of natural hazards is rather fragmented. This paper provides a concise overview of the current state of knowledge in Europe on the costing of mitigation measures for the reduction of natural hazard risks (droughts, floods, storms and induced coastal hazards as well as alpine hazards) and identifies knowledge gaps and related research recommendations. Furthermore, it provides a taxonomy of related mitigation options, classifying them into nine categories: (1) management plans, land-use planning, and climate adaptation; (2) hazard modification; (3) infrastructure; (4) mitigation measures (stricto sensu); (5) communication in advance of events; (6) monitoring and early warning systems; (7) emergency response and evacuation; (8) financial incentives; and (9) risk transfer (including insurance). It is found that the costing of mitigation measures in European and in other countries has almost exclusively focused on estimating direct costs. A cost assessment framework that addresses a range of costs, possibly informed by multiple stakeholders, would provide more accurate estimates and could provide better guidance to decision makers. (C) 2014 American Society of Civil Engineers.
Costing natural hazards
(2014)
Comparing thermal wave function methods for multi-configuration time-dependent Hartree simulations
(2014)
We compare two methods for creating stochastic temperature wave functions that can be used for Multi-Configuration Time-Dependent Hartree (MCTDH) simulations. In the first method, the MCTDH coefficients are chosen randomly, while the other method uses a single Hartree product of random single-particle functions (SPFs). We find that using random SPFs dramatically improves convergence for a model system for surface sticking.
The study presents the first evidence of metamorphosed submarine ashes in the Edough Massif, in northeastern Algeria. It occurs below the greenschist-facies Tellian units that represent the thrusted Mesozoic to Eocene passive paleomargin of northern Africa deposited on thinned continental crust. The metamorphic complex consists of tectonically superposed units composed of gneisses (lower unit) and micaschists (upper unit). At the Cap de Garde, these units enclose an "intermediate unit" composed of micaschists and meter-thick layers of marbles, which are sometimes intercalated with amphibolites. The latter occur as discontinuous small lenses and layers. The amphibolites are parallel to the primary bedding of the marbles and the main foliation. Chemical markers and field observations indicate that they are metamorphic equivalents of basic igneous rocks. The lenticular character, low thickness and multiple intercalations with marine sediments and the unusual high lithium concentrations suggest subaqueous near-source basaltic ash-fall deposits in a marine environment. (C) 2014 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.
Magmatism forming the Central Anatolian Volcanic Province of Cappadocia, central Turkey, records the last phase of Neotethyan subduction after similar to 11 Ma. Thirteen large calc-alkaline ignimbrite sheets form marker bands within the volcano-sedimentary succession (the Urgup Formation) and provide a robust chronostratigraphy for paleoecologic evaluation of the interleaved paleosols. This paper evaluates the chronologic record in the context of the radiometric, magnetostratigraphic and lithostratigraphic controls. Previous inconsistencies relating primarily to K/Ar evidence were reason for the initiation of an integrated study which includes Ar-40/Ar-39 dating, palaeomagnetic and stratigraphic evidence. The newly determined Ar-40/Ar-39-ages (Lepetit, 2010) are in agreement with Ar/Ar and U/Pb data meanwhile published by Pauquette and Le Pennec (2012) and Aydar et al. (2012). The Ar-40/Ar-39-ages restrict the end of the Urgup Formation to the late Miocene. The paleosol sequence enclosed by the ignimbrites is thus restricted to the late Miocene, the most intense formation of pedogene calcretes correlating with the Messinian Salinity Crisis.
Sedimentary basin fills along the windward flanks of orogenic plateaus are valuable archives of paleoenvironmental change with the potential to resolve the history of surface uplift and orographic barrier formation. The intermontane basins of the southern Central Andes contain thick successions of sedimentary material that are commonly interbedded with datable volcanic ashes. We relate variations in the hydrogen isotopic composition of hydrated volcanic glass (delta D-g) of Neogene to Quaternary fills in the semiarid intermontane Humahuaca Basin (Eastern Cordillera, northwest Argentina) to spatiotemporal changes in topography and associated orographic effects. delta D values from volcanic glass in the basin strata (-117 parts per thousand to -98 parts per thousand) show two main trends that accompany observed tectonosedimentary events in the study area. Between 6.0 and 3.5 Ma, delta D-g values decrease by similar to 17 parts per thousand; this is associated with surface uplift in the catchment area. After 3.5 Ma, delta D-g values show abrupt deuterium enrichment, which we associate with (1) the attainment of threshold elevations for blocking moisture transport in the basin-bounding ranges to the east, and (2) the onset of semiarid conditions in the basin. Such orographic barriers throughout the eastern flanks of the Central Andes have impeded moisture transport into the orogen interior; this has likely helped maintain aridity and internal drainage conditions on the adjacent Andean Plateau.
A valley-filling ignimbrite re-exposed through subsequent river incision at the southern margin of the Andean (Puna) plateau preserves pristine geological evidence of pre-late Miocene palaeotopography in the north western Argentine Andes. Our new Ar-40/(39) Ar dating of the Las Papas Ignimbrites yields a plateau age of 9.24 +/- 0.03 Ma, indicating valley-relief and orographic-barrier conditions comparable to the present-day. A later infill of Plio-Pleistocene coarse conglomerates has been linked to wetter conditions, but resulted in no additional net incision of the Las Papas valley, considering that the base of the ignimbrite remains unexposed in the valley bottom. Our observations indicate that at least 550 m of local plateau margin relief (and likely > 2 km) existed by 9 Ma at the southern Puna margin, which likely aided the efficiency of the orographic barrier to rainfall along the eastern and south eastern flanks of the Puna and causes aridity in the plateau interior.
This study explores the potential for regions to shift to a local food supply using food self-sufficiency (FSS) as an indicator. We considered a region food self-sufficient when its total calorie production is enough to meet its demand. For future scenarios, we considered population growth, dietary changes, improved feed conversion efficiency, climate change, and crop yield increments. Starting at the 5' resolution, we investigated FSS from the lowest administrative levels to continents. Globally, about 1.9 billion people are self-sufficient within their 5' grid, while about 1 billion people from Asia and Africa require cross-continental agricultural trade in 2000. By closing yield gaps, these regions can achieve FSS, which also reduces international trade and increases a self-sufficient population in a 5' grid to 2.9 billion. The number of people depending on international trade will vary between 1.5 and 6 billion by 2050. Climate change may increase the need for international agricultural trade by 4% to 16%.
Climate impacts on human livelihoods: where uncertainty matters in projections of water availability
(2014)
Climate change will have adverse impacts on many different sectors of society, with manifold consequences for human livelihoods and well-being. However, a systematic method to quantify human well-being and livelihoods across sectors is so far unavailable, making it difficult to determine the extent of such impacts. Climate impact analyses are often limited to individual sectors (e.g. food or water) and employ sector-specific target measures, while systematic linkages to general livelihood conditions remain unexplored. Further, recent multi-model assessments have shown that uncertainties in projections of climate impacts deriving from climate and impact models, as well as greenhouse gas scenarios, are substantial, posing an additional challenge in linking climate impacts with livelihood conditions. This article first presents a methodology to consistently measure what is referred to here as AHEAD (Adequate Human livelihood conditions for wEll-being And Development). Based on a trans-disciplinary sample of concepts addressing human well-being and livelihoods, the approach measures the adequacy of conditions of 16 elements. We implement the method at global scale, using results from the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) to show how changes in water availability affect the fulfilment of AHEAD at national resolution. In addition, AHEAD allows for the uncertainty of climate and impact model projections to be identified and differentiated. We show how the approach can help to put the substantial inter-model spread into the context of country-specific livelihood conditions by differentiating where the uncertainty about water scarcity is relevant with regard to livelihood conditions - and where it is not. The results indicate that livelihood conditions are compromised by water scarcity in 34 countries. However, more often, AHEAD fulfilment is limited through other elements. The analysis shows that the water-specific uncertainty ranges of the model output are outside relevant thresholds for AHEAD for 65 out of 111 countries, and therefore do not contribute to the overall uncertainty about climate change impacts on livelihoods. In 46 of the countries in the analysis, water-specific uncertainty is relevant to AHEAD. The AHEAD method presented here, together with first results, forms an important step towards making scientific results more applicable for policy decisions.
Water is an essential input to the majority of human activities. Often, access to sufficient water resources is limited by quality and infrastructure aspects, rather than by resource availability alone, and each activity has different requirements regarding the nature of these aspects. This paper develops an integrated approach to assess the adequacy of water resources for the three major water users: the domestic, agricultural and industrial sectors. Additionally, we include environmental water requirements. We first outline the main determinants of water adequacy for each sector. Subsequently, we present an integrated approach using fuzzy logic, which allows assessing sector-specific as well as overall water adequacy. We implement the approach in two case study settings to exemplify the main features of the approach. Using results from two climate models and two forcing RCPs (representative concentration pathways), as well as population projections, we further assess the impacts of climate change in combination with population growth on the adequacy of water resources. The results provide an important step forward in determining the most relevant factors, impeding adequate access to water, which remains an important challenge in many regions of the world. The methodology allows one to directly identify the factors that are most decisive in determining the adequacy of water in each region, pointing towards the most efficient intervention points to improve conditions. Our findings underline the fact that, in addition to water volumes, water quality is a limitation for all sectors and, especially for the environmental sector, high levels of pollution are a threat to water adequacy.
Changes in the Atlantic overturning circulation have a strong influence on European temperatures, North American sea level and other climate phenomena worldwide. A meaningful assessment of associated societal impacts needs to be based on the full range of its possible future evolution. This requires capturing both the uncertainty in future warming pathways and the inherently long-term response of the ocean circulation. While probabilistic projections of the global mean and regional temperatures exist, process-based probabilistic assessments of large-scale dynamical systems such as the Atlantic overturning are still missing. Here we present such an assessment and find that a reduction of more than 50 % in Atlantic overturning strength by the end of the 21 (s t) century is within the likely range under an unmitigated climate change scenario (RCP8.5). By combining linear response functions derived from comprehensive climate simulations with the full range of possible future warming pathways, we provide probability estimates of overturning changes by the year 2100. A weakening of more than 25 % is found to be very unlikely under a climate protection scenario (RCP2.6), but likely for unmitigated climate change. The method is able to reproduce the modelled recovery caused by climatic equilibration under climate protection scenarios which provides confidence in the approach. Within this century, a reduction of the Atlantic overturning is a robust climatic phenomena that intensifies with global warming and needs to be accounted for in global adaptation strategies.
Coastal flood damage and adaptation costs under 21st century sea-level rise are assessed on a global scale taking into account a wide range of uncertainties in continental topography data, population data, protection strategies, socioeconomic development and sea-level rise. Uncertainty in global mean and regional sea level was derived from four different climate models from the Coupled Model Intercomparison Project Phase 5, each combined with three land-ice scenarios based on the published range of contributions from ice sheets and glaciers. Without adaptation, 0.2-4.6% of global population is expected to be flooded annually in 2100 under 25-123 cm of global mean sea-level rise, with expected annual losses of 0.3-9.3% of global gross domestic product. Damages of this magnitude are very unlikely to be tolerated by society and adaptation will be widespread. The global costs of protecting the coast with dikes are significant with annual investment and maintenance costs of US$ 12-71 billion in 2100, but much smaller than the global cost of avoided damages even without accounting for indirect costs of damage to regional production supply. Flood damages by the end of this century are much more sensitive to the applied protection strategy than to variations in climate and socioeconomic scenarios as well as in physical data sources (topography and climate model). Our results emphasize the central role of long-term coastal adaptation strategies. These should also take into account that protecting large parts of the developed coast increases the risk of catastrophic consequences in the case of defense failure.
Making confident statements about the evolution of an ice-sheet shelf system with a numerical model requires the capability to reproduce the migration of the grounding line. Here we show that the shallow-ice approximation/shallow-shelf approximation hybrid-type Parallel Ice Sheet Model (PISM), with its recent improvements, is capable of modeling the grounding line motion in a perturbed ice-sheet shelf system. The model is set up according to the three-dimensional Marine Ice-Sheet Model Intercomparison Project (MISMIP3d), and simulations are carried out across a broad range of spatial resolutions. Using (1) a linear interpolation of the grounding line with locally interpolated basal friction and (2) an improved driving-stress computation across the grounding line, the reversibility of the grounding line (i.e. its retreat after an advance forced by a local perturbation of basal resistance) is captured by the model even at medium and low resolutions (Delta x > 10 km). The transient model response is qualitatively similar to that of higher-order models but reveals a higher initial sensitivity to perturbations on very short timescales. Our findings support the application of PISM to the Antarctic ice sheet from regional up to continental scales and on relatively low spatial resolutions.
The Isabena catchment (445 km(2)), Spain, features highly diverse spatial heterogeneity in land use, lithology and rainfall. Consequently, the relative contribution in terms of water and sediment yield varies immensely between its subcatchments, and also temporally. This study presents the synthesis of similar to 2.5 years of monitoring rainfall, discharge and suspended sediment concentration (SSC) in the five main subcatchments of the Isabena and its outlet.
Continuous discharge at the subcatchment outlets, nine tipping bucket rainfall and automatic SSC samplers (complemented by manual samples), were collected from June 2011 until November 2013. The water stage records were converted to discharge using a rating curve derived with Bayesian regression. For reconstructing sediment yields, the data from the intermittent SSC sampling needed to be interpolated. We employed non-parametric multivariate regression (Quantile Regression Forests, QRF) using the discharge and rainfall data plus different aggregation levels of these as ancillary predictors. The subsequent Monte Carlo simulations allowed the determination of monthly sediment yields and their uncertainty.
The Isabena catchment shows high erosion dynamics with great variability in space and time, with stark contrasts even between adjacent subcatchments. The natural conditions make water and sediment monitoring and instrumentation very challenging; the measurement of discharge is particularly prone to considerable uncertainties. The QRF method employed for reconstructing sedigraphs and monthly yields proved well suited for the task.
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.
Some of the largest and most erosive floods on Earth result from the failure of glacial dams. While potentially cataclysmic ice dams are recognized to have repeatedly formed along ice-sheet margins, much less is known about the frequency and longevity of ice dams caused by mountain glaciers, and their impact on landscape evolution. Here we present field observations and results from cosmogenic nuclide dating that allow reconstructing a > 100-k.y.-long history of glacial damming in the Shyok Valley, eastern Karakoram (South Asia). Our field observations provide evidence that Asia's second-longest glacier, the Siachen, once extended for over 180 km and blocked the Shyok River during the penultimate glacial period, leading to upstream deposition of a more than 400-m-thick fluvio-lacustrine valley fill. Be-10-depth profile modeling indicates that glacial damming ended with the onset of the Eemian interglacial and that the Shyok River subsequently incised the valley fill at an average rate of similar to 4-7 m k.y.(-1). Comparison with contemporary ice-dammed lakes in the Karakoram and elsewhere suggests recurring outburst floods during the aggradation period, while over 25 cycles of fining-upward lake deposits within the valley fill indicate impounding of floods from farther upstream. Despite prolonged damming, the net effect of this and probably earlier damming episodes by the Siachen Glacier is dominated by glacial erosion in excess of fluvial incision, as evidenced by a pronounced overdeepening that follows the glaciated valley reach. Strikingly similar overdeepened valleys at all major confluences of the Shyok and Indus Rivers with Karakoram tributaries indicate that glacial dams and subsequent outburst floods have been widespread and frequent in this region during the Quaternary. Our study suggests that the interaction of Karakoram glaciers with the Shyok and Indus Rivers promoted valley incision and headward erosion into the western margin of the Tibetan Plateau.
Sr isotope records from marginal marine basins track the mixing between seawater and local continental runoff, potentially recording the effects of sea level, tectonic, and climate forcing in marine fossils and sediments. Our 110 new Sr-87/Sr-86 analyses on oyster and foraminifera samples from six late Miocene stratigraphic sections in southern Turkey, Crete, and Sicily show that Sr-87/Sr-86 fell below global seawater values in the basins several million years before the Messinian Salinity Crisis, coinciding with tectonic uplift and basin shallowing. 87Sr/86Sr from more centrally located basins (away from the Mediterranean coast) drop below global seawater values only during the Messinian Salinity Crisis. In addition to this general trend, 55 new Sr-87/Sr-86 analyses from the astronomically tuned Lower Evaporites in the central Apennines (Italy) allow us to explore the effect of glacio-eustatic sea level and precipitation changes on Sr-87/Sr-86. Most variation in our data can be explained by changes in sea level, with greatest negative excursions from global seawater values occurring during relative sea level lowstands, which generally coincided with arid conditions in the Mediterranean realm. We suggest that this greater sensitivity to lowered sea level compared with higher runoff could relate to the inverse relationship between Sr concentration and river discharge. Variations in the residence time of groundwater within the karst terrain of the circum-Mediterranean region during arid and wet phases may help to explain the single (robust) occurrence of a negative excursion during a sea level highstand, but this explanation remains speculative without more detailed paleoclimatic data for the region.
Erosion in the Himalaya is responsible for one of the greatest mass redistributions on Earth and has fueled models of feedback loops between climate and tectonics. Although the general trends of erosion across the Himalaya are reasonably well known, the relative importance of factors controlling erosion is less well constrained. Here we present 25 Be-10-derived catchment-averaged erosion rates from the Yamuna catchment in the Garhwal Himalaya, northern India. Tributary erosion rates range between similar to 0.1 and 0.5mmyr(-1) in the Lesser Himalaya and similar to 1 and 2mmyr(-1) in the High Himalaya, despite uniform hillslope angles. The erosion-rate data correlate with catchment-averaged values of 5 km radius relief, channel steepness indices, and specific stream power but to varying degrees of nonlinearity. Similar nonlinear relationships and coefficients of determination suggest that topographic steepness is the major control on the spatial variability of erosion and that twofold to threefold differences in annual runoff are of minor importance in this area. Instead, the spatial distribution of erosion in the study area is consistent with a tectonic model in which the rock uplift pattern is largely controlled by the shortening rate and the geometry of the Main Himalayan Thrust fault (MHT). Our data support a shallow dip of the MHT underneath the Lesser Himalaya, followed by a midcrustal ramp underneath the High Himalaya, as indicated by geophysical data. Finally, analysis of sample results from larger main stem rivers indicates significant variability of Be-10-derived erosion rates, possibly related to nonproportional sediment supply from different tributaries and incomplete mixing in main stem channels.
Injection of nanoscale zero-valent iron (nZVI) has recently gained great interest as emerging technology for in-situ remediation of chlorinated organic compounds from groundwater systems. Zero-valent iron (ZVI) is able to reduce organic compounds and to render it to less harmful substances. The use of nanoscale particles instead of granular or microscale particles can increase dechlorination rates by-orders of magnitude due to its high surface area. However, classical nZVI appears to be hampered in its environmental application by its limited mobility. One approach is colloid supported transport of nZVI, where the nZVI gets transported by a Mobile colloid. In this study transport properties of activated carbon colloid supported nZVI (c-nZVI; d(50) = 2.4 mu m) are investigated in column tests using columns of 40 cm length, which were filled with porous media. A suspension was pumped through the column under different physicochemical conditions (addition of a polyanionic stabilizer and changes in pH and ionic strength). Highest observed breakthrough was 62% of the injected concentration in glass beads with addition of stabilizer. Addition of mono- and bivalent salt, e.g. more than 0.5 mM/L CaCl2, can decrease mobility and changes in pH to values below six can inhibit mobility at all. Measurements of colloid sizes and zeta potentials show changes in the mean particle size by a factor of ten and an increase of zeta potential from -62 mV to -80 mV during the transport experiment. However, results suggest potential applicability of c-nZVI under field conditions. (C) 2014 Elsevier B.V. All rights reserved.
We present the results of biogeochemical and mineralogical analyses on a sediment core that covers the Holocene sedimentation history of the climatically sensitive, closed, saline, and alkaline Lonar Lake in the core monsoon zone in central India. We compare our results of C/N ratios, stable carbon and nitrogen isotopes, grain-size, as well as amino acid derived degradation proxies with climatically sensitive proxies of other records from South Asia and the North Atlantic region. The comparison reveals some more or less contemporaneous climate shifts. At Lonar Lake, a general long term climate transition from wet conditions during the early Holocene to drier conditions during the late Holocene, delineating the insolation curve, can be reconstructed. In addition to the previously identified periods of prolonged drought during 4.6-3.9 and 2.0-0.6 cal ka that have been attributed to temperature changes in the Indo Pacific Warm Pool, several additional phases of shorter term climate alteration superimposed upon the general climate trend can be identified. These correlate with cold phases in the North Atlantic region. The most pronounced climate deteriorations indicated by our data occurred during 62-5.2,4.6-3.9, and 2.0-0.6 cal ka BP. The strong dry phase between 4.6 and 3.9 cal ka BP at Lonar Lake corroborates the hypothesis that severe climate deterioration contributed to the decline of the Indus Civilisation about 3.9 ka BP. (C) 2014 Elsevier B.V. All rights reserved.
Background and Aims Dynamic processes occurring at the soil-root interface crucially influence soil physical, chemical and biological properties at a local scale around the roots, and are technically challenging to capture in situ. This study presents a novel multi-imaging approach combining fluorescence and neutron radiography that is able to simultaneously monitor root growth, water content distribution, root respiration and root exudation.
Methods Germinated seeds of white lupins (Lupinus albus) were planted in boron-free glass rhizotrons. After 11 d, the rhizotrons were wetted from the bottom and time series of fluorescence and neutron images were taken during the subsequent day and night cycles for 13 d. The following day (i.e. 25 d after planting) the rhizotrons were again wetted from the bottom and the measurements were repeated. Fluorescence sensor foils were attached to the inner sides of the glass and measurements of oxygen and pH were made on the basis of fluorescence intensity. The experimental set-up allowed for simultaneous fluorescence imaging and neutron radiography.
Key Results The interrelated patterns of root growth and distribution in the soil, root respiration, exudation and water uptake could all be studied non-destructively and at high temporal and spatial resolution. The older parts of the root system with greater root-length density were associated with fast decreases of water content and rapid changes in oxygen concentration. pH values around the roots located in areas with low soil water content were significantly lower than the rest of the root system.
Conclusions The results suggest that the combined imaging set-up developed here, incorporating fluorescence intensity measurements, is able to map important biogeochemical parameters in the soil around living plants with a spatial resolution that is sufficiently high enough to relate the patterns observed to the root system.
We used inverse modelling techniques and soil moisture measured by the cosmic-ray neutron sensing (CRS) to estimate root-zone soil hydraulic properties at the field scale. A HYDRUS-1D model was developed for inverse modelling and calibrated with parameter estimation software (PEST) using a global optimizer. Integral CRS measurements recorded from a sunflower farm in Germany comprised the model input. Data were transformed to soil water storage to enable direct model calibration with a HYDRUS soil-water balance. Effective properties at the CRS scale were compared against local measurements and other inversely estimated soil properties from independent soil moisture profiles. Moreover, CRS-scale soil properties were tested on the basis of how field soil moisture (vertical distribution) and soil water storage were reproduced. This framework provided good estimates of effective soil properties at the CRS scale. Simulated soil moisture at different depths at the CRS scale agreed with field observations. Moreover, simulated soil water storage at the CRS scale compared well with calculations from point-scale profiles, despite their different support volumes. The CRS-scale soil properties estimated with the inverse model were within the range of variation of properties identified from all inverse simulations at the local scale. This study demonstrates the potential of CRS for inverse estimation of soil hydraulic properties.
Nanoscale zero-valent iron (nZVI) has recently gained great interest in the scientific community as in situ reagent for installation of permeable reactive barriers in aquifer systems, since nZVI is highly reactive with chlorinated compounds and may render them to harmless substances. However, nZVI has a high tendency to agglomerate and sediment; therefore it shows very limited transport ranges. One new approach to overcome the limited transport of nZVI in porous media is using a suited carrier colloid. In this study we tested mobility of a carbon colloid supported nZVI particle "Carbo-Iron Colloids" (CIC) with a mean size of 0.63 mu m in a column experiment of 40 cm length and an experiment in a two-dimensional (2D) aquifer test system with dimensions of 110x40x5 cm. Results show a breakthrough maximum of 82 % of the input concentration in the column experiment and 58 % in the 2D-aquifer test system. Detected residuals in porous media suggest a strong particle deposition in the first centimeters and few depositions in the porous media in the further travel path. Overall, this suggests a high mobility in porous media which might be a significant enhancement compared to bare or polyanionic stabilized nZVI.
Members of the acI lineage of Actinobacteria are the most abundant microorganisms in most freshwater lakes; however, our understanding of the keys to their success and their role in carbon and nutrient cycling in freshwater systems has been hampered by the lack of pure cultures and genomes. We obtained draft genome assemblies from 11 single cells representing three acI tribes (acI-A1, acI-A7, acI-B1) from four temperate lakes in the United States and Europe. Comparative analysis of acI SAGs and other available freshwater bacterial genomes showed that acI has more gene content directed toward carbohydrate acquisition as compared to Polynucleobacter and LD12 Alphaproteobacteria, which seem to specialize more on carboxylic acids. The acI genomes contain actinorhodopsin as well as some genes involved in anaplerotic carbon fixation indicating the capacity to supplement their known heterotrophic lifestyle. Genome-level differences between the acI-A and acI-B clades suggest specialization at the clade level for carbon substrate acquisition. Overall, the acI genomes appear to be highly streamlined versions of Actinobacteria that include some genes allowing it to take advantage of sunlight and N-rich organic compounds such as polyamines, di-and oligopeptides, branched-chain amino acids and cyanophycin. This work significantly expands the known metabolic potential of the cosmopolitan freshwater acI lineage and its ecological and genetic traits.
Sedimentary proxy records constitute a significant portion of the recorded evidence that allows us to investigate paleoclimatic conditions and variability. However, uncertainties in the dating of proxy archives limit our ability to fix the timing of past events and interpret proxy record intercomparisons. While there are various age-modeling approaches to improve the estimation of the age-depth relations of archives, relatively little focus has been placed on the propagation of the age (and radiocarbon calibration) uncertainties into the final proxy record.
We present a generic Bayesian framework to estimate proxy records along with their associated uncertainty, starting with the radiometric age-depth and proxy-depth measurements, and a radiometric calibration curve if required. We provide analytical expressions for the posterior proxy probability distributions at any given calendar age, from which the expected proxy values and their uncertainty can be estimated. We illustrate our method using two synthetic data sets and then use it to construct the proxy records for groundwater inflow and surface erosion from Lonar lake in central India.
Our analysis reveals interrelations between the uncertainty of the proxy record over time and the variance of proxies along the depth of the archive. For the Lonar lake proxies, we show that, rather than the age uncertainties, it is the proxy variance combined with calibration uncertainty that accounts for most of the final uncertainty. We represent the proxy records as probability distributions on a precise, error-free timescale that makes further time series analyses and intercomparisons of proxies relatively simple and clear. Our approach provides a coherent understanding of age uncertainties within sedimentary proxy records that involve radiometric dating. It can be potentially used within existing age modeling structures to bring forth a reliable and consistent framework for proxy record estimation.
This paper employs a complex network approach to determine the topology and evolution of the network of extreme precipitation that governs the organization of extreme rainfall before, during, and after the Indian Summer Monsoon (ISM) season. We construct networks of extreme rainfall events during the ISM (June-September), post-monsoon (October-December), and pre-monsoon (March-May) periods from satellite-derived (Tropical Rainfall Measurement Mission, TRMM) and rain-gauge interpolated (Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources, APHRODITE) data sets. The structure of the networks is determined by the level of synchronization of extreme rainfall events between different grid cells throughout the Indian subcontinent. Through the analysis of various complex-network metrics, we describe typical repetitive patterns in North Pakistan (NP), the Eastern Ghats (EG), and the Tibetan Plateau (TP). These patterns appear during the pre-monsoon season, evolve during the ISM, and disappear during the post-monsoon season. These are important meteorological features that need further attention and that may be useful in ISM timing and strength prediction.
Recurrence-plot-based recurrence networks are an approach used to analyze time series using a complex networks theory. In both approaches - recurrence plots and recurrence networks -, a threshold to identify recurrent states is required. The selection of the threshold is important in order to avoid bias of the recurrence network results. In this paper, we propose a novel method to choose a recurrence threshold adaptively. We show a comparison between the constant threshold and adaptive threshold cases to study period-chaos and even period-period transitions in the dynamics of a prototypical model system. This novel method is then used to identify climate transitions from a lake sediment record.
Deforestation is a prominent anthropogenic cause of erosive overland flow and slope instability, boosting rates of soil erosion and concomitant sediment flux. Conventional methods of gauging or estimating post-logging sediment flux often focus on annual timescales but overlook potentially important process response on shorter intervals immediately following timber harvest. We resolve such dynamics with non-parametric quantile regression forests (QRF) based on high-frequency (3 min) discharge measurements and sediment concentration data sampled every 30-60 min in similar-sized (similar to 0.1 km(2)) forested Chilean catchments that were logged during either the rainy or the dry season. The method of QRF builds on the random forest algorithm, and combines quantile regression with repeated random sub-sampling of both cases and predictors. The algorithm belongs to the family of decision-tree classifiers, which allow quantifying relevant predictors in high-dimensional parameter space. We find that, where no logging occurred, similar to 80% of the total sediment load was transported during extremely variable runoff events during only 5% of the monitoring period. In particular, dry-season logging dampened the relative role of these rare, extreme sediment-transport events by increasing load efficiency during more efficient moderate events. We show that QRFs outperform traditional sediment rating curves (SRCs) in terms of accurately simulating short-term dynamics of sediment flux, and conclude that QRF may reliably support forest management recommendations by providing robust simulations of post-logging response of water and sediment fluxes at high temporal resolution.
Subsurface stormflow is thought to occur mainly in humid environments with steep terrains. However, in semi-arid areas, preferential flow through macropores can also result in a significant contribution of subsurface stormflow to catchment runoff for varying catchment conditions. Most hydrological models neglect this important subsurface preferential flow. Here, we use the process-oriented hydrological model Hillflow-3D, which includes a macropore flow approach, to simulate rainfall-runoff in the semi-arid Parapunos catchment in Spain, where macropore flow was observed in previous research. The model was extended for this study to account for sorptivity under very dry soil conditions. The results of the model simulations with and without macropore flow are compared. Both model versions give reasonable results for average rainfall situations, although the approach with the macropore concept provides slightly better results. The model results for scenarios of extreme rainfall events (>13.3mm30min(-1)) however show large differences between the versions with and without macropores. These model results compared with measured rainfall-runoff data show that the model with the macropore concept is better. Our conclusion is that preferential flow is important in controlling surface runoff in case of specific, high intensity rainfall events. Therefore, preferential flow processes must be included in hydrological models where we know that preferential flow occurs. Hydrological process models with a less detailed process description may fit observed average events reasonably well but can result in erroneous predictions for more extreme events. Copyright (c) 2013 John Wiley & Sons, Ltd.
Two lines of research are combined in this study: first, the development of tools for the temporal disaggregation of precipitation, and second, some newer results on the exponential scaling of heavy short-term precipitation with temperature, roughly following the Clausius-Clapeyron (CC) relation. Having no extra temperature dependence, the traditional disaggregation schemes are shown to lack the crucial CC-type temperature dependence. The authors introduce a proof-of-concept adjustment of an existing disaggregation tool, the multiplicative cascade model of Olsson, and show that, in principal, it is possible to include temperature dependence in the disaggregation step, resulting in a fairly realistic temperature dependence of the CC type. They conclude by outlining the main calibration steps necessary to develop a full-fledged CC disaggregation scheme and discuss possible applications.
The 933 km(2) Bengue catchment in northeastern Brazil is characterized by distinct rainy and dry seasons. Precipitation is stored in variously sized reservoirs, which is essential for the local population. In this study, we used TerraSAR-X SM(HH) data for an one-year monitoring of seasonal changes in the reservoir areas from July 2011 to July 2012. The monitoring was based on acquisitions in the ascending pass direction, complemented by occasional descending-pass images. To detect water surface areas, a histogram analysis followed by a global threshold classification was performed, and the results were validated using in situ GPS data. Distinguishing between small reservoirs and similar looking dark areas was difficult. Therefore, we tested several approaches for identifying misclassified areas. An analysis of the surface area dynamics of the reservoirs indicated high spatial and temporal heterogeneities and a large decrease in the total water surface area of the reservoirs in the catchment by approximately 30% within one year.
Knowledge of sediment sources is a prerequisite for sustainable management practices and may furthermore improve our understanding of water and sediment fluxes. Investigations have shown that a number of characteristic soil properties can be used as "fingerprints" to trace back the sources of river sediments. Spectral properties have recently been successfully used as such characteristics in fingerprinting studies. Despite being less labour-intensive than geochemical analyses, for example, spectroscopy allows measurements of small amounts of sediment material (> 60 mg), thus enabling inexpensive analyses even of intra-event variability. The focus of this study is on the examination of spectral properties of fluvial sediment samples to detect changes in source contributions, both between and within individual flood events.
Sediment samples from the following three different origins were collected in the Isabena catchment (445 km(2)) in the central Spanish Pyrenees: (1) soil samples from the main potential source areas, (2) stored fine sediment from the channel bed once each season in 2011 and (3) suspended sediment samples during four flood events in autumn 2011 and spring 2012 at the catchment outlet as well as at several subcatchment outlets. All samples were dried and measured for spectral properties in the laboratory using an ASD spectroradiometer. Colour parameters and physically based features (e.g. organic carbon, iron oxide and clay content) were calculated from the spectra. Principal component analyses (PCA) were applied to all three types of samples to determine natural clustering of samples, and a mixing model was applied to determine source contributions.
We found that fine sediment stored in the river bed seems to be mainly influenced by grain size and seasonal variability, while sampling location-and thus the effect of individual tributaries or subcatchments-seem to be of minor importance. Suspended sediment sources were found to vary between, as well as within, flood events; although badlands were always the major source. Forests and grasslands contributed little (< 10 %), and other sources (not further determinable) contributed up to 40 %. The analyses further suggested that sediment sources differ among the subcatchments and that subcatchments comprising relatively large proportions of badlands contributed most to the four flood events analyzed.
Spectral fingerprints provide a rapid and cost-efficient alternative to conventional fingerprint properties. However, a combination of spectral and conventional fingerprint properties could potentially permit discrimination of a larger number of source types.
Ecohydrology analyses the interactions of biotic and abiotic aspects of our ecosystems and landscapes. It is a highly diverse discipline in terms of its thematic and methodical research foci. This article gives an overview of current German ecohydrological research approaches within plant-animal-soil-systems, meso-scale catchments and their river networks, lake systems, coastal areas and tidal rivers. It discusses their relevant spatial and temporal process scales and different types of interactions and feedback dynamics between hydrological and biotic processes and patterns. The following topics are considered key challenges: innovative analysis of the interdisciplinary scale continuum, development of dynamically coupled model systems, integrated monitoring of coupled processes at the interface and transition from basic to applied ecohydrological science to develop sustainable water and land resource management strategies under regional and global change.
This study aims to compare impacts of climate change on streamflow in four large representative African river basins: the Niger, the Upper Blue Nile, the Oubangui and the Limpopo. We set up the eco-hydrological model SWIM (Soil and Water Integrated Model) for all four basins individually. The validation of the models for four basins shows results from adequate to very good, depending on the quality and availability of input and calibration data.
For the climate impact assessment, we drive the model with outputs of five bias corrected Earth system models of Coupled Model Intercomparison Project Phase 5 (CMIP5) for the representative concentration pathways (RCPs) 2.6 and 8.5. This climate input is put into the context of climate trends of the whole African continent and compared to a CMIP5 ensemble of 19 models in order to test their representativeness. Subsequently, we compare the trends in mean discharges, seasonality and hydrological extremes in the 21st century. The uncertainty of results for all basins is high. Still, climate change impact is clearly visible for mean discharges but also for extremes in high and low flows. The uncertainty of the projections is the lowest in the Upper Blue Nile, where an increase in streamflow is most likely. In the Niger and the Limpopo basins, the magnitude of trends in both directions is high and has a wide range of uncertainty. In the Oubangui, impacts are the least significant. Our results confirm partly the findings of previous continental impact analyses for Africa. However, contradictory to these studies we find a tendency for increased streamflows in three of the four basins (not for the Oubangui). Guided by these results, we argue for attention to the possible risks of increasing high flows in the face of the dominant water scarcity in Africa. In conclusion, the study shows that impact intercomparisons have added value to the adaptation discussion and may be used for setting up adaptation plans in the context of a holistic approach.
Identification and characterization of growing large-scale en-echelon fractures in a salt mine
(2014)
The spatiotemporal seismicity of acoustic emission (AE) events recorded in the Morsleben salt mine is investigated. Almost a year after backfilling of the cavities from 2003, microevents are distributed with distinctive stripe shapes above cavities at different depth levels. The physical forces driving the creation of these stripes are still unknown. This study aims to find the active stripes and track fracture developments over time by combining two different temporal and spatial clustering techniques into a single methodological approach. Anomalous seismicity parameters values like sharp b-value changes for two active stripes are good indicators to explain possible stress accumulation at the stripe tips. We identify the formation of two new seismicity stripes and show that the AE activities in active clusters are migrated mostly unidirectional to eastward and upward. This indicates that the growth of underlying macrofractures is controlled by the gradient of extensional stress. Studying size distribution characteristic in terms of frequency-magnitude distribution and b-value in active phase and phase with constant seismicity rate show that deviations from the Gutenberg-Richter power law can be explained by the inclusion of different activity phases: (1) the inactive period before the formation of macrofractures, which is characterized by a deficit of larger events (higher b-values) and (2) the period of fracture growth characterized by the occurrence of larger events (smaller b-values).
A spatially localized seismic sequence originated few tens of kilometres offshore the Mediterranean coast of Spain, close to the Ebro river delta, starting on 2013 September 5, and lasting at least until 2013 October. The sequence culminated in a maximal moment magnitude M-w 4.3 earthquake, on 2013 October 1. The most relevant seismogenic feature in the area is the Fosa de Amposta fault system, which includes different strands mapped at different distances to the coast, with a general NE-SW orientation, roughly parallel to the coastline. However, no significant known historical seismicity has involved this fault system in the past. The epicentral region is also located near the offshore platform of the Castor project, where gas is conducted through a pipeline from mainland and where it was recently injected in a depleted oil reservoir, at about 2 km depth. We analyse the temporal evolution of the seismic sequence and use full waveform techniques to derive absolute and relative locations, estimate depths and focal mechanisms for the largest events in the sequence (with magnitude mbLg larger than 3), and compare them to a previous event (2012 April 8, mbLg 3.3) taking place in the same region prior to the gas injection. Moment tensor inversion results show that the overall seismicity in this sequence is characterized by oblique mechanisms with a normal fault component, with a 30A degrees low-dip angle plane oriented NNE-SSW and a subvertical plane oriented NW-SE. The combined analysis of hypocentral location and focal mechanisms could indicate that the seismic sequence corresponds to rupture processes along shallow low-dip surfaces, which could have been triggered by the gas injection in the reservoir, and excludes the activation of the Amposta fault, as its known orientation is inconsistent with focal mechanism results. An alternative scenario includes the iterated triggering of a system of steep faults oriented NW-SE, which were identified by prior marine seismics investigations.
Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterized by low signal-to-noise ratio, automated location methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this paper, we propose a picking free earthquake location method based on the use of the short-term-average/long-term-average (STA/LTA) traces at different stations as observed data. For the P phases, we use the STA/LTA traces of the vertical energy function, whereas for the S phases, we use the STA/LTA traces of a second characteristic function, which is obtained using the principal component analysis technique. In order to locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a 3-D grid, we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial coordinates of the seismic event. A pilot application was performed in the Campania-Lucania region (southern Italy) using a seismic network (Irpinia Seismic Network) with an aperture of about 150 km. We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < M-L < 2.7. A subset of these locations were compared with accurate manual locations refined by using a double-difference technique. Our results indicate a good agreement with manual locations. Moreover, our method is noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals.
Permafrost-affected soils are among the most obvious ecosystems in which current microbial controls on organic matter decomposition are changing as a result of global warming. Warmer conditions in polygonal tundra will lead to a deepening of the seasonal active layer, provoking changes in microbial processes and possibly resulting in exacerbated carbon degradation under increasing anoxic conditions. To identify current microbial assemblages in carbon rich, water saturated permafrost environments, four polygonal tundra sites were investigated on Herschel Island and the Yukon Coast, Western Canadian Arctic. Ion Torrent sequencing of bacterial and archaeal 16S rRNA amplicons revealed the presence of all major microbial soil groups and indicated a local, vertical heterogeneity of the polygonal tundra soil community with increasing depth. Microbial diversity was found to be highest in the surface layers, decreasing towards the permafrost table. Quantitative PCR analysis of functional genes involved in carbon and nitrogen-cycling revealed a high functional potential in the surface layers, decreasing with increasing active layer depth. We observed that soil properties driving microbial diversity and functional potential varied in each study site. These results highlight the small-scale heterogeneity of geomorphologically comparable sites, greatly restricting generalizations about the fate of permafrost-affected environments in a warming Arctic.
Flood hazard projections under climate change are typically derived by applying model chains consisting of the following elements: "emission scenario - global climate model - downscaling, possibly including bias correction hydrological model - flood frequency analysis". To date, this approach yields very uncertain results, due to the difficulties of global and regional climate models to represent precipitation. The implementation of such model chains requires major efforts, and their complexity is high.
We propose for the Mekong River an alternative approach which is based on a shortened model chain: "emission scenario - global climate model - non-stationary flood frequency model". The underlying idea is to use a link between the Western Pacific monsoon and local flood characteristics: the variance of the monsoon drives a non-stationary flood frequency model, yielding a direct estimate of flood probabilities. This approach bypasses the uncertain precipitation, since the monsoon variance is derived from large-scale wind fields which are better represented by climate models. The simplicity of the monsoon-flood link allows deriving large ensembles of flood projections under climate change. We conclude that this is a worthwhile, complementary approach to the typical model chains in catchments where a substantial link between climate and floods is found.
We investigate the usefulness of complex flood damage models for predicting relative damage to residential buildings in a spatial and temporal transfer context. We apply eight different flood damage models to predict relative building damage for five historic flood events in two different regions of Germany. Model complexity is measured in terms of the number of explanatory variables which varies from 1 variable up to 10 variables which are singled out from 28 candidate variables. Model validation is based on empirical damage data, whereas observation uncertainty is taken into consideration. The comparison of model predictive performance shows that additional explanatory variables besides the water depth improve the predictive capability in a spatial and temporal transfer context, i.e., when the models are transferred to different regions and different flood events. Concerning the trade-off between predictive capability and reliability the model structure seem more important than the number of explanatory variables. Among the models considered, the reliability of Bayesian network-based predictions in space-time transfer is larger than for the remaining models, and the uncertainties associated with damage predictions are reflected more completely.
Flood estimation and flood management have traditionally been the domain of hydrologists, water resources engineers and statisticians, and disciplinary approaches abound. Dominant views have been shaped; one example is the catchment perspective: floods are formed and influenced by the interaction of local, catchment-specific characteristics, such as meteorology, topography and geology. These traditional views have been beneficial, but they have a narrow framing. In this paper we contrast traditional views with broader perspectives that are emerging from an improved understanding of the climatic context of floods. We come to the following conclusions: (1) extending the traditional system boundaries (local catchment, recent decades, hydrological/hydraulic processes) opens up exciting possibilities for better understanding and improved tools for flood risk assessment and management. (2) Statistical approaches in flood estimation need to be complemented by the search for the causal mechanisms and dominant processes in the atmosphere, catchment and river system that leave their fingerprints on flood characteristics. (3) Natural climate variability leads to time-varying flood characteristics, and this variation may be partially quantifiable and predictable, with the perspective of dynamic, climate-informed flood risk management. (4) Efforts are needed to fully account for factors that contribute to changes in all three risk components (hazard, exposure, vulnerability) and to better understand the interactions between society and floods. (5) Given the global scale and societal importance, we call for the organization of an international multidisciplinary collaboration and data-sharing initiative to further understand the links between climate and flooding and to advance flood research.
In aquatic systems, natural organic matter (NOM) and in particular humic substances effectively absorb the ultraviolet (UV)/visible light spectrum of solar radiation and act as a photoprotective filter for organisms. Simultaneously, UV contributes to the generation of potentially harmful reactive oxygen species (ROS). Dose-response experiments were conducted on cyanobacteria and green algae with hydrogen peroxide (H2O2) as a long-lived representative of ROS. Delayed fluorescence (DF) decay kinetics was used as a non-invasive tool to follow changes of phytoplankton activity in real time. In order to investigate phototoxicity and photoprotection by NOM on phytoplankton, we exposed algae to UV-pre-irradiated NOM and direct UV excitation. Cyanobacteria responded to H2O2 concentrations as low as 10(-7) M, while green algae were 2 orders of magnitude less sensitive. UV irradiation of medium with NOM generated H2O2 concentrations of 1.5 x 10(-7) to 3.6 x 10(-7) M. When exposed to these concentrations, only the DF of cyanobacteria led to a measurable effect while that of green algae did not change. The addition of NOM protected all phytoplankton from direct UV irradiation, but cyanobacteria benefitted less. From this we conclude that UV-irradiated water enriched with NOM can adversely affect the physiology of cyanobacteria, but not of green algae, which might control phytoplankton composition and species-specific activities.
In aquatic environments, community dynamics of bacteria, especially actively growing bacteria (AGB), are tightly linked with dissolved organic matter (DOM) quantity and quality. We analyzed the community dynamics of DNA-synthesizing and accordingly AGB by linking an improved bromodeoxyuridine immunocytochemistry approach with fluorescence-activated cell sorting (BrdU-FACS). FACS-sorted cells of even oligotrophic ecosystems in winter were characterized by 16S rRNA gene analysis. In incubation experiments, we examined community shifts of AGB in response to the addition of N-acetyl-glucosamine (NAG), one of the most abundant aminosugars in aquatic systems. Our improved BrdU-FACS analysis revealed that AGB winter communities of oligotrophic Lake Stechlin (northeastern Germany) substantially differ from those of total bacteria and consist of Alpha-, Beta-, Gamma-, Deltaproteobacteria, Actinobacteria, Candidatus OP10 and Chloroflexi. AGB populations with different BrdU-fluorescence intensities and cell sizes represented different phylotypes suggesting that single-cell growth potential varies at the taxon level. NAG incubation experiments demonstrated that a variety of widespread taxa related to Alpha-, Beta-, Gammaproteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Planctomycetes, Spirochaetes, Verrucomicrobia and Chloroflexi actively grow in the presence of NAG. The BrdU-FACS approach enables detailed phylogenetic studies of AGB and, thus, to identify those phylotypes which are potential key players in aquatic DOM cycling.
In this paper we quantify the sediment dynamics in the formerly glaciated Zielbach catchment in the Italian Alps from the end of the Last Glacial Maximum (LGM) until today. As a basis for our quantification, we use the stratigraphic record offered by a 3.5 km(2) large fan that we explore with a seismic survey, stratigraphic analyses of drillhole material, and C-14 ages measured on organic matter encountered in these drillings. In addition, we calculate past denudation rate variability in the fan deposits using concentrations of cosmogenic Be-10. We merge this information into a scenario of how the sediment flux has changed through time and how this variability can be related to climatic variations, framed within well-known paraglacial models. The results document a highly complex natural system. From the LGM to the very early Holocene, ice-melted discharge and climate variability promoted a high sediment flux (sedimentation rate up to 40 mm/yr). This flux then dramatically decreased toward interglacial values (0.8 mm/yr at 5-4 calibrated kyr B.P.). However, in contrast to the trend of classic paraglacial models, the flux recorded at Zielbach shows secondary peaks at 6.5 ka and 2.5 ka, with values of 13 mm/yr and 1.5 mm/yr, respectively. Paleo-denudation rates also decrease from similar to 33 mm/yr at the beginning of the Holocene to 0.42 mm/yr at 5 ka, with peaks of similar to 6 mm/yr and 1.1 mm/yr at 6.5 ka and 2.5 ka. High-amplitude climate change is the most likely cause of the secondary peaks, but anthropogenic activities may have contributed as well. The good correlation between paleo-sedimentation and paleo-denudation rates suggests that the majority of the deglaciated material destocked from the Zielbach catchment is stored in the alluvial fan.
This article presents comparisons among the five ground-motion models described in other articles within this special issue, in terms of data selection criteria, characteristics of the models and predicted peak ground and response spectral accelerations. Comparisons are also made with predictions from the Next Generation Attenuation (NGA) models to which the models presented here have similarities (e.g. a common master database has been used) but also differences (e.g. some models in this issue are nonparametric). As a result of the differing data selection criteria and derivation techniques the predicted median ground motions show considerable differences (up to a factor of two for certain scenarios), particularly for magnitudes and distances close to or beyond the range of the available observations. The predicted influence of style-of-faulting shows much variation among models whereas site amplification factors are more similar, with peak amplification at around 1s. These differences are greater than those among predictions from the NGA models. The models for aleatory variability (sigma), however, are similar and suggest that ground-motion variability from this region is slightly higher than that predicted by the NGA models, based primarily on data from California and Taiwan.
Marine seismology usually relies on temporary deployments of stand alone seismic ocean bottom stations (OBS), which are initialized and synchronized on ship before deployment and re-synchronized and stopped on ship after recovery several months later. In between, the recorder clocks may drift and float at unknown rates. If the clock drifts are large or not linear and cannot be corrected for, seismological applications will be limited to methods not requiring precise common timing. Therefore, for example, array seismological methods, which need very accurate timing between individual stations, would not be applicable for such deployments.
We use an OBS test-array of 12 stations and 75 km aperture, deployed for 10 months in the deep sea (4.5-5.5 km) of the mid-eastern Atlantic. The experiment was designed to analyse the potential of broad-band array seismology at the seafloor. After recovery, we identified some stations which either show unusual large clock drifts and/or static time offsets by having a large difference between the internal clock and the GPS-signal (skew).
We test the approach of ambient noise cross-correlation to synchronize clocks of a deep water OBS array with km-scale interstation distances. We show that small drift rates and static time offsets can be resolved on vertical components with a standard technique. Larger clock drifts (several seconds per day) can only be accurately recovered if time windows of one input trace are shifted according to the expected drift between a station pair before the cross-correlation. We validate that the drifts extracted from the seismometer data are linear to first order. The same is valid for most of the hydrophones. Moreover, we were able to determine the clock drift at a station where no skew could be measured. Furthermore, we find that instable apparent drift rates at some hydrophones, which are uncorrelated to the seismometer drift recorded at the same digitizer, indicate a malfunction of the hydrophone.
Crustal earthquake swarms are an expression of intensive cracking and rock damaging over periods of days, weeks or month in a small source region in the crust. They are caused by longer lasting stress changes in the source region. Often, the localized stressing of the crust is associated with fluid or gas migration, possibly in combination with pre-existing zones of weaknesses. However, verifying and quantifying localized fluid movement at depth remains difficult since the area affected is small and geophysical prospecting methods often cannot reach the required resolution.
We apply a simple and robust method to estimate the velocity ratio between compressional (P) and shear (S) waves (upsilon(P)/upsilon(S)-ratio) in the source region of an earthquake swarm. The upsilon(P)/upsilon(S)-ratio may be unusual small if the swarm is related to gas in a porous or fractured rock. The method uses arrival time difference between P and S waves observed at surface seismic stations, and the associated double differences between pairs of earthquakes. An advantage is that earthquake locations are not required and the method seems lesser dependent on unknown velocity variations in the crust outside the source region. It is, thus, suited for monitoring purposes.
Applications comprise three natural, mid-crustal (8-10 km) earthquake swarms between 1997 and 2008 from the NW-Bohemia swarm region. We resolve a strong temporal decrease of upsilon(P)/upsilon(S) before and during the main activity of the swarm, and a recovery of upsilon(P)/upsilon(S) to background levels at the end of the swarms. The anomalies are interpreted in terms of the Biot-Gassman equations, assuming the presence of oversaturated fluids degassing during the beginning phase of the swarm activity.
Hydrosedimentological studies conducted in the semiarid Upper Jaguaribe Basin, Brazil, enabled the identification of the key processes controlling sediment connectivity at different spatial scales (10(0)-10(4) km(2)).
Water and sediment fluxes were assessed from discharge, sediment concentrations and reservoir siltation measurements. Additionally, mathematical modelling (WASA-SED model) was used to quantify water and sediment transfer within the watershed.
Rainfall erosivity in the study area was moderate (4600 MJ mm ha(-1) h(-1) year(-1)), whereas runoff depths (16-60 mm year(-1)), and therefore the sediment transport capacity, were low. Consequently, similar to 60 % of the eroded sediment was deposited along the landscape, regardless of the spatial scale. The existing high-density reservoir network (contributing area of 6 km(2) per reservoir) also limits sediment propagation, retaining up to 47 % of the sediment at the large basin scale. The sediment delivery ratio (SDR) decreased with the spatial scale; on average, 41 % of the eroded sediment was yielded from the hillslopes, while for the whole 24,600-km(2) basin, the SDR was reduced to 1 % downstream of a large reservoir (1940-hm(3) capacity).
Hydrological behaviour in the Upper Jaguaribe Basin represents a constraint on sediment propagation; low runoff depth is the main feature breaking sediment connectivity, which limits sediment transference from the hillslopes to the drainage system. Surface reservoirs are also important barriers, but their relative importance to sediment retention increases with scale, since larger contributing areas are more suitable for the construction of dams due to higher hydrological potential.
In a study from 2008, Lariviere and colleagues showed, for the field of natural sciences and engineering, that the median age of cited references is increasing over time. This result was considered counterintuitive: with the advent of electronic search engines, online journal issues and open access publications, one could have expected that cited literature is becoming younger. That study has motivated us to take a closer look at the changes in the age distribution of references that have been cited in water resources journals since 1965. Not only could we confirm the findings of Lariviere and colleagues. We were also able to show that the aging is mainly happening in the oldest 10-25% of an average reference list. This is consistent with our analysis of top-cited papers in the field of water resources. Rankings based on total citations since 1965 consistently show the dominance of old literature, including text books and research papers in equal shares. For most top-cited old-timers, citations are still growing exponentially. There is strong evidence that most citations are attracted by publications that introduced methods which meanwhile belong to the standard toolset of researchers and practitioners in the field of water resources. Although we think that this trend should not be overinterpreted as a sign of stagnancy, there might be cause for concern regarding how authors select their references. We question the increasing citation of textbook knowledge as it holds the risk that reference lists become overcrowded, and that the readability of papers deteriorates.
We suggest a new clustering approach to classify focal mechanisms from large moment tensor catalogues, with the purpose of automatically identify families of earthquakes with similar source geometry, recognize the orientation of most active faults, and detect temporal variations of the rupture processes. The approach differs in comparison to waveform similarity methods since clusters are detected even if they occur in large spatial distances. This approach is particularly helpful to analyse large moment tensor catalogues, as in microseismicity applications, where a manual analysis and classification is not feasible. A flexible algorithm is here proposed: it can handle different metrics, norms, and focal mechanism representations. In particular, the method can handle full moment tensor or constrained source model catalogues, for which different metrics are suggested. The method can account for variable uncertainties of different moment tensor components. We verify the method with synthetic catalogues. An application to real data from mining induced seismicity illustrates possible applications of the method and demonstrate the cluster detection and event classification performance with different moment tensor catalogues. Results proof that main earthquake source types occur on spatially separated faults, and that temporal changes in the number and characterization of focal mechanism clusters are detected. We suggest that moment tensor clustering can help assessing time dependent hazard in mines.
Knowledge of the origin of suspended sediment is important for improving our understanding of sediment dynamics and thereupon support of sustainable watershed management. An direct approach to trace the origin of sediments is the fingerprinting technique. It is based on the assumption that potential sediment sources can be discriminated and that the contribution of these sources to the sediment can be determined on the basis of distinctive characteristics (fingerprints). Recent studies indicate that visible-near-infrared (VNIR) and shortwave-infrared (SWIR) reflectance characteristics of soil may be a rapid, inexpensive alternative to traditional fingerprint properties (e.g. geochemistry or mineral magnetism).
To further explore the applicability of VNIR-SWIR spectral data for sediment tracing purposes, source samples were collected in the Isabena watershed, a 445 km(2) dryland catchment in the central Spanish Pyrenees. Grab samples of the upper soil layer were collected from the main potential sediment source types along with in situ reflectance spectra. Samples were dried and sieved, and artificial mixtures of known proportions were produced for algorithm validation. Then, spectral readings of potential source and artificial mixture samples were taken in the laboratory. Colour coefficients and physically based parameters were calculated from in situ and laboratory-measured spectra. All parameters passing a number of prerequisite tests were subsequently applied in discriminant function analysis for source discrimination and mixing model analyses for source contribution assessment.
The three source types (i.e. badlands, forest/grassland and an aggregation of other sources, including agricultural land, shrubland, unpaved roads and open slopes) could be reliably identified based on spectral parameters. Laboratory-measured spectral fingerprints permitted the quantification of source contribution to artificial mixtures, and introduction of source heterogeneity into the mixing model decreased accuracies for some source types. Aggregation of source types that could not be discriminated did not improve mixing model results. Despite providing similar discrimination accuracies as laboratory source parameters, in situ derived source information was found to be insufficient for contribution modelling.
The laboratory mixture experiment provides valuable insights into the capabilities and limitations of spectral fingerprint properties. From this study, we conclude that combinations of spectral properties can be used for mixing model analyses of a restricted number of source groups, whereas more straightforward in situ measured source parameters do not seem suitable. However, modelling results based on laboratory parameters also need to be interpreted with care and should not rely on the estimates of mean values only but should consider uncertainty intervals as well.
Many Mediterranean drylands are characterized by strong erosion in headwater catchments, where connectivity processes play an important role in the redistribution of water and sediments. Sediment connectivity describes the ease with which sediment can move through a catchment. The spatial and temporal characterization of connectivity patterns in a catchment enables the estimation of sediment contribution and transfer paths. Apart from topography, vegetation cover is one of the main factors driving sediment connectivity. This is particularly true for the patchy vegetation cover typical of many dryland environments. Several connectivity measures have been developed in the last few years. At the same time, advances in remote sensing have enabled an improved catchment-wide estimation of ground cover at the subpixel level using hyperspectral imagery.
The objective of this study was to assess the sediment connectivity for two adjacent subcatchments (similar to 70 km(2)) of the Isabena River in the Spanish Pyrenees in contrasting seasons using a quantitative connectivity index based on fractional vegetation cover and topography data. The fractional cover of green vegetation, non-photosynthetic vegetation, bare soil and rock were derived by applying a multiple endmember spectral mixture analysis approach to the hyperspectral image data. Sediment connectivity was mapped using the index of connectivity, in which the effect of land cover on runoff and sediment fluxes is expressed by a spatially distributed weighting factor. In this study, the cover and management factor (C factor) of the Revised Universal Soil Loss Equation (RUSLE) was used as a weighting factor. Bi-temporal C factor maps were derived by linking the spatially explicit fractional ground cover and vegetation height obtained from the airborne data to the variables of the RUSLE subfactors.
The resulting connectivity maps show that areas behave very differently with regard to connectivity, depending on the land cover and on the spatial distribution of vegetation abundances and topographic barriers. Most parts of the catchment show higher connectivity values in August as compared to April. The two subcatchments show a slightly different connectivity behaviour that reflects the different land cover proportions and their spatial configuration.
The connectivity estimation can support a better understanding of processes controlling the redistribution of water and sediments from the hillslopes to the channel network at a scale appropriate for land management. It allows hot spot areas of erosion to be identified and the effects of erosion control measures, as well as different land management scenarios, to be studied.
Response spectra are of fundamental importance in earthquake engineering and represent a standard measure in seismic design for the assessment of structural performance. However, unlike Fourier spectral amplitudes, the relationship of response spectral amplitudes to seismological source, path, and site characteristics is not immediately obvious and might even be considered counterintuitive for high oscillator frequencies. The understanding of this relationship is nevertheless important for seismic-hazard analysis. The purpose of the present study is to comprehensively characterize the variation of response spectral amplitudes due to perturbations of the causative seismological parameters. This is done by calculating the absolute parameter sensitivities (sensitivity coefficients) defined as the partial derivatives of the model output with respect to its input parameters. To derive sensitivities, we apply algorithmic differentiation (AD). This powerful approach is extensively used for sensitivity analysis of complex models in meteorology or aerodynamics. To the best of our knowledge, AD has not been explored yet in the seismic-hazard context. Within the present study, AD was successfully implemented for a proven and extensively applied simulation program for response spectra (Stochastic Method SIMulation [SMSIM]) using the TAPENADE AD tool. We assess the effects and importance of input parameter perturbations on the shape of response spectra for different regional stochastic models in a quantitative way. Additionally, we perform sensitivity analysis regarding adjustment issues of groundmotion prediction equations.
Military areas are valuable habitats and refuges for rare and endangered plants and animals. We developed a new approach applying innovative methods of hyperspectral remote sensing to bridge the existing gap between remote sensing technology and the demands of the nature conservation community. Remote sensing has already proven to be a valuable monitoring instrument. However, the approaches lack the consideration of the demands of applied nature conservation which includes the legal demands of the EU Habitat Directive. Following the idea of the Vital Signs Monitoring in the USA, we identified a subset of the highest priority monitoring indicators for our study area. We analyzed continuous spectral response curves and tested the measurability of N=19 indicators on the basis of complexity levels aggregated from extensive vegetation assemblages. The spectral differentiability for the floristic as well as faunistic indicators revealed values up to 100% accuracy. We point out difficulties when it comes to distinguishing faunistic habitat requirements of several species adapted to dry open landscapes, which in this case results in OVERALL ACCURACY of 67, 87-95, and 35% in the error matrix. In summary, we provide an applicable and feasible method to facilitating monitoring military areas by hyperspectral remote sensing in the following. (C) 2014 Elsevier Ltd. All rights reserved.
Pollen records from large lakes have been used for quantitative palaeoclimate reconstruction, but the influences that lake size (as a result of species-specific variations in pollen dispersal patterns that smaller pollen grains are more easily transported to lake centre) and taphonomy have on these climatic signals have not previously been systematically investigated. We introduce the concept of pollen source area to pollen-based climate calibration using the north-eastern Tibetan Plateau as our study area. We present a pollen data set collected from large lakes in the arid to semi-arid region of central Asia. The influences that lake size and the inferred pollen source areas have on pollen compositions have been investigated through comparisons with pollen assemblages in neighbouring lakes of various sizes. Modern pollen samples collected from different parts of Lake Donggi Cona (in the north-eastern part of the Tibetan Plateau) reveal variations in pollen assemblages within this large lake, which are interpreted in terms of the species-specific dispersal and depositional patterns for different types of pollen, and in terms of fluvial input components. We have estimated the pollen source area for each lake individually and used this information to infer modern climate data with which to then develop a modern calibration data set, using both the multivariate regression tree (MRT) and weighted-averaging partial least squares (WA-PLS) approaches. Fossil pollen data from Lake Donggi Cona have been used to reconstruct the climate history of the north-eastern part of the Tibetan Plateau since the Last Glacial Maximum (LGM). The meanannual precipitation was quantitatively reconstructed using WA-PLS: extremely dry conditions are found to have dominated the LGM, with annual precipitation of around 100 mm, which is only 32% of present-day precipitation. A gradually increasing trend in moisture conditions during the Late Glacial is terminated by an abrupt reversion to a dry phase that lasts for about 1000 yr and coincides with "Heinrich event 1" in the North Atlantic region. Subsequent periods corresponding to the Bolling/Allerod interstadial, with annual precipitation (P-ann) of about 350 mm, and the Younger Dryas event (about 270 mm P-ann) are followed by moist conditions in the early Holocene, with annual precipitation of up to 400 mm. A drier trend after 9 cal. ka BP is followed by a second wet phase in the middle Holocene, lasting until 4.5 cal. ka BP. Relatively steady conditions with only slight fluctuations then dominate the late Holocene, resulting in the present climatic conditions. The climate changes since the LGM have been primarily driven by deglaciation and fluctuations in the intensity of the Asian summer monsoon that resulted from changes in the Northern Hemisphere summer solar insolation, as well as from changes in the North Atlantic climate through variations in the circulation patterns and intensity of the westerlies.
We present pollen-based reconstructions of the spatio-temporal dynamics of northern European regional vegetation abundance through the Holocene. We apply the Regional Estimates of VEgetation Abundance from Large Sites (REVEALS) model using fossil pollen records from eighteen sites within five modern biomes in the region. The eighteen sites are classified into four time-trajectory types on the basis of principal components analysis of both the REVEALS-based vegetation estimates (RVs) and the pollen percentage (PPs). The four trajectory types are more clearly separated for RVs than PPs. Further, the timing of major Holocene shifts, rates of compositional change, and diversity indices (turnover and evenness) differ between RVs and PPs. The differences are due to the reduction by REVEALS of biases in fossil pollen assemblages caused by different basin size, and inter-taxonomic differences in pollen productivity and dispersal properties. For example, in comparison to the PPs, the RVs show an earlier increase in Corylus and Ulmus in the early-Holocene and a more pronounced increase in grassland and deforested areas since the mid-Holocene. The results suggest that the influence of deforestation and agricultural activities on plant composition and abundance from Neolithic times was stronger than previously inferred from PPs. Relative to PPs, RVs show a more rapid compositional change, a largest decrease in turnover, and less variable evenness in most of northern Europe since 5200 cal yr BP. All these changes are primarily related to the strong impact of human activities on the vegetation. This study demonstrates that RV-based estimates of diversity indices, timing of shifts, and rates of change in reconstructed vegetation provide new insights into the timing and magnitude of major human distribution on Holocene regional, vegetation, feature that are critical in the assessment of human impact on vegetation, land-cover, biodiversity, and climate in the past.
The Tibetan Plateau is a region that is highly sensitive to recent global warming, but the complexity and heterogeneity of its mountainous landscape can result in variable responses. In addition, the scarcity and brevity of regional instrumental and palaeoecological records still hamper our understanding of past and present patterns of environmental change. To investigate how the remote, high-alpine environments of the Nianbaoyeze Mountains, eastern Tibetan Plateau, are affected by climate change and human activity over the last similar to 600 years, we compared regional tree-ring studies with pollen and diatom remains archived in the dated sediments of Dongerwuka Lake (33.22A degrees N, 101.12A degrees E, 4,307 m a.s.l.). In agreement with previous studies from the eastern Tibetan Plateau, a strong coherence between our two juniper-based tree-ring chronologies from the Nianbaoyeze and the Anemaqin Mountains was observed, with pronounced cyclical variations in summer temperature reconstructions. A positive directional trend to warmer summer temperatures in the most recent decades, was, however, not observed in the tree-ring record. Likewise, our pollen and diatom spectra showed minimal change over the investigated time period. Although modest, the most notable change in the diatom relative abundances was a subtle decrease in the dominant planktonic Cyclotella ocellata and a concurrent increase in small, benthic fragilarioid taxa in the similar to 1820s, suggesting higher ecosystem variability. The pollen record subtly indicates three periods of increased cattle grazing activity (similar to 1400-1480 AD, similar to 1630-1760 AD, after 1850 AD), but shows generally no significant vegetation changes during past similar to 600 years. The minimal changes observed in the tree-ring, diatom and pollen records are consistent with the presence of localised cooling centres that are evident in instrumental and tree-ring data within the southeastern and eastern Tibetan Plateau. Given the minor changes in regional temperature records, our complacent palaeoecological profiles suggest that climatically induced ecological thresholds have not yet been crossed in the Nianbaoyeze Mountains region.
(Paleo-)climatologists are challenged to identify mechanisms that cause the observed abrupt Holocene monsoon events despite the fact that monsoonal circulation is assumed to be driven by gradual insolation changes. Here we provide proxy and model evidence to show that moisture-advection feedback can lead to a non-linear relationship between sea-surface and continental temperatures and monsoonal precipitation. A pollen record from Lake Ximencuo (Nianbaoyeze Mountains) indicates that vegetation from the eastern margin of the Tibetan Plateau was characterized by alpine deserts and glacial flora after the Last Glacial Maximum (LGM) (21-15.5 cal kyr BP), by alpine meadows during the Late Glacial (15.5-10.4 cal kyr BP) and second half of the Holocene (5.0 cal kyr BP to present) and by mixed forests during the first half of the Holocene (10.4-5.0 cal kyr BP). The application of pollen-based transfer functions yields an abrupt temperature increase at 10.4 cal kyr BP and a decrease at 5.0 cal kyr BP of about 3 degrees C. By applying endmember modeling to grain-size data from the same sediment core we infer that frequent fluvial events (probably originating from high-magnitude precipitation events) were more common in the early and mid Holocene. We assign the inferred exceptional strong monsoonal circulation to the initiation of moisture-advection feedback, a result supported by a simple model that reproduces this feedback pattern over the same time period. (C) 2014 Published by Elsevier B.V.
Lacustrine sediments have been widely used to investigate past climatic and environmental changes on millennial to seasonal time scales. Sedimentary archives of lakes in mountainous regions may also record non-climatic events such as earthquakes. We argue herein that a set of 64 annual laminae couplets reconciles a stratigraphically inconsistent accelerator mass spectrometry (AMS) C-14 chronology in a similar to 4-m-long sediment core from Lake Mengda, in the north-eastern Tibetan Plateau. The laminations suggest the lake was formed by a large landslide, triggered by the 1927 Gulang earthquake (M = 8.0). The lake sediment sequence can be separated into three units based on lithologic, sedimentary, and isotopic characteristics. Starting from the bottom of the sequence, these are: (1) unweathered, coarse, sandy valley-floor deposits or landslide debris that pre-date the lake, (2) landslide-induced, fine-grained soil or reworked landslide debris with a high organic content, and (3) lacustrine sediments with low organic content and laminations. These annual laminations provide a high-resolution record of anthropogenic and environmental changes during the twentieth century, recording enhanced sediment input associated with two phases of construction activities. The high mean sedimentation rates of up to 4.8 mm year(-1) underscore the potential for reconstructing such distinct sediment pulses in remote, forested, and seemingly undisturbed mountain catchments.
The widely reported paradox of methane oversaturation in oxygenated water challenges the prevailing paradigm that microbial methanogenesis only occurs under anoxic conditions. Using a combination of field sampling, incubation experiments, and modeling, we show that the recurring mid-water methane peak in Lake Stechlin, northeast Germany, was not dependent on methane input from the littoral zone or bottom sediment or on the presence of known micro-anoxic zones. The methane peak repeatedly overlapped with oxygen oversaturation in the seasonal thermocline. Incubation experiments and isotope analysis indicated active methane production, which was likely linked to photosynthesis and/or nitrogen fixation within the oxygenated water, whereas lessening of methane oxidation by light allowed accumulation of methane in the oxygen-rich upper layer. Estimated methane efflux from the surface water was up to 5 mmol m(-2) d(-1). Mid-water methane oversaturation was also observed in nine other lakes that collectively showed a strongly negative gradient of methane concentration within 0-20% dissolved oxygen (DO) in the bottom water, and a positive gradient within >= 20% DO in the upper water column. Further investigation into the responsible organisms and biochemical pathways will help improve our understanding of the global methane cycle.
A modern pollen dataset from China and Mongolia (18-52 degrees N, 74-132 degrees E) is investigated for its potential use in climate reconstructions. The dataset includes 2559 samples, 229 terrestrial pollen taxa and four climatic variables - mean annual precipitation (P-ann): 35-2091 mm, mean annual temperature (T-ann): -12.1-25.8 degrees C, mean temperature in the coldest month (Mt(co).): -33.8-21.7 degrees C, and mean temperature in the warmest month (Mt(wa)): 03-29.8 degrees C. Modern pollen-climate relationships are assessed using canonical correspondence analysis (CCA), Huisman-Olff-Fresco (HOF) models, the modern analogue technique (MAT), and weighted averaging partial least squares (WA-PLS). Results indicate that P-ann is the most important climatic determinant of pollen distribution and the most promising climate variable for reconstructions, as assessed by the coefficient of determination between observed and predicted environmental values (r(2)) and root mean square error of prediction (RMSEP). Mt(co) and Mt(wa) may be reconstructed too, but with caution. Samples from different depositional environments influence the performance of cross-validation differently, with samples from lake sediment-surfaces and moss polsters having the best fit with the lowest RMSEP. The better model performances of MAT are most probably caused by spatial autocorrelation. Accordingly, the WA-PLS models of this dataset are deemed most suitable for reconstructing past climate quantitatively because of their more reliable predictive power. (C) 2014 Elsevier B.V. All rights reserved.
AimFossil pollen spectra from lake sediments in central and western Mongolia have been used to interpret past climatic variations, but hitherto no suitable modern pollen-climate calibration set has been available to infer past climate changes quantitatively. We established such a modern pollen dataset and used it to develop a transfer function model that we applied to a fossil pollen record in order to investigate: (1) whether there was a significant moisture response to the Younger Dryas event in north-western Mongolia; and (2) whether the early Holocene was characterized by dry or wet climatic conditions.
LocationCentral and western Mongolia.
MethodsWe analysed pollen data from surface sediments from 90 lakes. A transfer function for mean annual precipitation (P-ann) was developed with weighted averaging partial least squares regression (WA-PLS) and applied to a fossil pollen record from Lake Bayan Nuur (49.98 degrees N, 93.95 degrees E, 932m a.s.l.). Statistical approaches were used to investigate the modern pollen-climate relationships and assess model performance and reconstruction output.
ResultsRedundancy analysis shows that the modern pollen spectra are characteristic of their respective vegetation types and local climate. Spatial autocorrelation and significance tests of environmental variables show that the WA-PLS model for P-ann is the most valid function for our dataset, and possesses the lowest root mean squared error of prediction.
Main conclusionsPrecipitation is the most important predictor of pollen and vegetation distributions in our study area. Our quantitative climate reconstruction indicates a dry Younger Dryas, a relatively dry early Holocene, a wet mid-Holocene and a dry late Holocene.
This study examines the course and driving forces of recent vegetation change in the Mongolian steppe. A sediment core covering the last 55years from a small closed-basin lake in central Mongolia was analyzed for its multi-proxy record at annual resolution. Pollen analysis shows that highest abundances of planted Poaceae and highest vegetation diversity occurred during 1977-1992, reflecting agricultural development in the lake area. A decrease in diversity and an increase in Artemisia abundance after 1992 indicate enhanced vegetation degradation in recent times, most probably because of overgrazing and farmland abandonment. Human impact is the main factor for the vegetation degradation within the past decades as revealed by a series of redundancy analyses, while climate change and soil erosion play subordinate roles. High Pediastrum (a green algae) influx, high atomic total organic carbon/total nitrogen (TOC/TN) ratios, abundant coarse detrital grains, and the decrease of C-13(org) and N-15 since about 1977 but particularly after 1992 indicate that abundant terrestrial organic matter and nutrients were transported into the lake and caused lake eutrophication, presumably because of intensified land use. Thus, we infer that the transition to a market economy in Mongolia since the early 1990s not only caused dramatic vegetation degradation but also affected the lake ecosystem through anthropogenic changes in the catchment area.
Subfossil Cladocera were sampled and examined from the surface sediments of 35 thermokarst lakes along a temperature gradient crossing the tree line in the Anabar-river basin in northwestern Yakutia, northeastern Siberia. The lakes were distributed through three environmental zones: typical tundra, southern tundra and forest tundra. All lakes were situated within the continuous permafrost zone. Our investigation showed that the cladoceran communities in the lakes of the Anabar region are diverse and abundant, as reflected by taxonomic richness, and high diversity and evenness indices (H = 1.89 +/- A 0.51; I = 0.8 +/- A 0.18). CONISS cluster analysis indicated that the cladoceran communities in the three ecological zones (typical tundra, southern tundra and forest-tundra) differed in their taxonomic composition and structure. Differences in the cladoceran assemblages were related to limnological features and geographical position, vegetation type, climate and water chemistry. The constrained redundancy analysis indicated that T-July, water depth and both sulphate (SO4 (2-)) and silica (Si4+) concentrations significantly (p a parts per thousand currency sign 0.05) explained variance in the cladoceran assemblage. T-July featured the highest percentage (17.4 %) of explained variance in the distribution of subfossil Cladocera. One of the most significant changes in the structure of the cladoceran communities in the investigated transect was the replacement of closely related species along the latitudinal and vegetation gradient. The results demonstrate the potential for a regional cladoceran-based temperature model for the Arctic regions of Russia, and for and Yakutia in particular.
Late Holocene glacier variations in westernmost Tibetan Plateau were studied based on the analysis of grain size, magnetic susceptibility, and elements from an 8.3m long distal glaciolacustrine sediment core of Kalakuli Lake. Our results show that there are four glacier expansion episodes occurring in 4200-3700calibrated years (cal years) B.P., 2950-2300cal years B.P., 1700-1070cal years B.P., and 570-100cal years B.P. and four glacier retreat periods of 3700-2950cal years B.P., 2300-1700cal years B.P., 1070-570cal years B.P., and 50cal years B.P.-present. The four glacier expansion episodes are generally in agreement with the glacier activities indicted by the moraines at Muztagh Ata and Kongur Shan, as well as with the late Holocene ice-rafting events in the North Atlantic. Over the last 2000years, our reconstructed glacier variations are in temporal agreement with reconstructed temperature from China and the Northern Hemisphere, indicating that glacier variations at centennial time scales are very sensitive to temperature in western Tibetan Plateau.
For a detailed characterization of near-surface environments, geophysical techniques are increasingly used to support more conventional point-based techniques such as borehole and direct-push logging. Because the underlying parameter relations are often complex, site-specific, or even poorly understood, a remaining challenging task is to link the geophysical parameter models to the actual geotechnical target parameters measured only at selected points. We propose a workflow based on nonparametric regression to establish functional relationships between jointly inverted geophysical parameters and selected geotechnical parameters as measured, for example, by different borehole and direct-push tools. To illustrate our workflow, we present field data collected to characterize a near-surface sedimentary environment Our field data base includes crosshole ground penetrating radar (GPR), seismic P-, and S-wave data sets collected between 25 m deep boreholes penetrating sand- and gravel dominated sediments. Furthermore, different typical borehole and direct-push logs are available. We perform a global joint inversion of traveltimes extracted from the crosshole geophysical data using a recently proposed approach based on particle swarm optimization. Our inversion strategy allows for generating consistent models of GPR, P-wave, and S-wave velocities including an appraisal of uncertainties. We analyze the observed complex relationships between geophysical velocities and target parameter logs using the alternating conditional expectation (ACE) algorithm. This nonparametric statistical tool allows us to perform multivariate regression analysis without assuming a specific functional relation between the variables. We are able to explain selected target parameters such as characteristic grain size values or natural gamma activity by our inverted geophysical data and to extrapolate these parameters to the inter-borehole plane covered by our crosshole experiments. We conclude that the ACE algorithm is a powerful tool to analyze a multivariate petrophysical data base and to develop an understanding of how a multi-parameter geophysical model can be linked and translated to selected geotechnical parameters.
Velocity models are essential to process two-and three-dimensional ground-penetrating radar (GPR) data. Furthermore, velocity information aids the interpretation of such data sets because velocity variations reflect important material properties such as water content. In many GPR applications, common midpoint (CMP) surveys are routinely collected to determine one-dimensional velocity models at selected locations. To analyse CMP data gathers, spectral velocity analyses relying on the normal-moveout (NMO) model are commonly employed. Using Dix's formula, the derived NMO velocities can be further converted to interval velocities which are needed for processing and interpretation. Because of the inherent assumptions and limitations of such approaches, we investigate and propose an alternative procedure based on the global inversion of reflection travel-times. We use a finite-difference solver of the Eikonal equation to accurately solve the forward problem in combination with particle swarm optimization (PSO) to find one-dimensional GPR velocity models explaining our data. Because PSO is a robust and efficient global optimization tool, our inversion approach includes generating an ensemble of representative solutions that allows us to analyse uncertainties in the model space. Using synthetic data examples, we test and evaluate our inversion approach to analyse CMP data collected across typical near-surface environments. Application to a field data set recorded at a well-constrained test site including a comparison to independent borehole and direct-push data, further illustrates the potential of the proposed approach, which includes a straightforward and understandable appraisal of non-uniqueness and uncertainty issues, respectively. We conclude that our methodology is a feasible and powerful tool to analyse GPR CMP data and allows practitioners and researchers to evaluate the reliability of CMP derived velocity models.
Vertical radar profiling (VRP) is a single-borehole geophysical technique, in which the receiver antenna is located within a borehole and the transmitter antenna is placed at one or various offsets from the borehole. Today, VRP surveying is primarily used to derive 1D velocity models by inverting the arrival times of direct waves. Using field data collected at a well-constrained test site in Germany, we evaluated a VRP workflow relying on the analysis of direct-arrival traveltimes and amplitudes as well as on imaging reflection events. To invert our VRP traveltime data, we used a global inversion strategy resulting in an ensemble of acceptable velocity models, and thus, it allowed us to appraise uncertainty issues in the estimated velocities as well as in porosity models derived via petrophysical translations. In addition to traveltime inversion, the analysis of direct-wave amplitudes and reflection events provided further valuable information regarding subsurface properties and architecture. The used VRP amplitude preprocessing and inversion procedures were adapted from raybased crosshole ground-penetrating radar (GPR) attenuation tomography and resulted in an attenuation model, which can be used to estimate variations in electrical resistivity. Our VRP reflection imaging approach relied on corridor stacking, which is a well-established processing sequence in vertical seismic profiling. The resulting reflection image outlines bounding layers and can be directly compared to surface-based GPR reflection profiling. Our results of the combined analysis of VRP, traveltimes, amplitudes, and reflections were consistent with independent core and borehole logs as well as GPR reflection profiles, which enabled us to derive a detailed hydro-stratigraphic model as needed, for example, to understand and model groundwater flow and transport.
Geophysical techniques offer the potential to tomographically image physical parameter variations in the ground in two or three dimensions. Due to the limited number and accuracy of the recorded data, geophysical model generation by inversion suffers ambiguity. Linking the model generation process of disparate data by jointly inverting two or more data sets allows for improved model reconstruction. Fully nonlinear inversion using optimization techniques searching the solution space of the inverse problem globally enables quantitative assessment of the ambiguity inherent to the model reconstruction. We used two different multiobjective particle swarm optimization approaches to jointly invert synthetic crosshole tomographic data sets comprising radar and P-wave traveltimes, respectively. Beginning with a nonlinear joint inversion founded on the principle of Pareto optimality and game theoretic concepts, we obtained a set of Pareto-optimal solutions comprising commonly structured radar and P-wave velocity models for low computational costs. However, the efficiency of the approach goes along with some risk of achieving a final model ensemble not adequately illustrating the ambiguity inherent to the model reconstruction process. Taking advantage of the results of the first approach, we inverted the database using a different nonlinear joint-inversion approach reducing the multiobjective optimization problem to a single-objective one. Computational costs were significantly higher, but the final models were obtained mutually independently allowing for objective appraisal of model parameter determination. Despite the high computational effort, the approach was found to be an efficient nonlinear joint-inversion formulation compared to what could be extracted from individual nonlinear inversions of both data sets.
According to Dooge (1986) intermediate-scale catchments are systems of organized complexity, being too organized and yet too small to be characterized on a statistical/conceptual basis, but too large and too heterogeneous to be characterized in a deterministic manner. A key requirement for building structurally adequate models precisely for this intermediate scale is a better understanding of how different forms of spatial organization affect storage and release of water and energy. Here, we propose that a combination of the concept of hydrological response units (HRUs) and thermodynamics offers several helpful and partly novel perspectives for gaining this improved understanding. Our key idea is to define functional similarity based on similarity of the terrestrial controls of gradients and resistance terms controlling the land surface energy balance, rainfall runoff transformation, and groundwater storage and release. This might imply that functional similarity with respect to these specific forms of water release emerges at different scales, namely the small field scale, the hillslope, and the catchment scale. We thus propose three different types of "functional units" - specialized HRUs, so to speak - which behave similarly with respect to one specific form of water release and with a characteristic extent equal to one of those three scale levels. We furthermore discuss an experimental strategy based on exemplary learning and replicate experiments to identify and delineate these functional units, and as a promising strategy for characterizing the interplay and organization of water and energy fluxes across scales. We believe the thermodynamic perspective to be well suited to unmask equifinality as inherent in the equations governing water, momentum, and energy fluxes: this is because several combinations of gradients and resistance terms yield the same mass or energy flux and the terrestrial controls of gradients and resistance terms are largely independent. We propose that structurally adequate models at this scale should consequently disentangle driving gradients and resistance terms, because this optionally allow sequifinality to be partly reduced by including available observations, e. g., on driving gradients. Most importantly, the thermodynamic perspective yields an energy-centered perspective on rainfall-runoff transformation and evapotranspiration, including fundamental limits for energy fluxes associated with these processes. This might additionally reduce equifinality and opens up opportunities for testing thermodynamic optimality principles within independent predictions of rainfall-runoff or land surface energy exchange. This is pivotal to finding out whether or not spatial organization in catchments is in accordance with a fundamental organizing principle.
To assess the ecological and economic implications of the redistributive land reform in semi-arid Namibia, we investigated to what extent land reform beneficiaries adjust herd size and herd composition according to environmental (rainfall, vegetation) and economic variables (herd size, financial assets, running costs). We performed model-based role-plays with Namibian land reform beneficiaries, simulating 10 years of rangeland management.
Our study revealed that the farmers surveyed mainly manage their herds according to their economic situation (herd size and account balance) but do not take environmental variability (rainfall and vegetation) into account. Further, our results indicate that, due to financial pressure, farmers are not able to apply their desired management strategies, and that owners of small farms face a higher risk of economic failure. However, farmers apply rather conservative and constant stocking rates and will thus, given the current economic limitations, likely not contribute to semi-arid savanna degradation.
We conclude that land reform beneficiaries need support to be able to apply straightforward and efficient management strategies. This could be achieved by facilitating cooperation between small farming businesses and by supporting initial investment in productive cattle herds at the time of redistribution of the land.
In Lutzito catchment on Barro Colorado Island, Panama, extraordinarily high suspended-sediment yields of 1-2Mgha-1year-1 were generated despite the dense forest cover coinciding with erosion-resistant soils. We hypothesized that ant mounding activity is an important zoogeomorphological mechanism in this area, providing relevant quantities of easily transportable material at the soil surface. To test this hypothesis, all ant mound material was collected collected for dry mass determination from thirty 4m2 plots installed in the study area every 1-3days during the 39-day sampling period. Additionally, three ground-nesting ant species responsible for mounds in the study area, Ectatomma ruidum, Trachymyrmex cornetzi and Strumigenys marginiventris, were identified. On the basis of the total of 1.38kg of material collected in the wet season of 2011, the estimate for the whole 8months wet season amounts to 725kgha-1. As this value is in the same order of magnitude as sediment output, it shows that ants may act as important ecosystem engineers and contribute to sediment production here by providing large quantities of fine-grained, readily erodible material at the soil surface for subsequent transport to the streambed. Copyright (c) 2014 John Wiley & Sons, Ltd.