Refine
Year of publication
- 2014 (228) (remove)
Document Type
- Article (196)
- Doctoral Thesis (19)
- Preprint (8)
- Review (4)
- Postprint (1)
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 (228) (remove)
The 2002 M-w 7.9 Denali Fault earthquake, Alaska, provides an unparalleled opportunity to investigate in quantitative detail the regional hillslope mass-wasting response to strong seismic shaking in glacierized terrain. We present the first detailed inventory of similar to 1580 coseismic slope failures, out of which some 20% occurred above large valley glaciers, based on mapping from multi-temporal remote sensing data. We find that the Denali earthquake produced at least one order of magnitude fewer landslides in a much narrower corridor along the fault ruptures than empirical predictions for an M 8 earthquake would suggest, despite the availability of sufficiently steep and dissected mountainous topography prone to frequent slope failure. In order to explore potential controls on the reduced extent of regional coseismic landsliding we compare our data with inventories that we compiled for two recent earthquakes in periglacial and formerly glaciated terrain, i.e. at Yushu, Tibet (M-w 6.9, 2010), and Aysen Fjord, Chile (2007 M-w 6.2). Fault movement during these events was, similarly to that of the Denali earthquake, dominated by strike-slip offsets along near-vertical faults. Our comparison returns very similar coseismic landslide patterns that are consistent with the idea that fault type, geometry, and dynamic rupture process rather than widespread glacier cover were among the first-order controls on regional hillslope erosional response in these earthquakes. We conclude that estimating the amount of coseismic hillslope sediment input to the sediment cascade from earthquake magnitude alone remains highly problematic, particularly if glacierized terrain is involved. (C) 2014 Elsevier Ltd. All rights reserved.
Soil in a changing world is subject to both anthropogenic and environmental stresses. Soil monitoring is essential to assess the magnitude of changes in soil variables and how they affect ecosystem processes and human livelihoods. However, we cannot always be sure which sampling design is best for a given monitoring task. We employed a rotational stratified simple random sampling (rotStRS) for the estimation of temporal changes in the spatial mean of saturated hydraulic conductivity (K-s) at three sites in central Panama in 2009, 2010 and 2011. To assess this design's efficiency we compared the resulting estimates of the spatial mean and variance for 2009 with those gained from stratified simple random sampling (StRS), which was effectively the data obtained on the first sampling time, and with an equivalent unexecuted simple random sampling (SRS). The poor performance of geometrical stratification and the weak predictive relationship between measurements of successive years yielded no advantage of sampling designs more complex than SRS. The failure of stratification may be attributed to the small large-scale variability of K-s. Revisiting previously sampled locations was not beneficial because of the large small-scale variability in combination with destructive sampling, resulting in poor consistency between revisited samples. We conclude that for our K-s monitoring scheme, repeated SRS is equally effective as rotStRS. Some problems of small-scale variability might be overcome by collecting several samples at close range to reduce the effect of small-scale variation. Finally, we give recommendations on the key factors to consider when deciding whether to use stratification and rotation in a soil monitoring scheme.
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.
The Dansgaard-Oeschger oscillations and Heinrich events described in North Atlantic sediments and Greenland ice are expressed in the climate of the tropics, for example, as documented in Arabian Sea sediments. Given the strength of this teleconnection, we seek to reconstruct its range of environmental impacts. We present geochemical and sedimentological data from core SO130-289KL from the Indus submarine slope spanning the last similar to 80 kyr. Elemental and grain size analyses consistently indicate that interstadials are characterized by an increased contribution of fluvial suspension from the Indus River. In contrast, stadials are characterized by an increased contribution of aeolian dust from the Arabian Peninsula. Decadal-scale shifts at climate transitions, such as onsets of interstadials, were coeval with changes in productivity-related proxies. Heinrich events stand out as especially dry and dusty events, indicating a dramatically weakened Indian summer monsoon, potentially increased winter monsoon circulation, and increased aridity on the Arabian Peninsula. This finding is consistent with other paleoclimate evidence for continental aridity in the northern tropics during these events. Our results strengthen the evidence that circum-North Atlantic temperature variations translate to hydrological shifts in the tropics, with major impacts on regional environmental conditions such as rainfall, river discharge, aeolian dust transport, and ocean margin anoxia.
Leaching of dissolved C in arable hummocky ground moraine soil landscapes is characterized by a spatial continuum of more or less erosion-affected Luvisols, Calcaric Regosols at exposed positions, and Colluvic Regosols in depressions. Our objective was to estimate the fluxes of dissolved C in four differently eroded soils as affected by erosion-induced pedological and soil structural alterations. In this model study, we considered landscape position effects by adapting the water table as the bottom boundary condition and erosion effects by using pedon-specific soil hydraulic properties. The one-dimensional vertical water movement was described with the Richards equation using HYDRUS-1D. Solute fluxes were obtained by combining calculated water fluxes with concentrations of dissolved organic and inorganic C (DOC and DIC, respectively) measured from soil solution extracted by suction cups at biweekly intervals. In the 3-yr period (2010-2012), DOC fluxes in the 2-m soil depth were similar at the three non-colluvic locations with -0.8 +/- 0.1 g m(-2) yr(-1) (i.e., outflow) but were 0.4 g m(-2) yr(-1) (i.e., input) in the depression. The DIC fluxes ranged from -10.2 g m(-2) yr(-1) for the eroded Luvisol, -9.2 g m(-2) yr(-1) for the Luvisol, and -6.1 g m(-2) yr(-1) for the Calcaric Regosol to 3.2 g m(-2) yr(-1) for the Colluvic Regosol. The temporal variations in DOC and DIC fluxes were controlled by water fluxes. The spatially distributed leaching results corroborate the hypothesis that the effects of soil erosion influence fluxes through modified hydraulic and transport properties and terrain-dependent boundary conditions.
After more than a decade of multidisciplinary studies of the Central American subduction zone mainly in the framework of two large research programmes, the US MARGINS program and the German Collaborative Research Center SFB 574, we here review and interpret the data pertinent to quantify the cycling of mineral-bound volatiles (H2O, CO2, Cl, S) through this subduction system. For input-flux calculations, we divide the Middle America Trench into four segments differing in convergence rate and slab lithological profiles, use the latest evidence for mantle serpentinization of the Cocos slab approaching the trench, and for the first time explicitly include subduction erosion of forearc basement. Resulting input fluxes are 40-62 (53) Tg/Ma/m H2O, 7.8-11.4 (9.3) Tg/Ma/m CO2, 1.3-1.9 (1.6) Tg/Ma/m Cl, and 1.3-2.1 (1.6) Tg/Ma/m S (bracketed are mean values for entire trench length). Output by cold seeps on the forearc amounts to 0.625-1.25 Tg/Ma/m H2O partly derived from the slab sediments as determined by geochemical analyses of fluids and carbonates. The major volatile output occurs at the Central American volcanic arc that is divided into ten arc segments by dextral strike-slip tectonics. Based on volcanic edifice and widespread tephra volumes as well as calculated parental magma masses needed to form observed evolved compositions, we determine long-term (10(5) years) average magma and K2O fluxes for each of the ten segments as 32-242 (106) Tg/Ma/m magma and 0.28-2.91 (1.38) Tg/Ma/m K2O (bracketed are mean values for entire Central American volcanic arc length). Volatile/K2O concentration ratios derived from melt inclusion analyses and petrologic modelling then allow to calculate volatile fluxes as 1.02-14.3 (6.2) Tg/Ma/m H2O, 0.02-0.45 (0.17) Tg/Ma/m CO2, and 0.07-0.34 (0.22) Tg/Ma/m Cl. The same approach yields long-term sulfur fluxes of 0.12-1.08 (0.54) Tg/Ma/m while present-day open-vent SO2-flux monitoring yields 0.06-2.37 (0.83) Tg/Ma/m S. Input-output comparisons show that the arc water fluxes only account for up to 40 % of the input even if we include an "invisible" plutonic component constrained by crustal growth. With 20-30 % of the H2O input transferred into the deeper mantle as suggested by petrologic modeling, there remains a deficiency of, say, 30-40 % in the water budget. At least some of this water is transferred into two upper-plate regions of low seismic velocity and electrical resistivity whose sizes vary along arc: one region widely envelopes the melt ascent paths from slab top to arc and the other extends obliquely from the slab below the forearc to below the arc. Whether these reservoirs are transient or steady remains unknown.
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.
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.
We assessed tropical montane cloud forest (TMCF) sensitivity to natural disturbance by drought, fire, and dieback with a 7300-year-long paleorecord. We analyzed pollen assemblages, charcoal accumulation rates, and higher plant biomarker compounds (average chain length [ACL] of n-alkanes) in sediments from Wai 'anapanapa, a small lake near the upper forest limit and the mean trade wind inversion ('IWI) in Hawai`i. The paleorecord of ACL suggests increased drought frequency and a lower awl elevation from 2555-1323 cal yr B.P. and 606-334 cal yr B.P. Charcoal began to accumulate and a novel fire regime was initiated ca. 880 cal yr B.P., followed by a decreased fire return interval at ca. 550 cal yr B.P. Diebacks occurred at 2931, 2161, 1162, and 306 cal yr B.P., and two of these were independent of drought or fire. Pollen assemblages indicate that on average species composition changed only 2.8% per decade. These dynamics, though slight, were significantly associated with disturbance. The direction of species composition change varied with disturbance type. Drought was associated with significantly more vines and lianas; fire was associated with an increase in the tree fern Sadleria and indicators of open, disturbed landscapes at the expense of epiphytic ferns; whereas stand-scale dieback was associated with an increase in the tree fern Cibotium. Though this cloud forest was dynamic in response to past disturbance, it has recovered, suggesting a resilient TMCF with no evidence of state change in vegetation type (e.g., grassland or shrubland).
Sedimentary proxies used to reconstruct marine productivity suffer from variable preservation and are sensitive to factors other than productivity. Therefore, proxy calibration is warranted. Here we map the spatial patterns of two paleoproductivity proxies, biogenic opal and barium fluxes, from a set of core-top sediments recovered in the Subarctic North Pacific. Comparisons of the proxy data with independent estimates of primary and export production, surface water macronutrient concentrations, and biological pCO(2) drawdown indicate that neither proxy shows a significant correlation with primary or export productivity for the entire region. Biogenic opal fluxes, when corrected for preservation using Th-230-normalized accumulation rates, show a good correlation with primary productivity along the volcanic arcs (tau = 0.71, p = 0.0024) and with export productivity throughout the western Subarctic North Pacific (tau = 0.71, p = 0.0107). Moderate and good correlations of biogenic barium flux with export production (tau = 0.57, p = 0.0022) and with surface water silicate concentrations (tau = 0.70, p = 0.0002) are observed for the central and eastern Subarctic North Pacific. For reasons unknown, however, no correlation is found in the western Subarctic North Pacific between biogenic barium flux and the reference data. Nonetheless, we show that barite saturation, uncertainty in the lithogenic barium corrections, and problems with the reference data sets are not responsible for the lack of a significant correlation between biogenic barium flux and the reference data. Further studies evaluating the factors controlling the variability of the biogenic constituents in the sediments are desirable in this region.
Due to increasing demands and competition for high quality groundwater resources in many parts of the world, there is an urgent need for efficient methods that shed light on the interplay between complex natural settings and anthropogenic impacts. Thus a new approach is introduced, that aims to identify and quantify the predominant processes or factors of influence that drive groundwater and lake water dynamics on a catchment scale. The approach involves a non-linear dimension reduction method called Isometric feature mapping (Isomap). This method is applied to time series of groundwater head and lake water level data from a complex geological setting in Northeastern Germany. Two factors explaining more than 95% of the observed spatial variations are identified: (1) the anthropogenic impact of a waterworks in the study area and (2) natural groundwater recharge with different degrees of dampening at the respective sites of observation. The approach enables a presumption-free assessment to be made of the existing geological conception in the catchment, leading to an extension of the conception. Previously unknown hydraulic connections between two aquifers are identified, and connections revealed between surface water bodies and groundwater. (C) 2014 Elsevier B.V. All rights reserved.
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.
The tropical warm pool waters surrounding Indonesia are one of the equatorial heat and moisture sources that are considered as a driving force of the global climate system. The climate in Indonesia is dominated by the equatorial monsoon system, and has been linked to El Niño-Southern Oscillation (ENSO) events, which often result in severe droughts or floods over Indonesia with profound societal and economic impacts on the populations living in the world's fourth most populated country. The latest IPCC report states that ENSO will remain the dominant mode in the tropical Pacific with global effects in the 21st century and ENSO-related precipitation extremes will intensify. However, no common agreement exists among climate simulation models for projected change in ENSO and the Australian-Indonesian Monsoon. Exploring high-resolution palaeoclimate archives, like tree rings or varved lake sediments, provide insights into the natural climate variability of the past, and thus helps improving and validating simulations of future climate changes. Centennial tree-ring stable isotope records | Within this doctoral thesis the main goal was to explore the potential of tropical tree rings to record climate signals and to use them as palaeoclimate proxies. In detail, stable carbon (δ13C) and oxygen (δ18O) isotopes were extracted from teak trees in order to establish the first well-replicated centennial (AD 1900-2007) stable isotope records for Java, Indonesia. Furthermore, different climatic variables were tested whether they show significant correlation with tree-ring proxies (ring-width, δ13C, δ18O). Moreover, highly resolved intra-annual oxygen isotope data were established to assess the transfer of the seasonal precipitation signal into the tree rings. Finally, the established oxygen isotope record was used to reveal possible correlations with ENSO events. Methodological achievements | A second goal of this thesis was to assess the applicability of novel techniques which facilitate and optimize high-resolution and high-throughput stable isotope analysis of tree rings. Two different UV-laser-based microscopic dissection systems were evaluated as a novel sampling tool for high-resolution stable isotope analysis. Furthermore, an improved procedure of tree-ring dissection from thin cellulose laths for stable isotope analysis was designed. The most important findings of this thesis are: I) The herein presented novel sampling techniques improve stable isotope analyses for tree-ring studies in terms of precision, efficiency and quality. The UV-laser-based microdissection serve as a valuable tool for sampling plant tissue at ultrahigh-resolution and for unprecedented precision. II) A guideline for a modified method of cellulose extraction from wholewood cross-sections and subsequent tree-ring dissection was established. The novel technique optimizes the stable isotope analysis process in two ways: faster and high-throughput cellulose extraction and precise tree-ring separation at annual to high-resolution scale. III) The centennial tree-ring stable isotope records reveal significant correlation with regional precipitation. High-resolution stable oxygen values, furthermore, allow distinguishing between dry and rainy season rainfall. IV) The δ18O record reveals significant correlation with different ENSO flavors and demonstrates the importance of considering ENSO flavors when interpreting palaeoclimatic data in the tropics. The findings of my dissertation show that seasonally resolved δ18O records from Indonesian teak trees are a valuable proxy for multi-centennial reconstructions of regional precipitation variability (monsoon signals) and large-scale ocean-atmosphere phenomena (ENSO) for the Indo-Pacific region. Furthermore, the novel methodological achievements offer many unexplored avenues for multidisciplinary research in high-resolution palaeoclimatology.
Recultivation of disturbed oil sand mining areas is an issue of increasing importance. Nevertheless only little is known about the fate of organic matter, cell abundances and microbial community structures during oil sand processing, tailings management and initial soil development on reclamation sites. Thus the focus of this work is on biogeochemical changes of mined oil sands through the entire process chain until its use as substratum for newly developing soils on reclamation sites. Therefore, oil sand, mature fine tailings (MFTs) from tailings ponds and drying cells and tailings sand covered with peat-mineral mix (PMM) as part of land reclamation were analyzed. The sample set was selected to address the question whether changes in the above-mentioned biogeochemical parameters can be related to oil sand processing or biological processes and how these changes influence microbial activities and soil development.
GC-MS analyses of oil-derived biomarkers reveal that these compounds remain unaffected by oil sand processing and biological activity. In contrast, changes in polycyclic aromatic hydrocarbon (PAH) abundance and pattern can be observed along the process chain. Especially naphthalenes, phenanthrenes and chrysenes are altered or absent on reclamation sites, Furthermore, root-bearing horizons on reclamation sites exhibit cell abundances at least ten times higher (10(8) to 10(9) cells g(-1)) than in oil sand and MFF samples (10(7) cells g(-1)) and show a higher diversity in their microbial community structure. Nitrate in the pore water and roots derived from the PMM seem to be the most important stimulants for microbial growth. The combined data show that the observed compositional changes are mostly related to biological activity and the addition of exogenous organic components (PMM), whereas oil extraction, tailings dewatering and compaction do not have significant influences on the evaluated compounds. Microbial community composition remains relatively stable through the entire process chain. (C) 2014 Elsevier B.V. All rights reserved.
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 International Union of Geological Sciences (JUGS) is evaluating whether there are additional geoscientific activities that would be beneficial in helping mitigate the impacts of tsunami. Public concerns about poor decisions and inaction, and advances in computing power and data mining call for new scientific approaches. Three fundamental requirements for mitigating impacts of natural hazards are defined. These are: (1) improvement of process-oriented understanding, (2) adequate monitoring and optimal use of data, and (3) generation of advice based on scientific, technical and socio-economic expertise. International leadership/coordination is also important.
To increase the capacity to predict and mitigate the impacts of tsunami and other natural hazards a broad consensus is needed. The main needs include the integration of systematic geological inputs - identifying and studying paleo-tsunami deposits for all subduction zones; optimising coverage and coordination of geodetic and seismic monitoring networks; underpinning decision making at national and international scales by developing appropriate mechanisms for gathering, managing and communicating authoritative scientific and technical advice information; international leadership for coordination and authoritative statements of best approaches. All these suggestions are reflected in the Sendai Agreement, the collective views of the experts at the International Workshop on Natural Hazards, presented later in this volume.
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.
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.
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.
Degrading permafrost can alter ecosystems, damage infrastructure, and release enough carbon dioxide (CO2) and methane (CH4) to influence global climate. The permafrost carbon feedback (PCF) is the amplification of surface warming due to CO2 and CH4 emissions from thawing permafrost. An analysis of available estimates PCF strength and timing indicate 120 +/- 85 Gt of carbon emissions from thawing permafrost by 2100. This is equivalent to 5.7 +/- 4.0% of total anthropogenic emissions for the Intergovernmental Panel on Climate Change (IPCC) representative concentration pathway (RCP) 8.5 scenario and would increase global temperatures by 0.29 +/- 0.21 degrees C or 7.8 +/- 5.7%. For RCP4.5, the scenario closest to the 2 degrees C warming target for the climate change treaty, the range of cumulative emissions in 2100 from thawing permafrost decreases to between 27 and 100 Gt C with temperature increases between 0.05 and 0.15 degrees C, but the relative fraction of permafrost to total emissions increases to between 3% and 11%. Any substantial warming results in a committed, long-term carbon release from thawing permafrost with 60% of emissions occurring after 2100, indicating that not accounting for permafrost emissions risks overshooting the 2 degrees C warming target. Climate projections in the IPCC Fifth Assessment Report (AR5), and any emissions targets based on those projections, do not adequately account for emissions from thawing permafrost and the effects of the PCF on global climate. We recommend the IPCC commission a special assessment focusing on the PCF and its impact on global climate to supplement the AR5 in support of treaty negotiation.
The biostratigraphy of Campanian-Maastrichtian carbonate platforms is largely based on the larger foraminiferal genus Orbitoides. However, while the taxonomy and the chronostratigraphic age of the younger species of this genus are well established, there are still many controversies on the earliest species. We have restudied their morphological characters using a large collection of samples from the type-localities and from continuous sections in the southern Pyrenees. Based on these new observations, the long forgotten species O. sanctae-pelagiae is reinstated, while O. dordoniensis is considered a junior synonym. Successive populations of O. hottingeri, O. sanctae-pelagiae and O. douvillei show gradual morphological changes in time marked by an increase in the size and complexity of the macrospheric embryonal apparatus, an increase of the size of the adult specimens of both generations and the progressive appearance and development of true lateral chamberlets. The Font de les Bagasses Unit in the southern Pyrenees preserves a high-resolution archive of the evolution of the earliest Orbitoides. Strontium isotope stratigraphy indicates that the oldest species, O. hottingeri, made its first appearance in the earliest Campanian, close to the Santonian-Campanian boundary, and was replaced by O. sanctae-pelagiae at a level closely corresponding to the boundary between the Placenticeras bidorsatum and Menabites delawarensis ammonite zones. (C) 2014 Elsevier Ltd. All rights reserved.
The water-level record from the 300 m deep paleo-lake Suguta (Northern Kenya Rift) during the African Humid Period (AHP, 15-5 ka BP) helps to explain decadal to centennial intensity variations in the West African Monsoon (WAM) and the Indian Summer Monsoon (ISM). This water-level record was derived from three different sources: (1) grain size variations in radiocarbon dated and reservoir corrected lacustrine sediments, (2) the altitudes and ages of paleo-shorelines within the basin, and (3) the results of hydro-balance modeling, providing important insights into the character of water level variations (abrupt or gradual) in the amplifier paleo-Lake Suguta. The results of these comprehensive analyses suggest that the AHP highstand in the Suguta Valley was the direct consequence of a northeastwards shift in the Congo Air Boundary (CAB), which was in turn caused by an enhanced atmospheric pressure gradient between East Africa and India during a northern hemisphere insolation maximum. Rapidly decreasing water levels of up to 90 m over less than a hundred years are best explained by changes in solar irradiation either reducing the East African-Indian atmospheric pressure gradient and preventing the CAB from reaching the study area, or reducing the overall humidity in the atmosphere, or a combination of both these effects. In contrast, although not well documented in our record we hypothesize a gradual end of the AHP despite an abrupt change in the source of precipitation when a decreasing pressure gradient between Asia and Africa prevented the CAB from reaching the Suguta Valley. The abruptness was probably buffered by a contemporaneous change in precession producing an insolation maximum at the equator during October. Whether or not this is the case, the water-level record from the Suguta Valley demonstrates the importance of both orbitally-controlled insolation variations and short-term changes in solar irradiation as factors affecting the significant water level variations in East African rift lakes.
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.
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.
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).
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.
Deep seepage estimation is important for water balance investigations of groundwater and the vadose zone. A simplified Buckingham-Darcy method to assess time series of deep seepage fluxes was proposed by Schindler and Muller (1998). In the method dynamics of water fluxes are calculated by a soil hydraulic conductivity function. Measured soil moistures and matric heads are used as input data. Resulting time series of flux dynamics are scaled to realistic absolute levels by calibrating the method with the areal water balance. An assumption of the method is that water fluxes at different positions exhibit identical dynamics although their absolute values can differ. The aim of this study was to investigate uncertainties of that method depending on the particle size distribution and textural heterogeneity in non-layered soils. We performed a numerical experiment using the two-dimensional Richards Equation. A basic model of transient water fluxes beneath the root and capillary zone was setup and used to simulate time series of soil moisture, matric head, and seepage fluxes for 4221 different cases of particle size distribution and intensities of textural heterogeneity. Soil hydraulic parameters were predicted by the pedotransfer function Rosetta. Textural heterogeneity was modeled with Miller and Miller scaling factors arranged in spatial random fields. Seepage fluxes were calculated with the Buckingham-Darcy method from simulated soil moisture and matric head time series and compared with simulated reference fluxes. The median of Root Mean Square Error was about 0.026 cm d(-1) and the median of maximum cross correlation was 0.96 when the method was calibrated adequately. The method's performance was mainly influenced by (i) the soil textural class and (ii) the time period used for flux calibration. It performed best in sandy loam while hotspots of errors occurred in sand and silty texture. Calibrating the method with time periods that exhibit high variance of seepage fluxes yielded the best performance. The geostatistical properties of the Miller and Miller scaling field influenced the performance only slightly. However, the Miller and Miller scaling procedure generated heterogeneous flow fields that were addressed as main reason for mismatches of simulated reference fluxes and fluxes obtained with the Buckingham-Darcy method.
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.
The Himalayan mountains are dissected by some of the deepest and most impressive gorges on Earth. Constraining the interplay between river incision and rock uplift is important for understanding tectonic deformation in this region. We report here the discovery of a deeply incised canyon of the Yarlung Tsangpo River, at the eastern end of the Himalaya, which is now buried under more than 500 meters of sediments. By reconstructing the former valley bottom and dating sediments at the base of the valley fill, we show that steepening of the Tsangpo Gorge started at about 2 million to 2.5 million years ago as a consequence of an increase in rock uplift rates. The high erosion rates within the gorge are therefore a direct consequence of rapid rock uplift.
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.
In hydrology, the storage-discharge relationship is a fundamental catchment property. Understanding what controls this relationship is at the core of catchment science. To date, there are no direct methods to measure water storage at catchment scales (10(1)-10(3)km(2)). In this study, we use direct measurements of terrestrial water storage dynamics by means of superconducting gravimetry in a small headwater catchment of the Regen River, Germany, to derive empirical storage-discharge relationships in nested catchments of increasing scale. Our results show that the local storage measurements are strongly related to streamflow dynamics at larger scales (> 100km(2); correlation coefficient=0.78-0.81), but at small scale, no such relationship exists (similar to 1km(2); correlation coefficients=-0.11). The geologic setting in the region can explain both the disconnection between local water storage and headwater runoff, and the connectivity between headwater storage and streams draining larger catchment areas. More research is required to understand what controls the form of the observed storage-discharge relationships at the catchment scale. This study demonstrates that high-precision gravimetry can provide new insights into the complex relationship between state and response of hydrological systems.
We have investigated the influence that megathrust earthquake slip has on the activation of splay faults using a 2-D finite element method (FEM), taking into account the effects of gravity and variations in the frictional strength properties of splay faults. We simulated both landward-dipping and seaward-dipping splay fault geometries, and imposed depth-variable slip distributions of subduction events. Our results indicate that the two types of splay fault exhibit a similar behavior, with variations in frictional properties along the faults affecting only the seismic magnitude. The triggering process is controlled by a critical depth. Megathrust slip concentrated at depths shallower than the critical depth will favor normal displacement, while megathrust slip concentrated at depths deeper than the critical depth is likely to result in reverse motion. Our results thus provide a useful tool for predicting the activation of secondary faults and may have direct implications for tsunami hazard research.
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.
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 spatiotemporal, kinematic, and source characteristics of induced seismicity occurring at different fluid injection rates are investigated to determine the predominant physical mechanisms responsible for induced seismicity at the northwestern part of The Geysers geothermal field, California. We analyze a relocated hypocenter catalog from a seismicity cluster where significant variations of the stress tensor orientation were previously observed to correlate with injection rates. We find that these stress tensor orientation changes may be related to increased pore pressure and the corresponding changes in poroelastic stresses at reservoir depth. Seismic events during peak injections tend to occur at greater distances from the injection well, preferentially trending parallel to the maximum horizontal stress direction. In contrast, at lower injection rates the seismicity tends to align in a different direction which suggests the presence of a local fault. During peak injection intervals, the relative contribution of strike-slip faulting mechanisms increases. Furthermore, increases in fluid injection rates also coincide with a decrease in b values. Our observations suggest that regardless of the injection stage, most of the induced seismicity results from thermal fracturing of the reservoir rock. However, during peak injection intervals, the increase in pore pressure may likewise be responsible for the induced seismicity. By estimating the thermal and hydraulic diffusivities of the reservoir, we confirm that the characteristic diffusion length for pore pressure is much greater than the corresponding length scale for temperature and also more consistent with the spatial extent of seismicity observed during different injection rates.
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.
Remote sensing analysis is a crucial tool for monitoring the extent of mine waste surfaces and their mineralogy in countries with a long mining history, such as South Africa, where gold and platinum have been produced for over 90 years. These mine waste sites have the potential to contain problematic trace element species (e. g., U, Pb, Cr). In our research, we aim to combine the mapping and monitoring capacities of multispectral and hyperspectral spaceborne sensors. This is done to assess the potential of existing multispectral and hyperspectral spaceborne sensors (OLI and Hyperion) and future missions, such as Sentinel-2 and EnMAP (Environmental Mapping and Analysis Program), for mapping the spatial extent of these mine waste surfaces. For this task we propose a new index, termed the iron feature depth (IFD), derived from Landsat-8 OLI data to map the 900-nm absorption feature as a potential proxy for monitoring the spatial extent of mine waste. OLI was chosen, because it represents the most suitable sensor to map the IFD over large areas in a multi-temporal manner due to its spectral band layout; its (183 km x 170 km) scene size and its revisiting time of 16 days. The IFD is in good agreement with primary and secondary iron-bearing minerals mapped by the Material Identification and Characterization Algorithm (MICA) from EO-1 Hyperion data and illustrates that a combination of hyperspectral data (EnMAP) for mineral identification with multispectral data (Sentinel-2) for repetitive area-wide mapping and monitoring of the IFD as mine waste proxy is a promising application for future spaceborne sensors. A maximum, absolute model error is used to assess the ability of existing and future multispectral sensors to characterize mine waste via its 900-nm iron absorption feature. The following sensor-signal similarity ranking can be established for spectra from gold mining material: EnMAP 100% similarity to the reference, ALI 97.5%, Sentinel-2 97%, OLI and ASTER 95% and ETM+ 91% similarity.
Many mountain belts sustain prolonged snow cover for parts of the year, although enquiries into rates of erosion in these landscapes have focused almost exclusively on the snow-free periods. This raises the question of whether annual snow cover contributes significantly to modulating rates of erosion in high-relief terrain. In this context, the sudden release of snow avalanches is a frequent and potentially relevant process, judging from the physical damage to subalpine forest ecosystems, and the amount of debris contained in avalanche deposits. To quantitatively constrain this visual impression and to expand the sparse literature, we sampled sediment concentrations of n = 28 river-spanning snow-avalanche deposits (snow bridges) in the area around Davos, eastern Swiss Alps, and inferred an orders-of-magnitude variability in specific fine sediment and organic carbon yields (1.8 to 830 t km(-2) yr(-1), and 0.04 to 131 tC km(-2) yr(-1), respectively). A Monte Carlo simulation demonstrates that, with a minimum of free parameters, such variability is inherent to the geometric scaling used for computing specific yields. Moreover, the widely applied method of linearly extrapolating plot scale sample data may be prone to substantial under- or overestimates. A comparison of our inferred yields with previously published work demonstrates the relevance of wet snow avalanches as prominent agents of soil erosion and transporters of biogeochemical constituents to mountain rivers. Given that a number of snow bridges persisted below the insulating debris cover well into the summer months, snow-avalanche deposits also contribute to regulating in-channel sediment and organic debris storage on seasonal timescales. Finally, our results underline the potential shortcomings of neglecting erosional processes in the winter and spring months in mountainous terrain subjected to prominent snow cover.
This paper presents a Bayesian non-parametric method based on Gaussian Process (GP) regression to derive ground-motion models for peak-ground parameters and response spectral ordinates. Due to its non-parametric nature there is no need to specify any fixed functional form as in parametric regression models. A GP defines a distribution over functions, which implicitly expresses the uncertainty over the underlying data generating process. An advantage of GP regression is that it is possible to capture the whole uncertainty involved in ground-motion modeling, both in terms of aleatory variability as well as epistemic uncertainty associated with the underlying functional form and data coverage. The distribution over functions is updated in a Bayesian way by computing the posterior distribution of the GP after observing ground-motion data, which in turn can be used to make predictions. The proposed GP regression models is evaluated on a subset of the RESORCE data base for the SIGMA project. The experiments show that GP models have a better generalization error than a simple parametric regression model. A visual assessment of different scenarios demonstrates that the inferred GP models are physically plausible.
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.
TopoToolbox is a MATLAB program for the analysis of digital elevation models (DEMs). With the release of version 2, the software adopts an object-oriented programming (OOP) approach to work with gridded DEMs and derived data such as flow directions and stream networks. The introduction of a novel technique to store flow directions as topologically ordered vectors of indices enables calculation of flow-related attributes such as flow accumulation similar to 20 times faster than conventional algorithms while at the same time reducing memory overhead to 33% of that required by the previous version. Graphical user interfaces (GUIs) enable visual exploration and interaction with DEMs and derivatives and provide access to tools targeted at fluvial and tectonic geomorphologists. With its new release, TopoToolbox has become a more memory-efficient and faster tool for basic and advanced digital terrain analysis that can be used as a framework for building hydrological and geomorphological models in MATLAB.
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.