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)
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.
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.
Trees and shrubs in tropical Africa use the C-3 cycle as a carbon fixation pathway during photosynthesis, while grasses and sedges mostly use the C-4 cycle. Leaf-wax lipids from sedimentary archives such as the long-chain n-alkanes (e.g., n-C-27 to n-C-33) inherit carbon isotope ratios that are representative of the carbon fixation pathway. Therefore, n-alkane delta C-13 values are often used to reconstruct past C-3/C-4 composition of vegetation, assuming that the relative proportions of C-3 and C-4 leaf waxes reflect the relative proportions of C-3 and C-4 plants. We have compared the delta C-13 values of n-alkanes from modern C-3 and C-4 plants with previously published values from recent lake sediments and provide a framework for estimating the fractional contribution (areal-based) of C-3 vegetation cover (f(C3)) represented by these sedimentary archives. Samples were collected in Cameroon, across a latitudinal transect that accommodates a wide range of climate zones and vegetation types, as reflected in the progressive northward replacement of C-3-dominated rain forest by C-4-dominated savanna. The C-3 plants analysed were characterised by substantially higher abundances of n-C-29 alkanes and by substantially lower abundances of n-C-33 alkanes than the C-4 plants. Furthermore, the sedimentary delta C-13 values of n-C-29 and n-C-31 alkanes from recent lake sediments in Cameroon (-37.4%) to 26.5%) were generally within the range of delta C-13 values for C-3 plants, even when from sites where C-4 plants dominated the catchment vegetation. In such cases simple linear mixing models fail to accurately reconstruct the relative proportions of C-3 and C-4 vegetation cover when using the delta C-13 values of sedimentary n-alkanes, overestimating the proportion of C-3 vegetation, likely as a consequence of the differences in plant wax production, preservation, transport, and/or deposition between C-3 and C-4 plants. We therefore tested a set of non-linear binary mixing models using delta C-13 values from both C-3 and C-4 vegetation as end-members. The non-linear models included a sigmoid function (sine-squared) that describes small variations in the f(C3) values as the minimum and maximum delta C-13 values are approached, and a hyperbolic function that takes into account the differences between C-3 and C-4 plants discussed above. Model fitting and the estimation of uncertainties were completed using the Monte Carlo algorithm and can be improved by future data addition. Models that provided the best fit with the observed delta C-13 values of sedimentary n-alkanes were either hyperbolic functions or a combination of hyperbolic and sine-squared functions. Such non-linear models may be used to convert delta C-13 measurements on sedimentary n-alkanes directly into reconstructions of C-3 vegetation cover. (C) 2014 Elsevier Ltd. All rights reserved.
Kinematic interaction of faults is an important issue for detailed seismic hazard assessments in seismically active regions. The Eastern Mosha Fault (EMF) and the North Tehran Fault (NTF) are two major active faults of the southern central Alborz mountains, located in proximity of Tehran (population similar to 9 million). We used field, geomorphological and paleoseismological data to explore the kinematic transition between the faults, and compare their short-term and long-term history of activity. We introduce the Niknamdeh segment of the NTF along which the strike-slip kinematics of EMF is transferred onto the NTF, and which is also responsible for the development of a pull-apart basin between the eastern segments of the NTF. The Ira trench site at the linkage zone between the two faults reveals the history of interaction between rock avalanches, active faulting and sag-pond development. The kinematic continuity between the EMF and NTF requires updating of seismic hazard models for the NTF, the most active fault adjacent to the Tehran Metropolis. Study of offsets of large-scale morphological features along the EMF, and comparison with estimated slip rates along the fault indicates that the EMF has started its left-lateral kinematics between 3.2 and 4.7 Ma. According to our paleoseismological data and the morphology of the nearby EMF and NTF, we suggest minimum and maximum values of about 1.8 and 3.0 mm/year for the left-lateral kinematics on the two faults in their linkage zone, averaged over Holocene time scales. Our study provides a partial interpretation, based on available data, for the fault activity in northeastern Tehran region, which may be completed with studies of other active faults of the region to evaluate a more realistic seismic hazard analysis for this heavily populated major city. (C) 2014 Elsevier B.V. All rights reserved.
Members of the acI lineage of Actinobacteria are the most abundant microorganisms in most freshwater lakes; however, our understanding of the keys to their success and their role in carbon and nutrient cycling in freshwater systems has been hampered by the lack of pure cultures and genomes. We obtained draft genome assemblies from 11 single cells representing three acI tribes (acI-A1, acI-A7, acI-B1) from four temperate lakes in the United States and Europe. Comparative analysis of acI SAGs and other available freshwater bacterial genomes showed that acI has more gene content directed toward carbohydrate acquisition as compared to Polynucleobacter and LD12 Alphaproteobacteria, which seem to specialize more on carboxylic acids. The acI genomes contain actinorhodopsin as well as some genes involved in anaplerotic carbon fixation indicating the capacity to supplement their known heterotrophic lifestyle. Genome-level differences between the acI-A and acI-B clades suggest specialization at the clade level for carbon substrate acquisition. Overall, the acI genomes appear to be highly streamlined versions of Actinobacteria that include some genes allowing it to take advantage of sunlight and N-rich organic compounds such as polyamines, di-and oligopeptides, branched-chain amino acids and cyanophycin. This work significantly expands the known metabolic potential of the cosmopolitan freshwater acI lineage and its ecological and genetic traits.
Inferring a ground-motion prediction equation (GMPE) for a region in which only a small number of seismic events has been observed is a challenging task. A response to this data scarcity is to utilise data from other regions in the hope that there exist common patterns in the generation of ground motion that can contribute to the development of a GMPE for the region in question. This is not an unreasonable course of action since we expect regional GMPEs to be related to each other. In this work we model this relatedness by assuming that the regional GMPEs occupy a common low-dimensional manifold in the space of all possible GMPEs. As a consequence, the GMPEs are fitted in a joint manner and not independent of each other, borrowing predictive strength from each other's regional datasets. Experimentation on a real dataset shows that the manifold assumption displays better predictive performance over fitting regional GMPEs independent of each other. (C) 2014 Elsevier Ltd. All rights reserved.
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.
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.
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.
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.
In situ UV laser spot Ar-40/Ar-39 analyses of distinct phengite types in eclogite-facies rocks from the Sesia-Lanzo Zone (Western Alps, Italy) were combined with SIMS boron isotope analyses as well as boron (B) and lithium (Li) concentration data to link geochronological information with constraints on fluid-rock interaction. In weakly deformed samples, apparent Ar-40/Ar-39 ages of phengite cores span a range of similar to 20 Ma, but inverse isochrons define two distinct main high-pressure (HP) phengite core crystallization periods of 88-82 and 77-74 Ma, respectively. The younger cores have on average lower B contents (similar to 36 mu g/g) than the older ones (similar to 43-48 mu g/g), suggesting that loss of B and resetting of the Ar isotopic system were related. Phengite cores have variable delta B-11 values (-18 parts per thousand to -10 parts per thousand), indicating the lack of km scale B homogenization during HP crystallization.
Overprinted phengite rims in the weakly deformed samples generally yield younger apparent Ar-40/Ar-39 ages than the respective cores. They also show variable effects of heterogeneous excess 40 Ar incorporation and Ar loss. One acceptable inverse isochron age of 77.1 +/- 1.1 Ma for rims surrounding older cores (82.6 +/- 0.6 Ma) overlaps with the second period of core crystallization. Compared to the phengite cores, all rims have lower B and Li abundances but similar delta B-11 values (-15 parts per thousand to -9 parts per thousand), reflecting internal redistribution of B and Li and internal fluid buffering of the B isotopic composition during rim growth. The combined observation of younger Ar-40/Ar-39 ages and boron loss, yielding comparable values of both parameters only in cores and rims of different samples, is best explained by a selective metasomatic overprint. In low permeability samples, this overprint caused recrystallization of phengite rims, whereas higher permeability in other samples led to complete recrystallization of phengite grains.
Strongly deformed samples from a several km long, blueschist-facies shear zone contain mylonitic phengite that forms a tightly clustered group of relatively young apparent Ar-40/Ar-39 ages (64.7-68.8 Ma), yielding an inverse isochron age of 65.0 +/- 3.0 Ma. Almost complete B and Li removal in mylonitic phengite is due to leaching into a fluid. The B isotopic composition is significantly heavier than in phengites from the weakly deformed samples, indicating an external control by a high-delta B-11 fluid (delta B-11 = + 7 +/- 4 parts per thousand). We interpret this result as reflecting phengite recrystallization related to deformation and associated fluid flow in the shear zone. This event also caused partial resetting of the Ar isotope system and further B loss in more permeable rocks of the adjacent unit. We conclude that geochemical evidence for pervasive or limited fluid flow is crucial for the interpretation of Ar-40/Ar-39 data in partially metasomatized rocks.
Velocity models are essential to process two-and three-dimensional ground-penetrating radar (GPR) data. Furthermore, velocity information aids the interpretation of such data sets because velocity variations reflect important material properties such as water content. In many GPR applications, common midpoint (CMP) surveys are routinely collected to determine one-dimensional velocity models at selected locations. To analyse CMP data gathers, spectral velocity analyses relying on the normal-moveout (NMO) model are commonly employed. Using Dix's formula, the derived NMO velocities can be further converted to interval velocities which are needed for processing and interpretation. Because of the inherent assumptions and limitations of such approaches, we investigate and propose an alternative procedure based on the global inversion of reflection travel-times. We use a finite-difference solver of the Eikonal equation to accurately solve the forward problem in combination with particle swarm optimization (PSO) to find one-dimensional GPR velocity models explaining our data. Because PSO is a robust and efficient global optimization tool, our inversion approach includes generating an ensemble of representative solutions that allows us to analyse uncertainties in the model space. Using synthetic data examples, we test and evaluate our inversion approach to analyse CMP data collected across typical near-surface environments. Application to a field data set recorded at a well-constrained test site including a comparison to independent borehole and direct-push data, further illustrates the potential of the proposed approach, which includes a straightforward and understandable appraisal of non-uniqueness and uncertainty issues, respectively. We conclude that our methodology is a feasible and powerful tool to analyse GPR CMP data and allows practitioners and researchers to evaluate the reliability of CMP derived velocity models.
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.
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.
The ratio of unsupported protactinium-231 to thorium-230 in marine sediments, (Pa/Th)(xs), is potentially sensitive to several processes of oceanographic and climatological interest: deep ocean circulation, marine biological productivity (as it relates to total particle flux) and particle composition (specifically, biogenic opal and authigenic Mn). In order to attribute variations in (Pa/Th)(xs) observed in sediment records to changes in specific processes through time, a better understanding of the chemical cycling of these elements in the modern ocean is necessary. To this end, a survey was undertaken of (Pa/Th)(xs) in surface sediments from the subarctic Pacific (SO202-INOPEX expedition) in combination with a Pacific-wide compilation of published data. Throughout the Pacific, (Pa/Th)(xs) is robustly correlated with the opal content of sediments. In the North and equatorial Pacific, simultaneous positive correlations with productivity indicators suggest that boundary scavenging and opal scavenging combine to enhance the removal of Pa in the eastern equatorial Pacific and subarctic Pacific. Deep ocean water mass ageing (>3.5 km) associated with the Pacific overturning appears to play a secondary role in determining the basin scale distribution of (Pa/Th)(xs). A basin-wide extrapolation of Pa removal is performed which suggests that the Pacific Pa budget is nearly in balance. We hypothesize that through time (Pa/Th)(xs) distributions in the Pacific could define the evolving boundaries of contrasting biogeographic provinces in the North Pacific, while the influence of hydrothermal scavenging of Pa potentially confounds this approach in the South Pacific.
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.
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.
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.
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.
(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.
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.
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.
The hydrology of western Central Asia is highly sensitive to climatic perturbations. In order to understand its long-term variability and to infer linkages between precipitation and atmospheric and oceanic systems, we conducted a thorough sedimentary and geochemical study on a composite core retrieved in lake Son Kul (central Kyrgyzstan). A multi-proxy approach was conducted on lake sediments based on grain size analyses, magnetic susceptibility, total organic carbon (TOC), total nitrogen (TN) and carbon and oxygen isotope analyses on bulk and biogenic materials (ostracoda and molluscs shells) at a resolution equivalent to ca 40 years, aiming to characterise the sequence of palaeolimnological changes in Son Kul.
As indicated by delta O-18 record of bulk carbonates, mainly consisting of aragonite, the Holocene hydrological balance was negative during most of time, suggesting an excess of evaporation (E) over precipitation (P). Limnological conditions fluctuated rapidly before 5000 cal yr BP indicating significant changes in regional hydrology and climate. In particular, the long-term negative hydrological balance was impeded by several short stages with marked increase of precipitation, lasting several decades to a few centuries (e.g., 8300-8200, 6900-6700, 6300-6100, 5500-5400, 5300-5200 and 3100 -3000 cal yr BP). Precipitation changes as inferred from 8180 data are also documented by increased minerogenic detritus and higher TOC. We propose that the seasonal pattern of precipitation varied transiently in western Central Asia during the Holocene, although evaporation changes may also account for the rapid changes observed in delta O-18 data. When the annual water balance was less critical (P <= E), the excess of water might be ascribed to increased precipitation during cold seasons mainly because winter precipitation has more negative delta O-18 than its summer equivalent. Conversely, when the annual water balance is negative (P E), the moisture was mainly delivered during the warm season, as between 5000 and 2000 cal yr BP. Our results thus imply that moisture sources could have changed as well during the Holocene. Moisture was delivered as today mainly during summer from the extended Caspian-Aral Basin and eastern Mediterranean, although Arctic and even North Atlantic seas might be important moisture sources when seasonal precipitation was dominated by winter precipitation. We hypothesise that warming Arctic and North Atlantic seas were important for the North Hemisphere circulation during the cold season. (C) 2014 Elsevier Ltd. All rights reserved.
Integrated biostratigraphic-chemostratigraphic studies provide evidence that the proto-North Atlantic realm witnessed major changes in carbonate platform production in the run-up of the Early Aptian oceanic anoxic event (OAE) 1a. Whereas pervasive growth of Lithocodium microencrusters represents an early harbinger of OAE1a-related environmental perturbation, the subsequent replacement of oligotrophic rudist-coral-nerineid by mesotrophic orbitolinid-oyster communities was clearly associated with the event itself. In order to test the supra-regional relevance of this major community replacement, two shallow-water sections in the southern Lusitanian Basin (Portugal) are investigated by means of geochemistry (carbon and oxygen isotopes), cement petrography and detailed sedimentological analysis. The focus is on a regional, prominent discontinuity surface (S4) at the transition between oligotrophic and mesotrophic carbonate platform production, which might indicate that the major biotic change could have been associated with a phase of non-sedimentation and possibly erosion. The studied sections (Sao Julia, Crismina) provide evidence that the major Early Aptian biotic turnover was preceded by numerous subordinate but significant changes in platform ecology, which mirrored a series of progressive short-term environmental changes in the course of OAE1. Several transient mass occurrences of orbitolinids indicate repeated phases of ecological stress arguably due to enhanced nutrient input and deepening. Small-scale sea-level changes at parasequence level below the major discontinuity surface are revealed by alternations of rudist assemblages dominated by clinger or recumbent forms as well as intercalated hardground and subaerial exposure stages. Expanded phases of subaerial exposure, however, can be largely ruled out following the geochemical and cement-petrographic data presented here. Enhanced continent-derived siliciclastic input characterising the lower orbitolinid-oyster dominated limestones is in support of a shift to more humid conditions during the middle Early Aptian. This is in line with palaeoclimatic data, which propose a southward movement of the mid-latitude arid climate belt during this time. The documented rapid replacement of oligotrophic assemblages by various environmental-stress adapted carbonate platform communities might be seen as explanation for ongoing Early Aptian proto-North Atlantic carbonate production during a time of widespread platform demise and drowning in the northern Tethyan realm. (C) 2014 Elsevier B.V. All rights reserved.
Biomass allometries and coarse root biomass distribution of mountain birch in southern Iceland
(2014)
Root systems are an important pool of biomass and carbon in forest ecosystems. However, most allometric studies on forest trees focus only on the aboveground components. When estimated, root biomass has most often been calculated by using a fixed conversion factor from aboveground biomass. In order to study the size-related development of the root system of native mountain birch (Betula pubescens Ehrh. ssp. czerepanovii), we collected the coarse root system of 25 different aged birch trees (stem diameter at 50 cm length between 0.2 and 14.1 cm) and characterized them by penetration depth (< 1 m) and root thickness. Based on this dataset, allometric functions for coarse roots (> 5 mm and > 2 mm), root stock, total belowground biomass and aboveground biomass components were calculated by a nonlinear and a linear fitting approach. The study showed that coarse root biomass of mountain birch was almost exclusively (> 95 weight-%) located in the top 30 cm, even in a natural old-growth woodland. By using a cross-validation approach, we found that the nonlinear fitting procedure performed better than the linear approach with respect to predictive power. In addition, our results underscore that general assumptions of fixed conversion factors lead to an underestimation of the belowground biomass. Thus, our results provide allometric functions for a more accurate root biomass estimation to be utilized in inventory reports and ecological studies.
The century-long debate over the origins of inner gorges that were repeatedly covered by Quaternary glaciers hinges upon whether the gorges are fluvial forms eroded by subaerial rivers, or subglacial forms cut beneath ice. Here we apply cosmogenic nuclide exposure dating to seven inner gorges along similar to 500 km of the former Fennoscandian ice sheet margin in combination with a new deglaciation map. We show that the timing of exposure matches the advent of ice-free conditions, strongly suggesting that gorges were cut by channelized subglacial meltwater while simultaneously being shielded from cosmic rays by overlying ice. Given the exceptional hydraulic efficiency required for meltwater channels to erode bedrock and evacuate debris, we deduce that inner gorges are the product of ice sheets undergoing intense surface melting. The lack of postglacial river erosion in our seven gorges implicates subglacial meltwater as a key driver of valley deepening on the Baltic Shield over multiple glacial cycles.
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.
In this paper we present the sedimentary features and the luminescence chronology for Nahal Mahanayeem Outlet (NMO), an archaeological open air site at the southern margin of the Hula Basin (Northern Jordan Rift Valley, Israel). The site is characterized by a lithic assemblage ascribed primarily to the Middle Paleolithic Mousterian tradition, and by an excellent preservation of floral and faunal remains.
Six geological units forming the stratigraphic sequence of the site were distinguished: (Unit 6) archaeologically sterile, light-colored limnic carbonates; (Unit 5) conglomerates of rounded basalt boulders and cobbles forming a hill-like topography; (Unit 4 and Unit 3) a sequence of similar dark silty sediments, attached to and overlaying the conglomerates, containing the archaeological horizons of the site; (Unit 2) a number of channels cutting into the top of Unit 3, filled with coarse sand and rounded basalt and limestone gravels of fluvial origin; and (Unit 1) a thin sand layer laid down by the present-day Jordan River covering another unconformity as a result of heavy machinery drainage operations in 1999.
The OSL age for Unit 6 yielded a minimum age older than 460 ka. Sedimentary features and the embedded fossils suggest that Unit 1 can be linked to the Early Pleistocene Gadot Chalk. Unit 5 represents a local geological feature and could be an indicator for a period of increased erosion with formation of coarse grained sediments. The archaeological horizons form the lower parts of Unit 4 and yielded OSL-ages between 55 and 65 ka, indicating an affiliation to the sediments called "Ashmura Formation" with an Upper Pleistocene age for the site. The channel fills of Unit 2 can be dated by the recovered artifacts. These range in age from the Upper Paleolithic (Aurignacian) to historic times. Unit 1 is recent.
The study of the complex NMO stratigraphy, combined with coherent OSL chronology, has enabled us to reconstruct parts of the geological history of the Hula Basin during the Late Pleistocene. It is this history that forms the background for the human migration and utilization of natural resources in the Upper Jordan Rift Valley. (C) 2013 Elsevier Ltd and INQUA. All rights reserved.
Tile drains strongly influence the water cycle in agricultural catchment in terms of water quantity and quality. The connectivity of preferential flow to tile drains can create shortcuts for rapid transport of solutes into surface waters. The leaching of pesticides can be linked to a set of main factors including, rainfall characteristics, soil moisture, chemical properties of the pesticides, soil properties, and preferential flow paths. The connectivity of the macropore system to the tile drain is crucial for pesticide leaching. Concurring influences of the main factors, threshold responses and the role of flow paths are still poorly understood. The objective of this study is to investigate these influences by a replica series of three irrigation experiments on a tile drain field site using natural and artificial tracers together with applied pesticides. We found a clear threshold behavior in the initialization of pesticide transport that was different between the replica experiments. Pre-event soil water contributed significantly to the tile drain flow, and creates a flow path for stored pesticides from the soil matrix to the tile drain. This threshold is controlled by antecedent soil moisture and precipitation characteristics, and the interaction between the soil matrix and preferential flow system. Fast transport of pesticides without retardation and the remobilization could be attributed to this threshold and the interaction between the soil matrix and the preferential flow system. Thus, understanding of the detailed preferential flow processes clearly enhances the understanding of pesticide leaching on event and long term scale, and can further improve risk assessment and modeling approaches. (C) 2014 Elsevier B.V. All rights reserved.
The Mohorovicic discontinuity, Moho for short, which marks the boundary between crust and mantle, is the main first-order structure within the lithosphere. Geodynamics and tectonic evolution determine its depth level and properties. Here, we present a map of the Moho in central Europe across the Teisseyre-Tornquist Zone, a region for which a number of previous studies are available. Our results are based on homogeneous and consistent processing of P-and S-receiver functions for the largest passive seismological data set in this region yet, consisting of more than 40 000 receiver functions from almost 500 station. Besides, we also provide new results for the crustal vP/vS ratio for the whole area.
Our results are in good agreement with previous, more localized receiver function studies, as well as with the interpretation of seismic profiles, while at the same time resolving a higher level of detail than previous maps covering the area, for example regarding the Eifel Plume region, Rhine Graben and northern Alps. The close correspondence with the seismic data regarding crustal structure also increases confidence in use of the data in crustal corrections and the imaging of deeper structure, for which no independent seismic information is available.
In addition to the pronounced, stepwise transition from crustal thicknesses of 30 km in Phanerozoic Europe to more than 45 beneath the East European Craton, we can distinguish other terrane boundaries based on Moho depth as well as average crustal v(P)/v(S) ratio and Moho phase amplitudes. The terranes with distinct crustal properties span a wide range of ages, from Palaeoproterozoic in Lithuania to Cenozoic in the Alps, reflecting the complex tectonic history of Europe. Crustal thickness and properties in the study area are also markedly influenced by tectonic overprinting, for example the formation of the Central European Basin System, and the European Cenozoic Rift System. In the areas affected by Cenozoic rifting and volcanism, thinning of the crust corresponds to lithospheric updoming reported in recent surface wave and S-receiver function studies, as expected for thermally induced deformation. The same correlation applies for crustal thickening, not only across the Trans-European Suture Zone, but also within the southern part of the Bohemian Massif.
The paper examines the quality of satellite-abased precipitation estimates for the lower Mahanadi River basin (eastern India). The considered data sets known as 3B42 and 3B42-RT (version 7/7A) are routinely produced by the tropical rainfall measuring mission (TRMM) from passive microwave and infrared recordings. While the 3B42-RT data are disseminated in real time, the gauge-aadjusted 3B42 data set is published with a delay of some months. The quality of the two products was assessed in a two-astep procedure. First, the correspondence between the remotely sensed precipitation rates and rain gauge data was evaluated at the subbasin scale. Second, the quality of the rainfall estimates was assessed by analysing their performance in the context of rainfall-arunoff simulation.
At sub-abasin level (4000 to 16 000 km(2)) the satellite-abased areal precipitation estimates were found to be moderately correlated with the gauge-abased counterparts (R-2 of 0.64-0.74 for 3B42 and 0.59-0.72 for 3B42-RT). Significant discrepancies between TRMM data and ground observations were identified at high-aintensity levels. The rainfall depth derived from rain gauge data is often not reflected by the TRMM estimates (hit rate < 0.6 for ground-abased intensities > 80 mm day(-1)). At the same time, the remotely sensed rainfall rates frequently exceed the gauge-abased equivalents (false alarm ratios of 0.2-0.6). In addition, the real-atime product 3B42-RT was found to suffer from a spatially consistent negative bias.
Since the regionalisation of rain gauge data is potentially associated with a number of errors, the above results are subject to uncertainty. Hence, a validation against independent information, such as stream flow, was essential. In this case study, the outcome of rainfall-arunoff simulation experiments was consistent with the above-mentioned findings. The best fit between observed and simulated stream flow was obtained if rain gauge data were used as model input (Nash-Sutcliffe index of 0.76-0.88 at gauges not affected by reservoir operation). This compares to the values of 0.71-0.78 for the gauge-djusted TRMM 3B42 data and 0.65-0.77 for the 3B42-RT real-atime data. Whether the 3B42-RT data are useful in the context of operational runoff prediction in spite of the identified problems remains a question for further research.
The goal of China's sloping land conversion programme (SLCP) is to combat soil erosion and to reduce rural poverty. An ex-ante assessment of possible SLCP impacts was conducted with a focus on rural sustainability, taking the drought-prone region of Guyuan in Western China as an example. The Framework for Participatory Impact Assessment (FoPIA) was used to conduct two complementary impact assessments, one assessing SLCP impacts at regional level and a second one assessing alternative forest management options, to explore possible trade-offs among the economic, social and environmental dimensions of sustainability. Regional stakeholders assessed the SLCP to be capable of reducing soil erosion but felt it negatively affected rural employment, and a further continuation of the Programme was advocated. Assessment of three forest management scenarios by scientists showed that an orientation towards energy forests is potentially beneficial to all three sustainability dimensions. Ecological forests had disproportionate positive impacts on environmental functions and adverse impact on the other two sustainability dimensions. Economic forests were assessed to serve primarily the economic and social sustainability dimensions, while environmental impacts were still tolerable. The FoPIA results were evaluated against the available literature on the SLCP. Overall, the assessment results appeared to be reasonable, but the results of the regional stakeholders appeared to be too optimistic compared with the more critical assessment of the scientists. The SLCP seems to have the potential to tackle soil erosion but requires integrated forest management to minimize the risk of water stress while contributing to economic and social benefits in Guyuan. Copyright (C) 2012 John Wiley & Sons, Ltd.
Quantitative geochemical modeling is today applied in a variety of geological environments from the petrogenesis of igneous rocks to radioactive waste disposal. In addition, the development of thermodynamic databases and computer programs to calculate equilibrium phase diagrams has greatly advanced our ability to model geodynamic processes. Combined with experimental data on elemental partitioning and isotopic fractionation, thermodynamic forward modeling unfolds enormous capacities that are far from exhausted.
In metamorphic petrology the combination of thermodynamic and trace element forward modeling can be used to study and to quantify processes at spatial scales from mu m to km. The thermodynamic forward models utilize Gibbs energy minimization to quantify mineralogical changes along a reaction path of a chemically open fluid/rock system. These results are combined with mass balanced trace element calculations to determine the trace element distribution between rock and melt/fluid during the metamorphic evolution. Thus, effects of mineral reactions, fluid-rock interaction and element transport in metamorphic rocks on the trace element and isotopic composition of minerals, rocks and percolating fluids or melts can be predicted.
Here we illustrate the capacities of combined thermodynamic-geochemical modeling based on two examples relevant to mass transfer during metamorphism. The first example focuses on fluid-rock interaction in and around a blueschist-facies shear zone in felsic gneisses, where fluid-induced mineral reactions and their effects on boron (B) concentrations and isotopic compositions in white mica are modeled. In the second example, fluid release from a subducted slab, the associated transport of B as well as variations in B concentrations and isotopic compositions in liberated fluids and residual rocks are modeled. We compare the modeled results of both examples to geochemical data of natural minerals and rocks and demonstrate that the combination of thermodynamic and geochemical models enables quantification of metamorphic processes and insights into element cycling that would have been unattainable if only one model approach was chosen. (C) 2014 Elsevier B.V. All rights reserved.
Quantifying volumes and rates of delivery of terrestrial sediment from island arcs to subduction zones is indispensable for refining estimates of the thickness of trench fills that may eventually control the location and timing of submarine landslides and tsunami-generating mega-earthquakes. Despite these motivating insights, knowledge about the rates of erosion and sediment export from the Japanese islands to their Pacific subduction zones remains patchy regardless of the increasing availability of highly resolved data on surface deformation, climate, geology, and topography. Traditionally, natural erosion rates across the island arc have been estimated from regression of topographic catchment metrics and reservoir sedimentation rates that were recorded over several years to decades. We review current research in this context, correct for a systematic bias in one of the most widely used predictions, and present new estimates of decadal to millennial-scale erosion rates of Japan's terrestrial inner forearc. We draw on several independent and unprecedented inventories of mass wasting, reservoir sedimentation, and concentrations of cosmogenic Be-10 in river sands. We find that natural Be-10-derived denudation rates of several mm yr(-1) in the Japanese Alps have been sustained over several centuries to millennia, and are, within error, roughly consistent with sediment yields inferred from artificial reservoir sedimentation. Local exceptions may likely result from release of sediment storage or regional landsliding episodes that trigger transient sediment pulses. Our synopsis further reveals that catchments draining Japan's eastern seaboard differ distinctly in their tectonic, lithological, topographic, and climatic characteristics between the Tohoku, Japanese Alps, and Nankai inner forearc segments, which is underscored by a marked asymmetric pattern of erosion rates along the island arc. Erosion rates are highest (up to at least 3 mm yr(-1)) in the Japanese Alps that mark the collision of two subduction zones, where high topographic relief, hillslope and bedrock-channel steepness foster rapid denudation by mass wasting. Comparable, if slightly lower, erosion rates characterise the Nankai inner forearc in southwest Japan, most likely due to higher typhoon-driven rainfall totals and variability rather than its high topographic relief. In contrast, our estimated erosion and flux rates are lowest in the Tohoku inner forearc catchments that feed sediment into the Japan Trench. We conclude that collisional mountain building of the Japanese Alps drives some of the highest erosion rates in the island arc despite similar uplift and precipitation controls in southwest Japan. We infer that, prior to extensive river damming, reservoir construction, and coastal works, the gross of Japan's total sediment export to the Pacific Ocean entered the accretionary margin of the Nankai Trough as opposed to the comparatively sediment-starved Japan Trench. Compared to documented contemporary rates of sediment flux from mountainous catchments elsewhere in the Pacific, the rivers draining Japan's inner forearc take an intermediate position despite high relief, steep slopes, very high seismicity, and frequent rainstorms. However, the average rates of millennial-scale denudation in the Japanese Alps particularly are amongst the highest reported worldwide.
Local mismatches between these late Holocene and modern rates emphasise the anthropogenic fingerprint on sediment retention that may have significantly reduced the island arc's mass flux to its subduction zones, as is the case elsewhere in east and southeast Asia. (C) 2014 Elsevier B.V. All rights reserved.
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.
Costing natural hazards
(2014)
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 geological feature in the Qaidam Basin known as the "Shell Bar" contains millions of freshwater clam shells buried in situ. Since the 1980s, this feature in the now hyper-arid basin has been interpreted to be lake deposits that provide evidence for a warmer and more humid climate than present during late marine isotope stage 3 (MIS 3). Global climate during late MIS 3 and the last glacial maximum, however, was cold and dry, with much lower sea levels. We re-investigated the feature geomorphologically and sedimentologically, and employed optically stimulated luminescence (OSL) dating to verify the chronology of the sediments. We interpret the Shell Bar to be a remnant of a river channel formed by a stream that ran across an exposed lake bed during a regressive lake phase. Deflation of the surrounding older, fine-grained lacustrine deposits has left the fluvial channel sediments topographically inverted, indicating the erosive nature of the landscape. Luminescence ages place the formation of the Shell Bar in MIS 5 (similar to 113-99 ka), much older than previous radiocarbon ages of < 40 ka BP, but place the paleoclimatic inferences more in accord with other regional and global climate proxy records. We present a brief review of the age differences derived from C-14 and OSL dating of some critical sections that were thought to represent a warmer and more humid climate than present during late MIS 3. We attribute the differences to underestimation of C-14 ages. We suggest that C-14 ages older than similar to 25 ka BP may require re-investigation, especially dates on samples from arid regions.
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.
We focus on the relation between seismic and total postseismic afterslip following the Maule M-w 8.8 earthquake on 2010 February 27 in central Chile. First, we calculate the cumulative slip released by aftershock seismicity. We do this by summing up the aftershock regions and slip estimated from scaling relations. Comparing the cumulative seismic slip with afterslip modelswe showthat seismic slip of individual aftershocks exceeds locally the inverted afterslip model from geodetic constraints. As the afterslip model implicitly contains the displacements from the aftershocks, this reflects the tendency of afterslip models to smear out the actual slip pattern. However, it also suggests that locally slip for a number of the larger aftershocks exceeds the aseismic slip in spite of the fact that the total equivalent moment of the afterslip exceeds the cumulative moment of aftershocks by a large factor. This effect, seen weakly for the Maule 2010 and also for the Tohoku 2011 earthquake, can be explained by taking into account the uncertainties of the seismicity and afterslip models. In spite of uncertainties, the hypocentral region of the Nias 2005 earthquake is suggested to release a large fraction of moment almost purely seismically. Therefore, these aftershocks are not driven solely by the afterslip but instead their slip areas have probably been stressed by interseismic loading and the mainshock rupture. In a second step, we divide the megathrust of the Maule 2010 rupture into discrete cells and count the number of aftershocks that occur within 50 km of the centre of each cell as a function of time. We then compare this number to a time-dependent afterslip model by defining the 'afterslip to aftershock ratio' (ASAR) for each cell as the slope of the best fitting line when the afterslip at time t is plotted against aftershock count. Although we find a linear relation between afterslip and aftershocks for most cells, there is significant variability in ASAR in both the downdip and along-strike directions of the megathrust. We compare the spatial distribution of ASAR with the spatial distribution of seismic coupling, coseismic slip and Bouguer gravity anomaly, and in each case we find no significant correlation.
Advance in geocomputation
(2014)
In general, a moderate drying trend is observed in mid-latitude arid Central Asia since the Mid-Holocene, attributed to the progressively weakening influence of the mid-latitude Westerlies on regional climate. However, as the spatio-temporal pattern of this development and the underlying climatic mechanisms are yet not fully understood, new high-resolution paleoclimate records from this region are needed. Within this study, a sediment core from Lake Son Kol (Central Kyrgyzstan) was investigated using sedimentological, (bio) geochemical, isotopic, and palynological analyses, aiming at reconstructing regional climate development during the last 6000 years. Biogeochemical data, mainly reflecting summer moisture conditions, indicate predominantly wet conditions until 4950 cal. yr BP, succeeded by a pronounced dry interval between 4950 and 3900 cal. yr BP. In the following, a return to wet conditions and a subsequent moderate drying trend until present times are observed. This is consistent with other regional paleoclimate records and likely reflects the gradual Late Holocene diminishment of the amount of summer moisture provided by the mid-latitude Westerlies. However, climate impact of the Westerlies was apparently not only restricted to the summer season but also significant during winter as indicated by recurrent episodes of enhanced allochthonous input through snowmelt, occurring before 6000 cal. yr BP and at 5100-4350, 3450-2850, and 1900-1500 cal. yr BP. The distinct similar to 1500year periodicity of these episodes of increased winter precipitation in Central Kyrgyzstan resembles similar cyclicities observed in paleoclimate records around the North Atlantic, likely indicating a hemispheric-scale climatic teleconnection and an impact of North Atlantic Oscillation (NAO) variability in Central Asia.
Magmatism forming the Central Anatolian Volcanic Province of Cappadocia, central Turkey, records the last phase of Neotethyan subduction after similar to 11 Ma. Thirteen large calc-alkaline ignimbrite sheets form marker bands within the volcano-sedimentary succession (the Urgup Formation) and provide a robust chronostratigraphy for paleoecologic evaluation of the interleaved paleosols. This paper evaluates the chronologic record in the context of the radiometric, magnetostratigraphic and lithostratigraphic controls. Previous inconsistencies relating primarily to K/Ar evidence were reason for the initiation of an integrated study which includes Ar-40/Ar-39 dating, palaeomagnetic and stratigraphic evidence. The newly determined Ar-40/Ar-39-ages (Lepetit, 2010) are in agreement with Ar/Ar and U/Pb data meanwhile published by Pauquette and Le Pennec (2012) and Aydar et al. (2012). The Ar-40/Ar-39-ages restrict the end of the Urgup Formation to the late Miocene. The paleosol sequence enclosed by the ignimbrites is thus restricted to the late Miocene, the most intense formation of pedogene calcretes correlating with the Messinian Salinity Crisis.
In aquatic systems, natural organic matter (NOM) and in particular humic substances effectively absorb the ultraviolet (UV)/visible light spectrum of solar radiation and act as a photoprotective filter for organisms. Simultaneously, UV contributes to the generation of potentially harmful reactive oxygen species (ROS). Dose-response experiments were conducted on cyanobacteria and green algae with hydrogen peroxide (H2O2) as a long-lived representative of ROS. Delayed fluorescence (DF) decay kinetics was used as a non-invasive tool to follow changes of phytoplankton activity in real time. In order to investigate phototoxicity and photoprotection by NOM on phytoplankton, we exposed algae to UV-pre-irradiated NOM and direct UV excitation. Cyanobacteria responded to H2O2 concentrations as low as 10(-7) M, while green algae were 2 orders of magnitude less sensitive. UV irradiation of medium with NOM generated H2O2 concentrations of 1.5 x 10(-7) to 3.6 x 10(-7) M. When exposed to these concentrations, only the DF of cyanobacteria led to a measurable effect while that of green algae did not change. The addition of NOM protected all phytoplankton from direct UV irradiation, but cyanobacteria benefitted less. From this we conclude that UV-irradiated water enriched with NOM can adversely affect the physiology of cyanobacteria, but not of green algae, which might control phytoplankton composition and species-specific activities.
During the last glacial period lake Kinneret (the Sea of Galilee) fluctuated between high and low water levels reflecting the hydrological conditions of the lake watershed. Here, we focus on the hydrology of the lake after its retreat from the last glacial MIS2 (similar to 27-25 ka BP) highest stand of similar to 170 m below mean sea level (m bsl) to the low stand of similar to 214 m bsl at similar to 24-21 ka BP. The limnological-hydrological history of this time interval is recovered from trench and borehole that were dug and drilled in the southwestern shore of the lake at Ohalo-II archeological site. Cyprideis torosa (Ostracoda) recovered from the trench yielded elemental, Sr-87/Sr-86 and delta O-18 isotope data that provide information on the shore environment during the low stand period. The Sr-87/Sr-86 and Sr/Ca ratios in the ostracods, varying between similar to 0.70789 and similar to 0.70815 and 0.0017 and 0.0030, respectively indicate contributions of waters from the last glacial lake and regional runoff. The increase in the Sr-87/Sr-86 ratios reflects the decreasing effect of the last glacial Lake Kinneret waters and enhanced contribution of local runoff that washed down dried mountain soils that were previously developed during the wet and vegetated glacial. The lake retreat at similar to 24 ka BP coincided with the Heinrich event H2 at the northern Atlantic. H2 was expressed by severe aridity in Lake Kinneret-Dead Sea watershed. The limnological-hydrological change at post H2 was accompanied by 2% decrease in the delta O-18 value from -2% to -4% reflecting the change in the composition of the east Mediterranean rain sources.
The last glacial lake Sr-87/Sr-86 ratio is similar to the Tiberias Spa saline waters and distinctly different from the modern Lake Kinneret fresh waters: Sr-87/Sr-86 similar to 0.70785 compared to similar to 0.70760, respectively. This difference is explained by enhanced contribution of Ca-chloride brines with high Sr-87/Sr-86 values to the last glacial lake and reduced Jordan River contribution due to cold freezing conditions at its headwaters, while the modern Lake Kinneret is more affected by low Sr-87/Sr-86 freshwater from the Jordan watershed.
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.
The strong present-day Asian monsoons are thought to have originated between 25 and 22 million years (Myr) ago, driven by Tibetan-Himalayan uplift. However, the existence of older Asian monsoons and their response to enhanced greenhouse conditions such as those in the Eocene period (55-34Myrago) are unknown because of the paucity of well-dated records. Here we show late Eocene climate records revealing marked monsoon-like patterns in rainfall and wind south and north of the Tibetan-Himalayan orogen. This is indicated by low oxygen isotope values with strong seasonality in gastropod shells and mammal teeth from Myanmar, and by aeolian dust deposition in northwest China. Our climate simulations support modern-like Eocene monsoonal rainfall and show that a reinforced hydrological cycle responding to enhanced greenhouse conditions counterbalanced the negative effect of lower Tibetan relief on precipitation. These strong monsoons later weakened with the global shift to icehouse conditions 34 Myr ago.