Refine
Year of publication
- 2014 (226) (remove)
Document Type
- Article (196)
- Doctoral Thesis (17)
- Preprint (8)
- Review (4)
- Postprint (1)
Is part of the Bibliography
- yes (226) (remove)
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 (226) (remove)
Surface displacement at volcanic edifices is related to subsurface processes associated with magma movements, fluid transfers within the volcano edifice and gravity-driven deformation processes. Understanding of associated ground displacements is of importance for assessment of volcanic hazards. For example, volcanic unrest is often preceded by surface uplift, caused by magma intrusion and followed by subsidence, after the withdrawal of magma. Continuous monitoring of the surface displacement at volcanoes therefore might allow the forecasting of upcoming eruptions to some extent. In geophysics, the measured surface displacements allow the parameters of possible deformation sources to be estimated through analytical or numerical modeling. This is one way to improve the understanding of subsurface processes acting at volcanoes. Although the monitoring of volcanoes has significantly improved in the last decades (in terms of technical advancements and number of monitored volcanoes), the forecasting of volcanic eruptions remains puzzling. In this work I contribute towards the understanding of the subsurface processes at volcanoes and thus to the improvement of volcano eruption forecasting. I have investigated the displacement field of Llaima volcano in Chile and of Tendürek volcano in East Turkey by using synthetic aperture radar interferometry (InSAR). Through modeling of the deformation sources with the extracted displacement data, it was possible to gain insights into potential subsurface processes occurring at these two volcanoes that had been barely studied before. The two volcanoes, although of very different origin, composition and geometry, both show a complexity of interacting deformation sources. At Llaima volcano, the InSAR technique was difficult to apply, due to the large decorrelation of the radar signal between the acquisition of images. I developed a model-based unwrapping scheme, which allows the production of reliable displacement maps at the volcano that I used for deformation source modeling. The modeling results show significant differences in pre- and post-eruptive magmatic deformation source parameters. Therefore, I conjecture that two magma chambers exist below Llaima volcano: a post-eruptive deep one and a shallow one possibly due to the pre-eruptive ascent of magma. Similar reservoir depths at Llaima have been confirmed by independent petrologic studies. These reservoirs are interpreted to be temporally coupled. At Tendürek volcano I have found long-term subsidence of the volcanic edifice, which can be described by a large, magmatic, sill-like source that is subject to cooling contraction. The displacement data in conjunction with high-resolution optical images, however, reveal arcuate fractures at the eastern and western flank of the volcano. These are most likely the surface expressions of concentric ring-faults around the volcanic edifice that show low magnitudes of slip over a long time. This might be an alternative mechanism for the development of large caldera structures, which are so far assumed to be generated during large catastrophic collapse events. To investigate the potential subsurface geometry and relation of the two proposed interacting sources at Tendürek, a sill-like magmatic source and ring-faults, I have performed a more sophisticated numerical modeling approach. The optimum source geometries show, that the size of the sill-like source was overestimated in the simple models and that it is difficult to determine the dip angle of the ring-faults with surface displacement data only. However, considering physical and geological criteria a combination of outward-dipping reverse faults in the west and inward-dipping normal faults in the east seem to be the most likely. Consequently, the underground structure at the Tendürek volcano consists of a small, sill-like, contracting, magmatic source below the western summit crater that causes a trapdoor-like faulting along the ring-faults around the volcanic edifice. Therefore, the magmatic source and the ring-faults are also interpreted to be temporally coupled. In addition, a method for data reduction has been improved. The modeling of subsurface deformation sources requires only a relatively small number of well distributed InSAR observations at the earth’s surface. Satellite radar images, however, consist of several millions of these observations. Therefore, the large amount of data needs to be reduced by several orders of magnitude for source modeling, to save computation time and increase model flexibility. I have introduced a model-based subsampling approach in particular for heterogeneously-distributed observations. It allows a fast calculation of the data error variance-covariance matrix, also supports the modeling of time dependent displacement data and is, therefore, an alternative to existing methods.
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.
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.
The purpose of this thesis is to develop an automated inversion scheme to derive point and finite source parameters for weak earthquakes, here intended with the unusual meaning of earthquakes with magnitudes at the limit or below the bottom magnitude threshold of standard source inversion routines. The adopted inversion approaches entirely rely on existing inversion software, the methodological work mostly targeting the development and tuning of optimized inversion flows. The resulting inversion scheme is tested for very different datasets, and thus allows the discussion on the source inversion problem at different scales. In the first application, dealing with mining induced seismicity, the source parameters determination is addressed at a local scale, with source-sensor distance of less than 3 km. In this context, weak seismicity corresponds to event below magnitude MW 2.0, which are rarely target of automated source inversion routines. The second application considers a regional dataset, namely the aftershock sequence of the 2010 Maule earthquake (Chile), using broadband stations at regional distances, below 300 km. In this case, the magnitude range of the target aftershocks range down to MW 4.0. This dataset is here considered as a weak seismicity case, since the analysis of such moderate seismicity is generally investigated only by moment tensor inversion routines, with no attempt to resolve source duration or finite source parameters. In this work, automated multi-step inversion schemes are applied to both datasets with the aim of resolving point source parameters, both using double couple (DC) and full moment tensor (MT) models, source duration and finite source parameters. A major result of the analysis of weaker events is the increased size of resulting moment tensor catalogues, which interpretation may become not trivial. For this reason, a novel focal mechanism clustering approach is used to automatically classify focal mechanisms, allowing the investigation of the most relevant and repetitive rupture features. The inversion of the mining induced seismicity dataset reveals the repetitive occurrence of similar rupture processes, where the source geometry is controlled by the shape of the mined panel. Moreover, moment tensor solutions indicate a significant contribution of tensile processes. Also the second application highlights some characteristic geometrical features of the fault planes, which show a general consistency with the orientation of the slab. The additional inversion for source duration allowed to verify the empirical correlation for moment normalized earthquakes in subduction zones among a decreasing rupture duration with increasing source depth, which was so far only observed for larger events.
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.
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 characterised by low signal-to-noise ratio, such 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 thesis I developed a picking free automated method based on the Short-Term-Average/Long-Term-Average (STA/LTA) traces at different stations as observed data. I used the STA/LTA of several characteristic functions in order to increase the sensitiveness to the P wave and the S waves. For the P phases we use the STA/LTA traces of the vertical energy function, while for the S phases, we use the STA/LTA traces of the horizontal energy trace and then a more optimized characteristic function which is obtained using the principal component analysis technique. The orientation of the horizontal components can be retrieved by robust and linear approach of waveform comparison between stations within a network using seismic sources outside the network (chapter 2). 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 three-dimensional grid we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial and temporal coordinates of the seismic event. Location uncertainties are then estimated by perturbing the STA/LTA parameters (i.e the length of both long and short time windows) and relocating each event several times. In order to test the location method I firstly applied it to a set of 200 synthetic events. Then we applied it to two different real datasets. A first one related to mining induced microseismicity in a coal mine in the northern Germany (chapter 3). In this case we successfully located 391 microseismic event with magnitude range between 0.5 and 2.0 Ml. To further validate the location method I compared the retrieved locations with those obtained by manual picking procedure. The second dataset consist in a pilot application performed in the Campania-Lucania region (southern Italy) using a 33 stations seismic network (Irpinia Seismic Network) with an aperture of about 150 km (chapter 4). We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < Ml < 2.7. A subset of these locations were compared with accurate locations retrieved by a manual location procedure based on the use of a double difference technique. In both cases results indicate good agreement with manual locations. Moreover, the waveform stacking location method results noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals.
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.
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.