Refine
Year of publication
- 2024 (4)
- 2023 (14)
- 2022 (68)
- 2021 (74)
- 2020 (97)
- 2019 (75)
- 2018 (82)
- 2017 (38)
- 2016 (9)
- 2015 (5)
- 2014 (15)
- 2013 (34)
- 2012 (21)
- 2011 (10)
- 2010 (6)
- 2009 (30)
- 2008 (19)
- 2007 (11)
- 2006 (24)
- 2005 (15)
- 2004 (20)
- 2003 (21)
- 2002 (13)
- 2001 (20)
- 2000 (18)
- 1999 (3)
- 1998 (10)
- 1997 (2)
- 1996 (2)
- 1995 (6)
- 1994 (2)
- 1993 (1)
- 1992 (2)
Document Type
- Article (537)
- Doctoral Thesis (116)
- Postprint (67)
- Monograph/Edited Volume (21)
- Other (9)
- Review (7)
- Habilitation Thesis (5)
- Master's Thesis (3)
- Part of a Book (2)
- Bachelor Thesis (1)
Language
- English (771) (remove)
Keywords
- climate change (35)
- Curriculum Framework (34)
- European values education (34)
- Europäische Werteerziehung (34)
- Lehrevaluation (34)
- Studierendenaustausch (34)
- Unterrichtseinheiten (34)
- curriculum framework (34)
- lesson evaluation (34)
- student exchange (34)
Institute
- Institut für Umweltwissenschaften und Geographie (771) (remove)
A conundrum of trends
(2022)
This comment is meant to reiterate two warnings: One applies to the uncritical use of ready-made (openly available) program packages, and one to the estimation of trends in serially correlated time series. Both warnings apply to the recent publication of Lischeid et al. about lake-level trends in Germany.
Indices of oscillatory behavior are conveniently obtained by projecting the fields in question into a phase space of a few (mostly just two) dimensions; empirical orthogonal functions (EOFs) or other, more dynamical, modes are typically used for the projection. If sufficiently coherent and in quadrature, the projected variables simply describe a rotating vector in the phase space, which then serves as the basis for predictions. Using the boreal summer intraseasonal oscillation (BSISO) as a test case, an alternative procedure is introduced: it augments the original fields with their Hilbert transform (HT) to form a complex series and projects it onto its (single) dominant EOF. The real and imaginary parts of the corresponding complex pattern and index are compared with those of the original (real) EOF. The new index explains slightly less variance of the physical fields than the original, but it is much more coherent, partly from its use of future information by the HT. Because the latter is in the way of real-time monitoring, the index can only be used in cases with predicted physical fields, for which it promises to be superior. By developing a causal approximation of the HT, a real-time variant of the index is obtained whose coherency is comparable to the noncausal version, but with smaller explained variance of the physical fields. In test cases the new index compares well to other indices of BSISO. The potential for using both indices as an alternative is discussed.
For three small, mountainous catchments in Germany two medium-range forecast systems are compared that predict precipitation for up to 5 days in advance. One system is composed of the global German weather service (DWD) model, GME, which is dynamically downscaled using the COSMO-EU regional model. The other system is an empirical (expanded) downscaling of the ECMWF model IFS. Forecasts are verified against multi-year daily observations, by applying standard skill scores to events of specified intensity. All event classes are skillfully predicted by the empirical system for up to five days lead time. For the available prediction range of one to two days it is superior to the dynamical system.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
A digital filter is introduced which treats the problem of predictability versus time averaging in a continuous, seamless manner. This seamless filter (SF) is characterized by a unique smoothing rule that determines the strength of smoothing in dependence on lead time. The rule needs to be specified beforehand, either by expert knowledge or by user demand. As a result, skill curves are obtained that allow a predictability assessment across a whole range of time-scales, from daily to seasonal, in a uniform manner. The SF is applied to downscaled SEAS5 ensemble forecasts for two focus regions in or near the tropical belt, the river basins of the Karun in Iran and the Sao Francisco in Brazil. Both are characterized by strong seasonality and semi-aridity, so that predictability across various time-scales is in high demand. Among other things, it is found that from the start of the water year (autumn), areal precipitation is predictable with good skill for the Karun basin two and a half months ahead; for the Sao Francisco it is only one month, longer-term prediction skill is just above the critical level.
This study follows up on a previous downscaling intercomparison for present climate. Using a larger set of eight methods the authors downscale atmospheric fields representing present (1981-2000) and future (2046-65) conditions, as simulated by six global climate models following three emission scenarios. Local extremes were studied at 20 locations in British Columbia as measured by the same set of 27 indices, ClimDEX, as in the precursor study. Present and future simulations give 2 x 3 x 6 x 8 x 20 x 27 = 155 520 index climatologies whose analysis in terms of mean change and variation is the purpose of this study. The mean change generally reinforces what is to be expected in a warmer climate: that extreme cold events become less frequent and extreme warm events become more frequent, and that there are signs of more frequent precipitation extremes. There is considerable variation, however, about this tendency, caused by the influence of scenario, climate model, downscaling method, and location. This is analyzed using standard statistical techniques such as analysis of variance and multidimensional scaling, along with an assessment of the influence of each modeling component on the overall variation of the simulated change. It is found that downscaling generally has the strongest influence, followed by climate model; location and scenario have only a minor influence. The influence of downscaling could be traced back in part to various issues related to the methods, such as the quality of simulated variability or the dependence on predictors. Using only methods validated in the precursor study considerably reduced the influence of downscaling, underpinning the general need for method verification.
Several overlapping crises which affected the EU during the past ten years have recently aggravated. Especially the progressing refugee crisis, the persisting financial crisis and geopolitical turmoil in the EU's neighbourhood contributed to the rise of anti-EU movements and diverse articulations of Euroscepticism. Although public opinion and mainstream political analysis have easily identified right-wing populism as one of the most important drivers, it is still doubtful if it can be equated with Euroscepticism without further ado. To date it is by no means clear how and where Euroscepticism exactly originates.
Soil is today considered a non-renewable resource on societal time scale, as the rate of soil loss is higher than the one of soil formation.
Soil formation is complex, can take several thousands of years and is influenced by a variety of factors, one of them is time. Oftentimes, there is the assumption of constant and progressive conditions for soil and/or profile development (i.e., steady-state). In reality, for most of the soils, their (co-)evolution leads to a complex and irregular soil development in time and space characterised by “progressive” and “regressive” phases.
Lateral transport of soil material (i.e., soil erosion) is one of the principal processes shaping the land surface and soil profile during “regressive” phases and one of the major environmental problems the world faces.
Anthropogenic activities like agriculture can exacerbate soil erosion. Thus, it is of vital importance to distinguish short-term soil redistribution rates (i.e., within decades) influenced by human activities differ from long-term natural rates. To do so, soil erosion (and denudation) rates can be determined by using a set of isotope methods that cover different time scales at landscape level.
With the aim to unravel the co-evolution of weathering, soil profile development and lateral redistribution on a landscape level, we used Pluthonium-239+240 (239+240Pu), Beryllium-10 (10Be, in situ and meteoric) and Radiocarbon (14C) to calculate short- and long-term erosion rates in two settings, i.e., a natural and an anthropogenic environment in the hummocky ground moraine landscape of the Uckermark, North-eastern Germany. The main research questions were:
1. How do long-term and short-term rates of soil redistributing processes differ?
2. Are rates calculated from in situ 10Be comparable to those of using meteoric 10Be?
3. How do soil redistribution rates (short- and long-term) in an agricultural and in a natural landscape compare to each other?
4. Are the soil patterns observed in northern Germany purely a result of past events (natural and/or anthropogenic) or are they imbedded in ongoing processes?
Erosion and deposition are reflected in a catena of soil profiles with no or almost no erosion on flat positions (hilltop), strong erosion on the mid-slope and accumulation of soil material at the toeslope position. These three characteristic process domains were chosen within the CarboZALF-D experimental site, characterised by intense anthropogenic activities. Likewise, a hydrosequence in an ancient forest was chosen for this study and being regarded as a catena strongly influenced by natural soil transport.
The following main results were obtained using the above-mentioned range of isotope methods available to measure soil redistribution rates depending on the time scale needed (e.g., 239+240Pu, 10Be, 14C):
1. Short-term erosion rates are one order of magnitude higher than long-term rates in agricultural settings.
2. Both meteoric and in situ 10Be are suitable soil tracers to measure the long-term soil redistribution rates giving similar results in an anthropogenic environment for different landscape positions (e.g., hilltop, mid-slope, toeslope)
3. Short-term rates were extremely low/negligible in a natural landscape and very high in an agricultural landscape – -0.01 t ha-1 yr-1 (average value) and -25 t ha-1 yr-1 respectively. On the contrary, long-term rates in the forested landscape are comparable to those calculated in the agricultural area investigated with average values of -1.00 t ha-1 yr-1 and -0.79 t ha-1 yr-1.
4. Soil patterns observed in the forest might be due to human impact and activities started after the first settlements in the region, earlier than previously postulated, between 4.5 and 6.8 kyr BP, and not a result of recent soil erosion.
5. Furthermore, long-term soil redistribution rates are similar independently from the settings, meaning past natural soil mass redistribution processes still overshadow the present anthropogenic erosion processes.
Overall, this study could make important contributions to the deciphering of the co-evolution of weathering, soil profile development and lateral redistribution in North-eastern Germany. The multi-methodological approach used can be challenged by the application in a wider range of landscapes and geographic regions.
Traditional drainage systems combining man-made channels and subsurface tile drains have been used since Roman times to control water excess in Mediterranean lowland regions, favouring adequate soil water regime for agriculture purposes. However, mechanization of agriculture, abandonment or land use changes lead to a progressive deterioration of these drains in the last decades. The effects of these structures on hydrological and sediment dynamics have been previously analyzed in a small Mediterranean lowland catchment (Can Revull, Mallorca, Spain, 1.4 km2) by establishing an integrated sediment budget with a multi-technique approach. Moreover, the recent advances in morphometric techniques enable the completion of this analysis by the accurate identification of active areas (i.e. sources, pathway links, and sinks) and improve the understanding of (de-)coupling mechanisms of water and sediment linkages. In this study, the Borselli's index of connectivity (IC; Cavalli et al. (2013)'s version) derived from a LiDAR-based high resolution DEM (>1 pt m−2; RMSE < 0.2 m) was used to evaluate the spatial patterns of sediment connectivity of the catchment under two different scenarios: (1) the current scenario, including an accurate representation of the 3800 m of artificial channels and levees (CS - Channelled Scenario), and (2) a hypothetical scenario in which these anthropogenic features were removed (US - Unchannelled Scenario). Design and configuration of the drainage system in Can Revull generated changes favouring lateral decoupling between different compartments, with hillslopes-floodplain and floodplain-channels relationships, showing a general decrease of IC values, and high longitudinal connectivity along the artificial channel network. Field observations corroborated these results: structures enabled rapid drainage of the water excess also promoting low surface runoff within the field crops, proving to be an effective management practice for erosion control in agricultural Mediterranean lowland catchments. By contrast, US demonstrated that the abandonment of the current agricultural practices and the subsequent destruction of the drainage system could lead the higher soil loss rates owning to more intense/effective processes of sediment connectivity.
This study investigates spatial patterns of Ks and tests the hypothesis of whether structural variance emerges from noise with increasing sampling precision. We analyzed point measurements of Ks along independent transects at sampling intervals of 25, 10, 1 and 0.25 m. The field area is a tropical rainforest catena (i.e. toposequence) characterized by systematic downslope changes in soil properties including color (red to yellow), mineralogy (kaolinite- illite to kaolinite) and texture (sandy clay to sand). Independent tramsects spanning the entire catena at lag intervals of 25 and 10 in reveal little to no spatial patterns in Ks; i.e. scatter plots are noisy and lack apparent spatial trends, and semivariograms suggest little to no autocorrelation in Ks. As sampling precision is increased (h = 1 and 0.25 m), spatial patterns emerge in Ks for the downslope areas, in which distinctive hydraulic boundaries in Ks correlate with relatively small-scale, topography-controlled soils with coarse textures (greater than or equal to 80% sand). For these areas, semivariograms of Ks and those of %sand and %clay exhibit similar spatial structure characterized by small nugget variances and large ranges, and nugget variance is reduced as sampling precision increases from 1 to 0.25 m. In the upslope, clay-rich locations along this toposequence, Ks exhibits few spatial patterns, irrespective of sampling scale. For these locations, scatter plots are noisy without apparent spatial trends, and semivariograms show almost complete nugget variance, suggesting little to no correlation in this hydraulic parameter at any scale. This study suggests that in the absence of coarse textures (greater than or equal to 80% sand), there is little predictability in Ks, even at sampling intervals of 0.25 m. We believe this lack of spatial structure is due to a predominance of small-scale processes such as biological activity that largely control Ks in this forested setting. (C) 2003 Elsevier B.V. All rights reserved
Flood loss modeling is an important component within flood risk assessments. Traditionally, stage-damage functions are used for the estimation of direct monetary damage to buildings. Although it is known that such functions are governed by large uncertainties, they are commonly applied - even in different geographical regions - without further validation, mainly due to the lack of real damage data. Until now, little research has been done to investigate the applicability and transferability of such damage models to other regions. In this study, the last severe flood event in the Austrian Lech Valley in 2005 was simulated to test the performance of various damage functions from different geographical regions in Central Europe for the residential sector. In addition to common stage-damage curves, new functions were derived from empirical flood loss data collected in the aftermath of recent flood events in neighboring Germany. Furthermore, a multi-parameter flood loss model for the residential sector was adapted to the study area and also evaluated with official damage data. The analysis reveals that flood loss functions derived from related and more similar regions perform considerably better than those from more heterogeneous data sets of different regions and flood events. While former loss functions estimate the observed damage well, the latter overestimate the reported loss clearly. To illustrate the effect of model choice on the resulting uncertainty of damage estimates, the current flood risk for residential areas was calculated. In the case of extreme events like the 300 yr flood, for example, the range of losses to residential buildings between the highest and the lowest estimates amounts to a factor of 18, in contrast to properly validated models with a factor of 2.3. Even if the risk analysis is only performed for residential areas, our results reveal evidently that a carefree model transfer in other geographical regions might be critical. Therefore, we conclude that loss models should at least be selected or derived from related regions with similar flood and building characteristics, as far as no model validation is possible. To further increase the general reliability of flood loss assessment in the future, more loss data and more comprehensive loss data for model development and validation are needed.
Fires are a fundamental part of the Earth System. In the last decades, they have been altering ecosystem structure, biogeochemical cycles and atmospheric composition with unprecedented rapidity. In this study, we implement a complex networks-based methodology to track individual fires over space and time. We focus on extreme fires-the 5% most intense fires-in the tropical forests of the Brazilian Legal Amazon over the period 2002-2019. We analyse the interannual variability in the number and spatial patterns of extreme forest fires in years with diverse climatic conditions and anthropogenic pressure to examine potential synergies between climate and anthropogenic drivers. We observe that major droughts, that increase forest flammability, co-occur with high extreme fire years but also that it is fundamental to consider anthropogenic activities to understand the distribution of extreme fires. Deforestation fires, fires escaping from managed lands, and other types of forest degradation and fragmentation provide the ignition sources for fires to ignite in the forests. We find that all extreme forest fires identified are located within a 0.5-km distance from forest edges, and up to 56% of them are within a 1-km distance from roads (which increases to 73% within 5 km), showing a strong correlation that defines spatial patterns of extreme fires.
Pollen records from Siberia are mostly absent in global or Northern Hemisphere synthesis works. Here we present a taxonomically harmonized and temporally standardized pollen dataset that was synthesized using 173 palynological records from Siberia and adjacent areas (northeastern Asia, 42-75 degrees N, 50-180 degrees E). Pollen data were taxonomically harmonized, i.e. the original 437 taxa were assigned to 106 combined pollen taxa. Age-depth models for all records were revised by applying a constant Bayesian age-depth modelling routine. The pollen dataset is available as count data and percentage data in a table format (taxa vs. samples), with age information for each sample. The dataset has relatively few sites covering the last glacial period between 40 and 11.5 ka (calibrated thousands of years before 1950 CE) particularly from the central and western part of the study area. In the Holocene period, the dataset has many sites from most of the area, with the exception of the central part of Siberia. Of the 173 pollen records, 81 % of pollen counts were downloaded from open databases (GPD, EPD, PANGAEA) and 10 % were contributions by the original data gatherers, while a few were digitized from publications. Most of the pollen records originate from peatlands (48 %) and lake sediments (33 %). Most of the records (83 %) have >= 3 dates, allowing the establishment of reliable chronologies. The dataset can be used for various purposes, including pollen data mapping (example maps for Larix at selected time slices are shown) as well as quantitative climate and vegetation reconstructions. The datasets for pollen counts and pollen percentages are available at https://doi.org/10.1594/PANGAEA.898616 (Cao et al., 2019a), also including the site information, data source, original publication, dating data, and the plant functional type for each pollen taxa.
Pollen records from Siberia are mostly absent in global or Northern Hemisphere synthesis works. Here we present a taxonomically harmonized and temporally standardized pollen dataset that was synthesized using 173 palynological records from Siberia and adjacent areas (northeastern Asia, 42-75 degrees N, 50-180 degrees E). Pollen data were taxonomically harmonized, i.e. the original 437 taxa were assigned to 106 combined pollen taxa. Age-depth models for all records were revised by applying a constant Bayesian age-depth modelling routine. The pollen dataset is available as count data and percentage data in a table format (taxa vs. samples), with age information for each sample. The dataset has relatively few sites covering the last glacial period between 40 and 11.5 ka (calibrated thousands of years before 1950 CE) particularly from the central and western part of the study area. In the Holocene period, the dataset has many sites from most of the area, with the exception of the central part of Siberia. Of the 173 pollen records, 81 % of pollen counts were downloaded from open databases (GPD, EPD, PANGAEA) and 10 % were contributions by the original data gatherers, while a few were digitized from publications. Most of the pollen records originate from peatlands (48 %) and lake sediments (33 %). Most of the records (83 %) have >= 3 dates, allowing the establishment of reliable chronologies. The dataset can be used for various purposes, including pollen data mapping (example maps for Larix at selected time slices are shown) as well as quantitative climate and vegetation reconstructions. The datasets for pollen counts and pollen percentages are available at https://doi.org/10.1594/PANGAEA.898616 (Cao et al., 2019a), also including the site information, data source, original publication, dating data, and the plant functional type for each pollen taxa.
Atmospheric dynamics of extreme discharge events from 1979 to 2016 in the southern Central Andes
(2020)
During the South-American Monsoon season, deep convective systems occur at the eastern flank of the Central Andes leading to heavy rainfall and flooding. We investigate the large- and meso-scale atmospheric dynamics associated with extreme discharge events (> 99.9th percentile) observed in two major river catchments meridionally stretching from humid to semi-arid conditions in the southern Central Andes. Based on daily gauge time series and ERA-Interim reanalysis, we made the following three key observations: (1) for the period 1940-2016 daily discharge exhibits more pronounced variability in the southern, semi-arid than in the northern, humid catchments. This is due to a smaller ratio of discharge magnitudes between intermediate (0.2 year return period) and rare events (20 year return period) in the semi-arid compared to the humid areas; (2) The climatological composites of the 40 largest discharge events showed characteristic atmospheric features of cold surges based on 5-day time-lagged sequences of geopotential height at different levels in the troposphere; (3) A subjective classification revealed that 80% of the 40 largest discharge events are mainly associated with the north-northeastward migration of frontal systems and 2/3 of these are cold fronts, i.e. cold surges. This work highlights the importance of cold surges and their related atmospheric processes for the generation of heavy rainfall events and floods in the southern Central Andes.
Apple replant disease (ARD) is a specific apple-related form of soil fertility loss due to unidentified causes and is also known as soil fatigue. The effect typically appears in monoculture production sites and leads to production decreases of up to 50%, even though the cultivation practice remains the same. However, an indication of replant disease is challenged by the lack of specification of the particular microbial group responsible for ARD. The objective of this study was to establish an algorithm for estimating growth suppression in orchards irrespective of the unknowns in the complex causal relationship by assessing plant-soil interaction in the orchard several years after planting. Based on a comparison between no-replant and replant soils, the Alternaria group (Ag) was identified as a soil-fungal population responding to replant with abundance. The trunk cross-sectional area (CSA) was found to be a practical and robust parameter representing below-ground and above-ground tree performance. Suppression of tree vigour was therefore calculated by dividing the two inversely related parameters, Q = ln(Ag)/CSA, as a function of soil-fungal proportions and plant responses at the single-tree level. On this basis, five clusters of tree vigour suppression (Q) were defined: (1) no tree vigour suppression/vital (0%), (2) escalating (- 38%), (3) strong (- 53%), (4) very strong (- 62%), and (5) critical (- 74%). By calculating Q at the level of the single tree, trees were clustered according to tree vigour suppression. The weighted frequency of clusters in the field allowed replant impact to be quantified at field level. Applied to a case study on sandy brown, dry diluvial soils in Brandenburg, Germany, the calculated tree vigour suppression was 46% compared to the potential tree vigour on no-replant soil in the same field. It is highly likely that the calculated growth suppression corresponds to ARD-impact This result is relevant for identifying functional changes in soil and for monitoring the economic effects of soil fatigue in apple orchards, particularly where long-period crop rotation or plot exchange are improbable.
Apple replant disease (ARD) is a specific apple-related form of soil fertility loss due to unidentified causes and is also known as soil fatigue. The effect typically appears in monoculture production sites and leads to production decreases of up to 50%, even though the cultivation practice remains the same. However, an indication of replant disease is challenged by the lack of specification of the particular microbial group responsible for ARD. The objective of this study was to establish an algorithm for estimating growth suppression in orchards irrespective of the unknowns in the complex causal relationship by assessing plant-soil interaction in the orchard several years after planting. Based on a comparison between no-replant and replant soils, the Alternaria group (Ag) was identified as a soil-fungal population responding to replant with abundance. The trunk cross-sectional area (CSA) was found to be a practical and robust parameter representing below-ground and above-ground tree performance. Suppression of tree vigour was therefore calculated by dividing the two inversely related parameters, Q = ln(Ag)/CSA, as a function of soil-fungal proportions and plant responses at the single-tree level. On this basis, five clusters of tree vigour suppression (Q) were defined: (1) no tree vigour suppression/vital (0%), (2) escalating (- 38%), (3) strong (- 53%), (4) very strong (- 62%), and (5) critical (- 74%). By calculating Q at the level of the single tree, trees were clustered according to tree vigour suppression. The weighted frequency of clusters in the field allowed replant impact to be quantified at field level. Applied to a case study on sandy brown, dry diluvial soils in Brandenburg, Germany, the calculated tree vigour suppression was 46% compared to the potential tree vigour on no-replant soil in the same field. It is highly likely that the calculated growth suppression corresponds to ARD-impact This result is relevant for identifying functional changes in soil and for monitoring the economic effects of soil fatigue in apple orchards, particularly where long-period crop rotation or plot exchange are improbable.
Nitrogen transformations in flowpaths leading from soils to streams in Amazon forest and pasture
(2009)
The paper presents the first study of heavy-mineral sand beaches from the Mediterranean coast of Annaba/Algeria. The studied beaches run along the basement outcrops of the Edough massif, which are mainly composed by micaschists, tourmaline-rich quartzo-feldspathic veins, gneisses, skarns and marbles. Sand samples were taken from three localities (Ain Achir, Plage-Militaire and El Nasr). The heavy-mineral fraction comprises between 74 and 91 vol%. The garnets of the beaches are almandine rich and tourmalines vary with respect to their location from schorl to dravite. Tourmaline at Ain Achir and the Plage-Militaire is schorlits, while at El Nasr beach dravite is ubiquitous. The World Shale Average normalised REE of the sands and the basement outcrops reveal: (i) Ain Achir beach: REE pattern of sand and the coastal rocks from the studied beaches reflects a multiple sources; (ii) Plage-Militaire: the sand and the coastal outcrops show similar LREE and a strong enrichment in HREE, suggesting the presence HREE-rich phases found as inclusions in staurolite; (iii) El Nasr: two types of sand patterns are found: one with flat REE pattern similar to the proximal rocks and other one enriched in HREE suggesting a mixed source.
Oases is a special ecosystem formed in arid climate and hungriness environment, in which resident, water and soil are the principal factor and exchanges of materials, energy and information are the main functional characteristics. The oases regions in central Asia are not only the basilic cradle of civilization of human beings, but also the important strategic places in world growing awareness of the potential benefits. We choose Keriya River Basin oases in south of Xinjiang as a case to study critical controlling of Oases Evolution, Based on the theories and methods used for environmental geology, physical geography, land resource research, and oases ecology. This study try to indicate the essential factors driving the oases ecosystem and the interactional dynamic mechanism in different scales and levels, confirm the optimal equilibrium aggregate of harmonious development between Population, Resources, Environment and Development, and establish the critical controlling pattern of sustainable development. We advance the indicator system to research the evolution of the PRED System of oases in Keriya River valley oases, in basis of the information derived from the field investigation and local materials. According to inquisitional result based on technical support of Geographic Information System (GIS) and Remote Sense (RS), the comparisons and analyses are carried out in land use at the upper reaches, vegetation change in the middle reaches, and desertification at the lower reaches, which narrates the regulations of Keriya River Valley oases land cover dynamic change. The main land cover types represent distinct characteristics of the local place. On the basis of field survey and statistical data, we use ARCINFO software to preprocess these data and the 2 TM satellite images. Through analyzing these images resulting from post- classification compare, we sums up the concrete quantificational dynamic distributed data of 13 land types covering a span of 15 years and regulation of the local ecological environment system. It finally points out that the trend of Keriya River Valley oases desertification expansion is mainly related to two important reasons: impact of natural environment and impact of human activities.
The impact of inclined faults on the hydrothermal field is assessed by adding simplified structural settings to synthetic models. This study is innovative in carrying out numerical simulations because it integrates the real 3-D nature of flow influenced by a fault in a porous medium, thereby providing a useful tool for complex geothermal modelling. The 3-D simulations for the coupled fluid flow and heat transport processes are based on the finite element method. In the model, one geological layer is dissected by a dipping fault. Sensitivity analyses are conducted to quantify the effects of the fault's transmissivity on the fluid flow and thermal field. Different fault models are compared with a model where no fault is present to evaluate the effect of varying fault transmissivity. The results show that faults have a significant impact on the hydrothermal field. Varying either the fault zone width or the fault permeability will result in relevant differences in the pressure, velocity and temperature field. A linear relationship between fault zone width and fluid velocity is found, indicating that velocities increase with decreasing widths. The faults act as preferential pathways for advective heat transport in case of highly transmissive faults, whereas almost no fluid may be transported through poorly transmissive faults.
Relating to the students
(2013)
With the growing size and use of night light time series from the Visible Infrared Imaging Radiometer Suite Day/Night Band (DNB), it is important to understand the stability of the dataset. All satellites observe differences in pixel values during repeat observations. In the case of night light data, these changes can be due to both environmental effects and changes in light emission. Here we examine the stability of individual locations of particular large scale light sources (e.g., airports and prisons) in the monthly composites of DNB data from April 2012 to September 2017. The radiances for individual pixels of most large light emitters are approximately normally distributed, with a standard deviation of typically 15-20% of the mean. Greenhouses and flares, however, are not stable sources. We observe geospatial autocorrelation in the monthly variations for nearby sites, while the correlation for sites separated by large distances is small. This suggests that local factors contribute most to the variation in the pixel radiances and furthermore that averaging radiances over large areas will reduce the total variation. A better understanding of the causes of temporal variation would improve the sensitivity of DNB to lighting changes.
Fashion design plays a significant role in Berlin's creative industries and for its start-up scene. Berlin has the highest concentration of designers in Germany, most of them working in small start-ups, while the spatial organisation of their production is stretched from the local level to the global network of fashion events, showing different entrepreneurial strategies within the production process. Different spatial structures of the production organisation are identifiable through which it is possible to discuss the role of Berlin in the production network of fashion designers and the kinds of relations holding between the city, designers, and their production network.
Continental rift systems form by propagation of isolated rift segments that interact, and eventually evolve into continuous zones of deformation. This process impacts many aspects of rifting including rift morphology at breakup, and eventual ocean-ridge segmentation. Yet, rift segment growth and interaction remain enigmatic. Here we present geological data from the poorly documented Ririba rift (South Ethiopia) that reveals how two major sectors of the East African rift, the Kenyan and Ethiopian rifts, interact. We show that the Ririba rift formed from the southward propagation of the Ethiopian rift during the Pliocene but this propagation was short-lived and aborted close to the Pliocene-Pleistocene boundary. Seismicity data support the abandonment of laterally offset, overlapping tips of the Ethiopian and Kenyan rifts. Integration with new numerical models indicates that rift abandonment resulted from progressive focusing of the tectonic and magmatic activity into an oblique, throughgoing rift zone of near pure extension directly connecting the rift sectors.
Channel transmission losses in drylands take place normally in extensive alluvial channels or streambeds underlain by fractured rocks. They can play an important role in streamflow rates, groundwater recharge, freshwater supply and channel-associated ecosystems. We aim to develop a process-oriented, semi-distributed channel transmission losses model, using process formulations which are suitable for data-scarce dryland environments and applicable to both hydraulically disconnected losing streams and hydraulically connected losing(/gaining) streams. This approach should be able to cover a large variation in climate and hydro-geologic controls, which are typically found in dryland regions of the Earth. Our model was first evaluated for a losing/gaining, hydraulically connected 30 km reach of the Middle Jaguaribe River (MJR), Ceara, Brazil, which drains a catchment area of 20 000 km(2). Secondly, we applied it to a small losing, hydraulically disconnected 1.5 km channel reach in the Walnut Gulch Experimental Watershed (WGEW), Arizona, USA. The model was able to predict reliably the streamflow volume and peak for both case studies without using any parameter calibration procedure. We have shown that the evaluation of the hypotheses on the dominant hydrological processes was fundamental for reducing structural model uncertainties and improving the streamflow prediction. For instance, in the case of the large river reach (MJR), it was shown that both lateral stream-aquifer water fluxes and groundwater flow in the underlying alluvium parallel to the river course are necessary to predict streamflow volume and channel transmission losses, the former process being more relevant than the latter. Regarding model uncertainty, it was shown that the approaches, which were applied for the unsaturated zone processes (highly nonlinear with elaborate numerical solutions), are much more sensitive to parameter variability than those approaches which were used for the saturated zone (mathematically simple water budgeting in aquifer columns, including backwater effects). In case of the MJR-application, we have seen that structural uncertainties due to the limited knowledge of the subsurface saturated system interactions (i.e. groundwater coupling with channel water; possible groundwater flow parallel to the river) were more relevant than those related to the subsurface parameter variability. In case of the WEGW application we have seen that the non-linearity involved in the unsaturated flow processes in disconnected dryland river systems (controlled by the unsaturated zone) generally contain far more model uncertainties than do connected systems controlled by the saturated flow. Therefore, the degree of aridity of a dryland river may be an indicator of potential model uncertainty and subsequent attainable predictability of the system.
In precipitation nowcasting, it is common to track the motion of precipitation in a sequence of weather radar images and to extrapolate this motion into the future. The total error of such a prediction consists of an error in the predicted location of a precipitation feature and an error in the change of precipitation intensity over lead time. So far, verification measures did not allow isolating the extent of location errors, making it difficult to specifically improve nowcast models with regard to location prediction. In this paper, we introduce a framework to directly quantify the location error. To that end, we detect and track scale-invariant precipitation features (corners) in radar images. We then consider these observed tracks as the true reference in order to evaluate the performance (or, inversely, the error) of any model that aims to predict the future location of a precipitation feature. Hence, the location error of a forecast at any lead time Delta t ahead of the forecast time t corresponds to the Euclidean distance between the observed and the predicted feature locations at t + Delta t. Based on this framework, we carried out a benchmarking case study using one year worth of weather radar composites of the German Weather Service. We evaluated the performance of four extrapolation models, two of which are based on the linear extrapolation of corner motion from t - 1 to t (LK-Lin1) and t - 4 to t (LK-Lin4) and the other two are based on the Dense Inverse Search (DIS) method: motion vectors obtained from DIS are used to predict feature locations by linear (DIS-Lin1) and Semi-Lagrangian extrapolation (DIS-Rot1). Of those four models, DIS-Lin1 and LK-Lin4 turned out to be the most skillful with regard to the prediction of feature location, while we also found that the model skill dramatically depends on the sinuosity of the observed tracks. The dataset of 376,125 detected feature tracks in 2016 is openly available to foster the improvement of location prediction in extrapolation-based nowcasting models.
Accurate weather observations are the keystone to many quantitative applications, such as precipitation monitoring and nowcasting, hydrological modelling and forecasting, climate studies, as well as understanding precipitation-driven natural hazards (i.e. floods, landslides, debris flow). Weather radars have been an increasingly popular tool since the 1940s to provide high spatial and temporal resolution precipitation data at the mesoscale, bridging the gap between synoptic and point scale observations. Yet, many institutions still struggle to tap the potential of the large archives of reflectivity, as there is still much to understand about factors that contribute to measurement errors, one of which is calibration. Calibration represents a substantial source of uncertainty in quantitative precipitation estimation (QPE). A miscalibration of a few dBZ can easily deteriorate the accuracy of precipitation estimates by an order of magnitude. Instances where rain cells carrying torrential rains are misidentified by the radar as moderate rain could mean the difference between a timely warning and a devastating flood.
Since 2012, the Philippine Atmospheric, Geophysical, and Astronomical Services Administration (PAGASA) has been expanding the country’s ground radar network. We had a first look into the dataset from one of the longest running radars (the Subic radar) after devastating week-long torrential rains and thunderstorms in August 2012 caused by the annual southwestmonsoon and enhanced by the north-passing Typhoon Haikui. The analysis of the rainfall spatial distribution revealed the added value of radar-based QPE in comparison to interpolated rain gauge observations. However, when compared with local gauge measurements, severe miscalibration of the Subic radar was found. As a consequence, the radar-based QPE would have underestimated the rainfall amount by up to 60% if they had not been adjusted by rain gauge observations—a technique that is not only affected by other uncertainties, but which is also not feasible in other regions of the country with very sparse rain gauge coverage.
Relative calibration techniques, or the assessment of bias from the reflectivity of two radars, has been steadily gaining popularity. Previous studies have demonstrated that reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and its successor, the Global Precipitation Measurement (GPM), are accurate enough to serve as a calibration reference for ground radars over low-to-mid-latitudes (± 35 deg for TRMM; ± 65 deg for GPM). Comparing spaceborne radars (SR) and ground radars (GR) requires cautious consideration of differences in measurement geometry and instrument specifications, as well as temporal coincidence. For this purpose, we implement a 3-D volume matching method developed by Schwaller and Morris (2011) and extended by Warren et al. (2018) to 5 years worth of observations from the Subic radar. In this method, only the volumetric intersections of the SR and GR beams are considered.
Calibration bias affects reflectivity observations homogeneously across the entire radar domain. Yet, other sources of systematic measurement errors are highly heterogeneous in space, and can either enhance or balance the bias introduced by miscalibration. In order to account for such heterogeneous errors, and thus isolate the calibration bias, we assign a quality index to each matching SR–GR volume, and thus compute the GR calibration bias as a qualityweighted average of reflectivity differences in any sample of matching SR–GR volumes. We exemplify the idea of quality-weighted averaging by using beam blockage fraction (BBF) as a quality variable. Quality-weighted averaging is able to increase the consistency of SR and GR observations by decreasing the standard deviation of the SR–GR differences, and thus increasing the precision of the bias estimates.
To extend this framework further, the SR–GR quality-weighted bias estimation is applied to the neighboring Tagaytay radar, but this time focusing on path-integrated attenuation (PIA) as the source of uncertainty. Tagaytay is a C-band radar operating at a lower wavelength and is therefore more affected by attenuation. Applying the same method used for the Subic radar, a time series of calibration bias is also established for the Tagaytay radar.
Tagaytay radar sits at a higher altitude than the Subic radar and is surrounded by a gentler terrain, so beam blockage is negligible, especially in the overlapping region. Conversely, Subic radar is largely affected by beam blockage in the overlapping region, but being an SBand radar, attenuation is considered negligible. This coincidentally independent uncertainty contributions of each radar in the region of overlap provides an ideal environment to experiment with the different scenarios of quality filtering when comparing reflectivities from the two ground radars. The standard deviation of the GR–GR differences already decreases if we consider either BBF or PIA to compute the quality index and thus the weights. However, combining them multiplicatively resulted in the largest decrease in standard deviation, suggesting that taking both factors into account increases the consistency between the matched samples.
The overlap between the two radars and the instances of the SR passing over the two radars at the same time allows for verification of the SR–GR quality-weighted bias estimation method. In this regard, the consistency between the two ground radars is analyzed before and after bias correction is applied. For cases when all three radars are coincident during a significant rainfall event, the correction of GR reflectivities with calibration bias estimates from SR overpasses dramatically improves the consistency between the two ground radars which have shown incoherent observations before correction. We also show that for cases where adequate SR coverage is unavailable, interpolating the calibration biases using a moving average can be used to correct the GR observations for any point in time to some extent. By using the interpolated biases to correct GR observations, we demonstrate that bias correction reduces the absolute value of the mean difference in most cases, and therefore improves the consistency between the two ground radars.
This thesis demonstrates that in general, taking into account systematic sources of uncertainty that are heterogeneous in space (e.g. BBF) and time (e.g. PIA) allows for a more consistent estimation of calibration bias, a homogeneous quantity. The bias still exhibits an unexpected variability in time, which hints that there are still other sources of errors that remain unexplored. Nevertheless, the increase in consistency between SR and GR as well as between the two ground radars, suggests that considering BBF and PIA in a weighted-averaging approach is a step in the right direction.
Despite the ample room for improvement, the approach that combines volume matching between radars (either SR–GR or GR–GR) and quality-weighted comparison is readily available for application or further scrutiny. As a step towards reproducibility and transparency in atmospheric science, the 3D matching procedure and the analysis workflows as well as sample data are made available in public repositories. Open-source software such as Python and wradlib are used for all radar data processing in this thesis. This approach towards open science provides both research institutions and weather services with a valuable tool that can be applied to radar calibration, from monitoring to a posteriori correction of archived data.
On April 29, 2017 at 0:56 UTC (2:56 local time), an M (W) = 2.8 earthquake struck the metropolitan area between Leipzig and Halle, Germany, near the small town of Markranstadt. The earthquake was felt within 50 km from the epicenter and reached a local intensity of I (0) = IV. Already in 2015 and only 15 km northwest of the epicenter, a M (W) = 3.2 earthquake struck the area with a similar large felt radius and I (0) = IV. More than 1.1 million people live in the region, and the unusual occurrence of the two earthquakes led to public attention, because the tectonic activity is unclear and induced earthquakes have occurred in neighboring regions. Historical earthquakes south of Leipzig had estimated magnitudes up to M (W) ae 5 and coincide with NW-SE striking crustal basement faults. We use different seismological methods to analyze the two recent earthquakes and discuss them in the context of the known tectonic structures and historical seismicity. Novel stochastic full waveform simulation and inversion approaches are adapted for the application to weak, local earthquakes, to analyze mechanisms and ground motions and their relation to observed intensities. We find NW-SE striking normal faulting mechanisms for both earthquakes and centroid depths of 26 and 29 km. The earthquakes are located where faults with large vertical offsets of several hundred meters and Hercynian strike have developed since the Mesozoic. We use a stochastic full waveform simulation to explain the local peak ground velocities and calibrate the method to simulate intensities. Since the area is densely populated and has sensitive infrastructure, we simulate scenarios assuming that a 12-km long fault segment between the two recent earthquakes is ruptured and study the impact of rupture parameters on ground motions and expected damage.
How fast the Northern Hemisphere (NH) forest biome tracks strongly warming climates is largely unknown. Regional studies reveal lags between decades and millennia. Here we report a conundrum: Deglacial forest expansion in the NH extra-tropics occurs approximately 4000 years earlier in a transient MPI-ESM1.2 simulation than shown by pollen-based biome reconstructions. Shortcomings in the model and the reconstructions could both contribute to this mismatch, leaving the underlying causes unresolved. The simulated vegetation responds within decades to simulated climate changes, which agree with pollen-independent reconstructions. Thus, we can exclude climate biases as main driver for differences. Instead, the mismatch points at a multi-millennial disequilibrium of the NH forest biome to the climate signal. Therefore, the evaluation of time-slice simulations in strongly changing climates with pollen records should be critically reassessed. Our results imply that NH forests may be responding much slower to ongoing climate changes than Earth System Models predict. <br /> Deglacial forest expansion in the Northern Hemisphere poses a conundrum: Model results agree with the climate signal but are several millennia ahead of reconstructed forest dynamics. The underlying causes remain unsolved.
A methodology is presented to assess the impact of reservoir silting oil water availability for semiarid environments, applied to seven representative watersheds in the state of Ceara, Brazil. Water yield is computed using stochastic modelling for several reliability levels and water yield reduction is quantified for the focus areas. The yield-volume elasticity concept, which indicates the relative yield reduction in terms of relative storage capacity of the reservoirs, is presented and applied. Results chow that storage capacity was reduced by 0.2% year(-1) due to silting, that the risk of water shortage almost doubled in less than 50 years for the most critical reservoir, and that reduction of storage capacity had three times more impact oil yield reduction than the increase in evaporation. Average 90% reliable yield-volume elasticity was 0.8, which means that the global water yield (Q(90)) in Ceara is expected to diminish yearly by 388 L s(-1) due to reservoir silting
Contemporary drought impact assessments have been constrained due to data availability, leading to an incomplete representation of impact trends. To address this, we present a novel method for the comprehensive and near-real-time monitoring of drought socio-economic impacts based on media reports. We tested its application using the case of the exceptional 2018/19 German drought. By employing text mining techniques, 4839 impact statements were identified, relating to livestock, agriculture, forestry, fires, recreation, energy and transport sectors. An accuracy of 95.6% was obtained for their automatic classification. Furthermore, high levels of performance in terms of spatial and temporal precision were found when validating our results against independent data (e.g. soil moisture, average precipitation, population interest in droughts, crop yield and forest fire statistics). The findings highlight the applicability of media data for rapidly and accurately monitoring the propagation of drought consequences over time and space. We anticipate our method to be used as a starting point for an impact-based early warning system.
1. Developing lesson plans and choosing strategies 2. The aims of the lesson plans in general 3. Strategies as a means to achieve theaims of the lesson plans 4. Evaluating the quality of lesson plans 5. Difficulties during lessons and adaptations afterwards 6. Student teachers’ overall feeling about their work 7. Using the strategies in future classes 8. Conclusion
Formation of epithermal Sn-Ag-(Zn) vein-type mineralization at the Pirquitas deposit, NW Argentina
(2018)
The Pirquitas Sn-Ag-(Zn) deposit in northwestern Argentina is thought to be an analogue to the Miocene polymetallic epithermal Sn-Ag deposits of the southern Bolivian Tin Belt, but little is known in detail about the origin and evolution of ore-forming fluids at Pirquitas. This paper reports on a microthermometric study of fluid inclusions in quartz, sphalerite, Ag-Sn sulfides, and Ag-rich sulfosalts using transmitted near infrared and visible light, combined with noble gas isotope analyses of fluids released from mineral separates. The study focused on the vein-hosted mineralization, which formed during two major mineralization events, whereby the first event I comprises two stages (I-1 and I-2). All studied minerals exclusively contain aqueous two-phase inclusions, indicating that the ore-forming fluids did not undergo two-phase phase separation (boiling). Salinity of fluid inclusions in I-1 quartz that precipitated along with pyrite and pyrrhotite ranges between 0 and 7.5 wt% NaCl equiv. and homogenization temperatures (Th) are between 233 and 370 degrees C. Stage I-2 is characterized by abundant Sn-Ag-Pb-Zn-sulfides and a variety of Ag-rich sulfosalts. Fluid inclusions in stage I-2 Ag-Sn sulfides have salinities up to 10.6 wt% NaCl equiv. and Th between 213 and 274 degrees C. The deposition of stage I-2 ore is likely related to a new pulse of saline magmatic fluids to the hydrothermal system. The mineralization event II deposited the richest Ag ores at Pirquitas. Colloform sphalerite and pyrargyrite deposited during event II contain two-phase aqueous fluid inclusions with homogenization temperatures between 190 and 252 degrees C and salinities between 0.9 and 4.3 wt% NaCl equiv. Noble gas concentrations and isotopic compositions of ore-hosted fluid inclusions were determined from crushing hand-picked ore minerals from both mineralization events. With one exception, all samples yielded He-3/He-4 ratios between 1.9 and 4.1 Ra, which is within the range of published data from the volcanic arc and somewhat higher than typical values of meteoric water-derived hot-springs in the region. This demonstrates a significant contribution of magmatic fluids to the Pirquitas mineralization although no intrusive rocks are exposed in the mine region. Taking the noble gas evidence for a magmatic fluid source, we interpret the trends of decreasing Th and salinity values in fluid inclusions from events I and II to represent waning of the magmatic-hydrothermal system and/or increased admixing of meteoric water to the magmatic fluids.
Trends in streamflow, rainfall and potential evapotranspiration (PET) time series, from 1970 to 2017, were assessed for five important hydrological basins in Southeastern Brazil. The concept of elasticity was also used to assess the streamflow sensitivity to changes in climate variables, for annual data and 5-, 10- and 20-year moving averages. Significant negative trends in streamflow and rainfall and significant increasing trend in PET were detected. For annual analysis, elasticity revealed that 1% decrease in rainfall resulted in 1.21-2.19% decrease in streamflow, while 1% increase in PET induced different reductions percentages in streamflow, ranging from 2.45% to 9.67%. When both PET and rainfall were computed to calculate the elasticity, results were positive for some basins. Elasticity analysis considering 20-year moving averages revealed that impacts on the streamflow were cumulative: 1% decrease in rainfall resulted in 1.83-4.75% decrease in streamflow, while 1% increase in PET induced 3.47-28.3% decrease in streamflow. This different temporal response may be associated with the hydrological memory of the basins. Streamflow appears to be more sensitive in less rainy basins. This study provides useful information to support strategic government decisions, especially when the security of water resources and drought mitigation are considered in face of climate change.
The growing worldwide impact of flood events has motivated the development and application of global flood hazard models (GFHMs). These models have become useful tools for flood risk assessment and management, especially in regions where little local hazard information is available. One of the key uncertainties associated with GFHMs is the estimation of extreme flood magnitudes to generate flood hazard maps. In this study, the 1-in-100 year flood (Q100) magnitude was estimated using flow outputs from four global hydrological models (GHMs) and two global flood frequency analysis datasets for 1350 gauges across the conterminous US. The annual maximum flows of the observed and modelled timeseries of streamflow were bootstrapped to evaluate the sensitivity of the underlying data to extrapolation. Results show that there are clear spatial patterns of bias associated with each method. GHMs show a general tendency to overpredict Western US gauges and underpredict Eastern US gauges. The GloFAS and HYPE models underpredict Q100 by more than 25% in 68% and 52% of gauges, respectively. The PCR-GLOBWB and CaMa-Flood models overestimate Q100 by more than 25% at 60% and 65% of gauges in West and Central US, respectively. The global frequency analysis datasets have spatial variabilities that differ from the GHMs. We found that river basin area and topographic elevation explain some of the spatial variability in predictive performance found in this study. However, there is no single model or method that performs best everywhere, and therefore we recommend a weighted ensemble of predictions of extreme flood magnitudes should be used for large-scale flood hazard assessment.
Study region: Tisza and Prut catchments, originating on the slopes of the Carpathian mountains. Study focus: The study reported here investigates (i) climate change impacts on flood risk in the region, and (ii) uncertainty related to hydrological modelling, downscaling techniques and climate projections. The climate projections used in the study were derived from five GCMs, downscaled either dynamically with RCMs or with the statistical downscaling model XDS. The resulting climate change scenarios were applied to drive the eco-hydrological model SWIM, which was calibrated and validated for the catchments in advance using observed climate and hydrological data. The changes in the 30-year flood hazards and 98 and 95 percentiles of discharge were evaluated for the far future period (2071-2100) in comparison with the reference period (1981-2010). New hydrological insights for the region: The majority of model outputs under RCP 4.5 show a small to strong increase of the 30-year flood level in the Tisza ranging from 4.5% to 62%, and moderate increase in the Prut ranging from 11% to 22%. The impact results under RCP 8.5 are more uncertain with changes in both directions due to high uncertainties in GCM-RCM climate projections, downscaling methods and the low density of available climate stations.
Landscapes in high northern latitudes are assumed to be highly sensitive to future global change, but the rates and long-term trajectories of changes are rather uncertain. In the boreal zone, fires are an important factor in climate-vegetation interactions and biogeochemical cycles. Fire regimes are characterized by small, frequent, low-intensity fires within summergreen boreal forests dominated by larch, whereas evergreen boreal forests dominated by spruce and pine burn large areas less frequently but at higher intensities. Here, we explore the potential of the monosaccharide anhydrides (MA) levoglucosan, mannosan and galactosan to serve as proxies of low-intensity biomass burning in glacial-to-interglacial lake sediments from the high northern latitudes. We use sediments from Lake El'gygytgyn (cores PG 1351 and ICDP 5011-1), located in the far north-east of Russia, and study glacial and interglacial samples of the last 430 kyr (marine isotope stages 5e, 6, 7e, 8, 11c and 12) that had different climate and biome configurations. Combined with pollen and non-pollen palynomorph records from the same samples, we assess how far the modern relationships between fire, climate and vegetation persisted during the past, on orbital to centennial timescales. We find that MAs attached to particulates were well-preserved in up to 430 kyr old sediments with higher influxes from low-intensity biomass burning in interglacials compared to glacials. MA influxes significantly increase when summergreen boreal forest spreads closer to the lake, whereas they decrease when tundra-steppe environments and, especially, Sphagnum peatlands spread. This suggests that low-temperature fires are a typical characteristic of Siberian larch forests also on long timescales. The results also suggest that low-intensity fires would be reduced by vegetation shifts towards very dry environments due to reduced biomass availability, as well as by shifts towards peatlands, which limits fuel dryness. In addition, we observed very low MA ratios, which we interpret as high contributions of galactosan and mannosan from biomass sources other than those currently monitored, such as the moss-lichen mats in the understorey of the summergreen boreal forest. Overall, sedimentary MAs can provide a powerful proxy for fire regime reconstructions and extend our knowledge of long-term natural fire-climate-vegetation feedbacks in the high northern latitudes.
Landscapes in high northern latitudes are assumed to be highly sensitive to future global change, but the rates and long-term trajectories of changes are rather uncertain. In the boreal zone, fires are an important factor in climate-vegetation interactions and biogeochemical cycles. Fire regimes are characterized by small, frequent, low-intensity fires within summergreen boreal forests dominated by larch, whereas evergreen boreal forests dominated by spruce and pine burn large areas less frequently but at higher intensities. Here, we explore the potential of the monosaccharide anhydrides (MA) levoglucosan, mannosan and galactosan to serve as proxies of low-intensity biomass burning in glacial-to-interglacial lake sediments from the high northern latitudes. We use sediments from Lake El'gygytgyn (cores PG 1351 and ICDP 5011-1), located in the far north-east of Russia, and study glacial and interglacial samples of the last 430 kyr (marine isotope stages 5e, 6, 7e, 8, 11c and 12) that had different climate and biome configurations. Combined with pollen and non-pollen palynomorph records from the same samples, we assess how far the modern relationships between fire, climate and vegetation persisted during the past, on orbital to centennial timescales. We find that MAs attached to particulates were well-preserved in up to 430 kyr old sediments with higher influxes from low-intensity biomass burning in interglacials compared to glacials. MA influxes significantly increase when summergreen boreal forest spreads closer to the lake, whereas they decrease when tundra-steppe environments and, especially, Sphagnum peatlands spread. This suggests that low-temperature fires are a typical characteristic of Siberian larch forests also on long timescales. The results also suggest that low-intensity fires would be reduced by vegetation shifts towards very dry environments due to reduced biomass availability, as well as by shifts towards peatlands, which limits fuel dryness. In addition, we observed very low MA ratios, which we interpret as high contributions of galactosan and mannosan from biomass sources other than those currently monitored, such as the moss-lichen mats in the understorey of the summergreen boreal forest. Overall, sedimentary MAs can provide a powerful proxy for fire regime reconstructions and extend our knowledge of long-term natural fire-climate-vegetation feedbacks in the high northern latitudes.
Climate change heavily impacts smallholder farming worldwide. Cross-scale vulnerability assessment has a high potential to identify nested measures for reducing vulnerability of smallholder farmers. Despite their high practical value, there are currently only limited examples of cross-scale assessments. The presented study aims at assessing the vulnerability of smallholder farmers in the Northeast of Brazil across three scales: regional, farm and field scale. In doing so, it builds on existing vulnerability indices and compares results between indices at the same scale and across scales. In total, six independent indices are tested, two at each scale. The calculated indices include social, economic and ecological indicators, based on municipal statistics, meteorological data, farm interviews and soil analyses. Subsequently, indices and overlapping indicators are normalized for intra- and cross-scale comparison. The results show considerable differences between indices across and within scales. They indicate different activities to reduce vulnerability of smallholder farmers. Major shortcomings arise from the conceptual differences between the indices. We therefore recommend the development of hierarchical indices, which are adapted to local conditions and contain more overlapping indicators for a better understanding of the nested vulnerabilities of smallholder farmers.
Climate change heavily impacts smallholder farming worldwide. Cross-scale vulnerability assessment has a high potential to identify nested measures for reducing vulnerability of smallholder farmers. Despite their high practical value, there are currently only limited examples of cross-scale assessments. The presented study aims at assessing the vulnerability of smallholder farmers in the Northeast of Brazil across three scales: regional, farm and field scale. In doing so, it builds on existing vulnerability indices and compares results between indices at the same scale and across scales. In total, six independent indices are tested, two at each scale. The calculated indices include social, economic and ecological indicators, based on municipal statistics, meteorological data, farm interviews and soil analyses. Subsequently, indices and overlapping indicators are normalized for intra- and cross-scale comparison. The results show considerable differences between indices across and within scales. They indicate different activities to reduce vulnerability of smallholder farmers. Major shortcomings arise from the conceptual differences between the indices. We therefore recommend the development of hierarchical indices, which are adapted to local conditions and contain more overlapping indicators for a better understanding of the nested vulnerabilities of smallholder farmers.
The objectives of the present investigation are (i) to study the effects of climate change on precipitation extremes and (ii) to assess the uncertainty in the climate projections. The investigation is performed on the Lech catchment, located in the Northern Limestone Alps. In order to estimate the uncertainty in the climate projections, two statistical downscaling models as well as a number of global and regional climate models were considered. The downscaling models applied are the Expanded Downscaling (XDS) technique and the Long Ashton Research Station Weather Generator (LARS-WG). The XDS model, which is driven by analyzed or simulated large-scale synoptic fields, has been calibrated using ECMWF-interim reanalysis data and local station data. LARS-WG is controlled through stochastic parameters representing local precipitation variability, which are calibrated from station data only. Changes in precipitation mean and variability as simulated by climate models were then used to perturb the parameters of LARS-WG in order to generate climate change scenarios. In our study we use climate simulations based on the A1B emission scenario. The results show that both downscaling models perform well in reproducing observed precipitation extremes. In general, the results demonstrate that the projections are highly variable. The choice of both the GCM and the downscaling method are found to be essential sources of uncertainty. For spring and autumn, a slight tendency toward an increase in the intensity of future precipitation extremes is obtained, as a number of simulations show statistically significant increases in the intensity of 90th and 99th percentiles of precipitation on wet days as well as the 5- and 20-yr return values.
The tropical peat swamp forests of South-East Asia are being rapidly converted to agricultural plantations of oil palm and Acacia creating a significant global “hot-spot” for CO2 emissions. However, the effect of this major perturbation has yet to be quantified in terms of global warming potential (GWP) and the Earth's radiative budget. We used a GWP analysis and an impulse-response model of radiative forcing to quantify the climate forcing of this shift from a long-term carbon sink to a net source of greenhouse gases (CO2 and CH4). In the GWP analysis, five tropical peatlands were sinks in terms of their CO2 equivalent fluxes while they remained undisturbed. However, their drainage and conversion to oil palm and Acacia plantations produced a dramatic shift to very strong net CO2-equivalent sources. The induced losses of peat carbon are ~20× greater than the natural CO2 sequestration rates. In contrast, a radiative forcing model indicates that the magnitude of this shift from a net cooling to warming effect is ultimately related to the size of an individual peatland's carbon pool. The continuous accumulation of carbon in pristine tropical peatlands produced a progressively negative radiative forcing (i.e., cooling) that ranged from −2.1 to −6.7 nW/m2 per hectare peatland by 2010 CE, referenced to zero at the time of peat initiation. Peatland conversion to plantations leads to an immediate shift from negative to positive trend in radiative forcing (i.e., warming). If drainage persists, peak warming ranges from +3.3 to +8.7 nW/m2 per hectare of drained peatland. More importantly, this net warming impact on the Earth's radiation budget will persist for centuries to millennia after all the peat has been oxidized to CO2. This previously unreported and undesirable impact on the Earth's radiative balance provides a scientific rationale for conserving tropical peatlands in their pristine state.
Natural earthquakes often have very few observable foreshocks which significantly complicates tracking potential preparatory processes. To better characterize expected preparatory processes before failures, we study stick-slip events in a series of triaxial compression tests on faulted Westerly granite samples. We focus on the influence of fault roughness on the duration and magnitude of recordable precursors before large stick-slip failure. Rupture preparation in the experiments is detectable over long time scales and involves acoustic emission (AE) and aseismic deformation events. Preparatory fault slip is found to be accelerating during the entire pre-failure loading period, and is accompanied by increasing AE rates punctuated by distinct activity spikes associated with large slip events. Damage evolution across the fault zones and surrounding wall rocks is manifested by precursory decrease of seismic b-values and spatial correlation dimensions. Peaks in spatial event correlation suggest that large slip initiation occurs by failure of multiple asperities. Shear strain estimated from AE data represents only a small fraction (< 1%) of total shear strain accumulated during the preparation phase, implying that most precursory deformation is aseismic. The relative contribution of aseismic deformation is amplified by larger fault roughness. Similarly, seismic coupling is larger for smooth saw-cut faults compared to rough faults. The laboratory observations point towards a long-lasting and continuous preparation process leading to failure and large seismic events. The strain partitioning between aseismic and observable seismic signatures depends on fault structure and instrument resolution.
Rock glaciers are permafrost or glacial landforms of debris and ice that deform under the influence of gravity. Recent estimates hold that, in the semiarid Chilean Andes for example, active rock glaciers store more water than glaciers. However, little is known about how many rock glaciers might decay because of global warming and how much this decay might contribute to water and sediment release. We investigated an inventory of >6500 rock glaciers in the Argentinian Andes, spanning the climatic gradient from the Desert Andes to cold-temperate Tierra del Fuego. We used active rock glaciers as a diagnostic of permafrost, assuming that the toes mark the 0 degrees C isotherm in climate scenarios for the twenty-first century and their impact on freezing conditions near the rock glacier toes. We find that, under future worst case warming, up to 95% of rock glaciers in the southern Desert Andes and in the Central Andes will rest in areas above 0 degrees C and that this freezing level might move up more than twice as much (similar to 500 m) as during the entire Holocene (similar to 200 m). Many active rock glaciers are already well below the current freezing level and exemplify how local controls may confound regional prognoses. A Bayesian Multifactor Analysis of Variance further shows that only in the Central Andes are the toes of active rock glaciers credibly higher than those of inactive ones. Elsewhere in the Andes, active and inactive rock glaciers occupy indistinguishable elevation bands, regardless of aspect, the formation mechanism, or shape of rock glaciers. The state of rock glacier activity predicts differences in elevations of toes to 140 m at best so that regional inference of the distribution of discontinuous permafrost from rock-glacier toes cannot be more accurate than this in the Argentinian Andes. We conclude that the Central Andes-where rock glaciers are largest, cover the most area, and have a greater density than glaciers-is likely to experience the most widespread disturbance to the thermal regime of the twenty-first century. (C) 2018 Elsevier B.V. All rights reserved.
Understanding hydrological processes is of fundamental importance for the Vietnamese national food security and the livelihood of the population in the Vietnamese Mekong Delta (VMD). As a consequence of sparse data in this region, however, hydrologic processes, such as the controlling processes of precipitation, the interaction between surface and groundwater, and groundwater dynamics, have not been thoroughly studied. The lack of this knowledge may negatively impact the long-term strategic planning for sustainable groundwater resources management and may result in insufficient groundwater recharge and freshwater scarcity. It is essential to develop useful methods for a better understanding of hydrological processes in such data-sparse regions. The goal of this dissertation is to advance methodologies that can improve the understanding of fundamental hydrological processes in the VMD, based on the analyses of stable water isotopes and monitoring data. The thesis mainly focuses on the controlling processes of precipitation, the mechanism of surface–groundwater interaction, and the groundwater dynamics. These processes have not been fully addressed in the VMD so far. The thesis is based on statistical analyses of the isotopic data of Global Network of Isotopes in Precipitation (GNIP), of meteorological and hydrological data from Vietnamese agencies, and of the stable water isotopes and monitoring data collected as part of this work.
First, the controlling processes of precipitation were quantified by the combination of trajectory analysis, multi-factor linear regression, and relative importance analysis (hereafter, a model‐based statistical approach). The validity of this approach is confirmed by similar, but mainly qualitative results obtained in other studies. The total variation in precipitation isotopes (δ18O and δ2H) can be better explained by multiple linear regression (up to 80%) than single-factor linear regression (30%). The relative importance analysis indicates that atmospheric moisture regimes control precipitation isotopes rather than local climatic conditions. The most crucial factor is the upstream rainfall along the trajectories of air mass movement. However, the influences of regional and local climatic factors vary in importance over the seasons. The developed model‐based statistical approach is a robust tool for the interpretation of precipitation isotopes and could also be applied to understand the controlling processes of precipitation in other regions.
Second, the concept of the two-component lumped-parameter model (LPM) in conjunction with stable water isotopes was applied to examine the surface–groundwater interaction in the VMD. A calibration framework was also set up to evaluate the behaviour, parameter identifiability, and uncertainties of two-component LPMs. The modelling results provided insights on the subsurface flow conditions, the recharge contributions, and the spatial variation of groundwater transit time. The subsurface flow conditions at the study site can be best represented by the linear-piston flow distribution. The contributions of the recharge sources change with distance to the river. The mean transit time (mTT) of riverbank infiltration increases with the length of the horizontal flow path and the decreasing gradient between river and groundwater. River water infiltrates horizontally mainly via the highly permeable aquifer, resulting in short mTTs (<40 weeks) for locations close to the river (<200 m). The vertical infiltration from precipitation takes place primarily via a low‐permeable overlying aquitard, resulting in considerably longer mTTs (>80 weeks). Notably, the transit time of precipitation infiltration is independent of the distance to the river. All these results are hydrologically plausible and could be quantified by the presented method for the first time. This study indicates that the highly complex mechanism of surface–groundwater interaction at riverbank infiltration systems can be conceptualized by exploiting two‐component LPMs. It is illustrated that the model concept can be used as a tool to investigate the hydrological functioning of mixing processes and the flow path of multiple water components in riverbank infiltration systems.
Lastly, a suite of time series analysis approaches was applied to examine the groundwater dynamics in the VMD. The assessment was focused on the time-variant trends of groundwater levels (GWLs), the groundwater memory effect (representing the time that an aquifer holds water), and the hydraulic response between surface water and multi-layer alluvial aquifers. The analysis indicates that the aquifers act as low-pass filters to reduce the high‐frequency signals in the GWL variations, and limit the recharge to the deep groundwater. The groundwater abstraction has exceeded groundwater recharge between 1997 and 2017, leading to the decline of groundwater levels (0.01-0.55 m/year) in all considered aquifers in the VMD. The memory effect varies according to the geographical location, being shorter in shallow aquifers and flood-prone areas and longer in deep aquifers and coastal regions. Groundwater depth, season, and location primarily control the variation of the response time between the river and alluvial aquifers. These findings are important contributions to the hydrogeological literature of a little-known groundwater system in an alluvial setting. It is suggested that time series analysis can be used as an efficient tool to understand groundwater systems where resources are insufficient to develop a physical-based groundwater model.
This doctoral thesis demonstrates that important aspects of hydrological processes can be understood by statistical analysis of stable water isotope and monitoring data. The approaches developed in this thesis can be easily transferred to regions in similar tropical environments, particularly those in alluvial settings. The results of the thesis can be used as a baseline for future isotope-based studies and contribute to the hydrogeological literature of little-known groundwater systems in the VMD.
Water resources from Central Asia’s mountain regions have a high relevance for the water supply of the water scarce lowlands. A good understanding of the water cycle in these mountain regions is therefore needed to develop water management strategies. Hydrological modeling helps to improve our knowledge of the regional water cycle, and it can be used to gain a better understanding of past changes or estimate future hydrologic changes in view of projected changes in climate. However, due to the scarcity of hydrometeorological data, hydrological modeling for mountain regions in Central Asia involves large uncertainties.
Addressing this problem, the first aim of this thesis was to develop hydrological modeling approaches that can increase the credibility of hydrological models in data sparse mountain regions. This was achieved by using additional data from remote sensing and atmospheric modeling. It was investigated whether spatial patterns from downscaled reanalysis data can be used for the interpolation of station-based precipitation data. This approach was compared to other precipitation estimates using a hydrologic evaluation based on hydrological modeling and a comparison of simulated and observed discharge, which demonstrated a generally good performance of this method. The study further investigated the value of satellite-derived snow cover data for model calibration. Trade-offs of good model performance in terms of discharge and snow cover were explicitly evaluated using a multiobjective optimization algorithm, and the results were contrasted with single-objective calibration and Monte Carlo simulations. The study clearly shows that the additional use of snow cover data improved the internal consistency of the hydrological model. In this context, it was further investigated for the first time how many snow cover scenes were required for hydrological model calibration.
The second aim of this thesis was the application of the hydrological model in order to investigate the causes of observed streamflow increases in two headwater catchments of the Tarim River over the recent decades. This simulation-based approach for trend attribution was complemented by a data-based approach. The hydrological model was calibrated to discharge and glacier mass balance data and considered changes in glacier geometry over time. The results show that in the catchment with a lower glacierization, increasing precipitation and temperature both contributed to the streamflow increases, while in the catchment with a stronger glacierization, increasing temperatures were identified as the dominant driver.
Despite its many challenges and limitations the concept of in situ upgrading of informal settlements has become one of the most favoured approaches to the housing crisis in the ‘Global South’. Due to its inherent principles of incremental in situ development, prevention of relocations, protection of local livelihoods and democratic participation and cooperation, this approach is often perceived to be more sustainable than other housing approaches that often rely on quantitative housing delivery and top down planning methodologies. While this study does not question the benefits of the in situ upgrading approach, it seeks to identify problems of its practical implementation within a specific national and local context. The study discusses the origin and importance of this approach on the basis of a review of international housing policy development and analyses the broader political and social context of the incorporation of this approach into South African housing policy. It further uses insights from a recent case study in Cape Town to determine complications and conflicts that can arise when applying in situ upgrading of informal settlements in a complex local context. On that basis benefits and limitations of the in situ upgrading approach are specified and prerequisites for its successful implementation formulated.
The topography of first-order catchments in a region of western Amazonia was found to exhibit distinctive, recurrent features: a steep, straight lower side slope, a flat or nearly flat terrace at an intermediate elevation between valley floor and interfluve, and an upper side slope connecting interfluve and intermediate terrace. A detailed survey of soil-saturated hydraulic conductivity (K sat)-depth relationships, involving 740 undisturbed soil cores, was conducted in a 0.75-ha first-order catchment. The sampling approach was stratified with respect to the above slope units. Exploratory data analysis suggested fourth-root transformation of batches from the 0–0.1 m depth interval, log transformation of batches from the subsequent 0.1 m depth increments, and the use of robust estimators of location and scale. The K sat of the steep lower side slope decreased from 46 to 0.1 mm/h over the overall sampling depth of 0.4 m. The corresponding decrease was from 46 to 0.1 mm/h on the intermediate terrace, from 335 to 0.01 mm/h on the upper side slope, and from 550 to 0.015 mm/h on the interfluve. A depthwise comparison of these slope units led to the formulation of several hypotheses concerning the link between K sat and topography.
Rainfall erosivities as defined by the R factor from the universal soil loss equation were determined for all events during a two-year period at the station La Cuenca in western Amazonia. Three methods based on a power relationship between rainfall amount and erosivity were then applied to estimate event and daily rainfall erosivities from the respective rainfall amounts. A test of the resulting regression equations against an independent data set proved all three methods equally adequate in predicting rainfall erosivity from daily rainfall amount. We recommend the Richardson model for testing in the Amazon Basin, and its use with the coefficient from La Cuenca in western Amazonia.
Throughfall measurements were made under primary terra firme rainforest in the Rio Pichis valley, in the Upper Amazon Basin of Peru. Based on 214 precipitation events over nearly 18 months, throughfall was estimated to be 83.1±8.8% of gross precipitation. Regression analysis of all events revealed that gross precipitation is the only significant explanatory variable; the use of one-burst events does not significantly improve the regression relationship. Gross precipitation is, however, a poor predictor of throughfall for small rainfall events. The two forest structure parameters, canopy capacity, S, and free throughfall coefficient, p, were determined to be 1.3±0.2 mm and 0.32±0.18 mm. Rainfall intensity was found to influence these parameters. New methods which attempt to minimize the influence of meteorologic variables are used to estimate the potential values of these canopy parameters.
Chemical fingerprints of hydrological compartments and flow paths at La Cuenca, western Amazonia
(1995)
A forested first-order catchment in western Amazonia was monitored for 2 years to determine the chemical fingerprints of precipitation, throughfall, overland flow, pipe flow, soil water, groundwater, and streamflow. We used five tracers (hydrogen, calcium, magnesium, potassium, and silica) to distinguish “fast” flow paths mainly influenced by the biological subsystem from “slow” flow paths in the geochemical subsystem. The former comprise throughfall, overland flow, and pipe flow and are characterized by a high potassium/silica ratio; the latter are represented by soil water and groundwater, which have a low potassium/silica ratio. Soil water and groundwater differ with respect to calcium and magnesium. The groundwater-controlled streamflow chemistry is strongly modified by contributions from fast flow paths during precipitation events. The high potassium/silica ratio of these flow paths suggests that the storm flow response at La Cuenca is dominated by event water.
Previous hydrometric studies demonstrated the prevalence of overland flow as a hydrological pathway in the tropical rain forest catchment of South Creek, northeast Queensland. The purpose of this study was to consider this information in a mixing analysis with the aim of identifying sources of, and of estimating their contribution to, storm flow during two events in February 1993. K and acid-neutralizing capacity (ANC) were used as tracers because they provided the best separation of the potential sources, saturation overland flow, soil water from depths of 0.3, 0.6, and 1.2 m, and hillslope groundwater in a two-dimensional mixing plot. It was necessary to distinguish between saturation overland flow, generated at the soil surface and following unchanneled pathways, and overland flow in incised pathways. This latter type of overland flow was a mixture of saturation overland flow (event water) with high concentrations of K and a low ANC, soil water (preevent water) with low concentrations of K and a low ANC, and groundwater (preevent water) with low concentrations of K and a high ANC. The same sources explained the streamwater chemistry during the two events with strongly differing rainfall and antecedent moisture conditions. The contribution of saturation overland flow dominated the storm flow during the first, high-intensity, 178-mm event, while the contribution of soil water reached 50% during peak flow of the second, low-intensity, 44-mm event 5 days later. This latter result is remarkably similar to soil water contributions to storm flow in mountainous forested catchments of the southeastern United States. In terms of event and preevent water the storm flow hydrograph of the high-intensity event is dominated by event water and that of the low-intensity event by preevent water. This study highlights the problems of applying mixing analyses to overland flow-dominated catchments and soil environments with a poorly developed vertical chemical zonation and emphasizes the need for independent hydrometric information for a complete characterization of watershed hydrology and chemistry.
Spinning up large-scale coupled surface-subsurface numerical models can be a time and resource consuming task. If an uninformed initial condition is chosen, the spin-up can easily require 20 years of repeated simulations on high-performance computing machines. In this paper we compare the classical approach of starting from a fixed shallow depth to groundwater (here 3 m) with three more informed approaches for the definition of initial conditions in the spin up. In the first of these three approaches, we start from a known-steady state groundwater table, calculated with a 2-D groundwater model and the yearly net recharge, and combine it with an unsaturated zone that assumes hydrostatic conditions. In the second approach, we start from the same groundwater table combined with vertical profiles in the unsaturated zone with uniform vertical flow identical to the groundwater recharge. In the third approach we calculate a dynamic steady state from a simplified subsurface model combining a transient 2-D groundwater model with a limited number of 1-D transient unsaturated zone columns on top. Results for spinning-up a 3-D Parflow-CLM model using the different initial conditions show that large gains can be made by considering states in groundwater and the vadose zone that are consistent, i.e. where groundwater recharge and the vertical flux in the vadose zone agree. By this, the spin-up time was reduced from about 10 years to about 3 years of simulated time. In the light of seasonal fluctuations of net recharge, using the transient approach showed more stable results.
The Postmasburg Manganese Field (PMF), Northern Cape Province, South Africa, once represented one of the largest sources of manganese ore worldwide. Two belts of manganese ore deposits have been distinguished in the PMF, namely the Western Belt of ferruginous manganese ores and the Eastern Belt of siliceous manganese ores. Prevailing models of ore formation in these two belts invoke karstification of manganese-rich dolomites and residual accumulation of manganese wad which later underwent diagenetic and low-grade metamorphic processes. For the most part, the role of hydrothermal processes and metasomatic alteration towards ore formation has not been adequately discussed. Here we report an abundance of common and some rare Al-, Na-, K- and Ba-bearing minerals, particularly aegirine, albite, microcline, banalsite, serandite-pectolite, paragonite and natrolite in Mn ores of the PMF, indicative of hydrothermal influence. Enrichments in Na, K and/or Ba in the ores are generally on a percentage level for most samples analysed through bulk-rock techniques. The presence of As-rich tokyoite also suggests the presence of As and V in the hydrothermal fluid. The fluid was likely oxidized and alkaline in nature, akin to a mature basinal brine. Various replacement textures, particularly of Na- and K- rich minerals by Ba-bearing phases, suggest sequential deposition of gangue as well as ore-minerals from the hydrothermal fluid, with Ba phases being deposited at a later stage. The stratigraphic variability of the studied ores and their deviation from the strict classification of ferruginous and siliceous ores in the literature, suggests that a re-evaluation of genetic models is warranted. New Ar-Ar ages for K-feldspars suggest a late Neoproterozoic timing for hydrothermal activity. This corroborates previous geochronological evidence for regional hydrothermal activity that affected Mn ores at the PMF but also, possibly, the high-grade Mn ores of the Kalahari Manganese Field to the north. A revised, all-encompassing model for the development of the manganese deposits of the PMF is then proposed, whereby the source of metals is attributed to underlying carbonate rocks beyond the Reivilo Formation of the Campbellrand Subgroup. The main process by which metals are primarily accumulated is attributed to karstification of the dolomitic substrate. The overlying Asbestos Hills Subgroup banded iron formation (BIF) is suggested as a potential source of alkali metals, which also provides a mechanism for leaching of these BIFs to form high-grade residual iron ore deposits.
Large earthquakes initiate chains of surface processes that last much longer than the brief moments of strong shaking. Most moderate‐ and large‐magnitude earthquakes trigger landslides, ranging from small failures in the soil cover to massive, devastating rock avalanches. Some landslides dam rivers and impound lakes, which can collapse days to centuries later, and flood mountain valleys for hundreds of kilometers downstream. Landslide deposits on slopes can remobilize during heavy rainfall and evolve into debris flows. Cracks and fractures can form and widen on mountain crests and flanks, promoting increased frequency of landslides that lasts for decades. More gradual impacts involve the flushing of excess debris downstream by rivers, which can generate bank erosion and floodplain accretion as well as channel avulsions that affect flooding frequency, settlements, ecosystems, and infrastructure. Ultimately, earthquake sequences and their geomorphic consequences alter mountain landscapes over both human and geologic time scales. Two recent events have attracted intense research into earthquake‐induced landslides and their consequences: the magnitude M 7.6 Chi‐Chi, Taiwan earthquake of 1999, and the M 7.9 Wenchuan, China earthquake of 2008. Using data and insights from these and several other earthquakes, we analyze how such events initiate processes that change mountain landscapes, highlight research gaps, and suggest pathways toward a more complete understanding of the seismic effects on the Earth's surface.
In 2010, project CoCoCo (incipient COntinent-COntinent COllision) recorded a 650 km long amphibian N-S wide-angle seismic profile, extending from the Anatolian plateau across southern Turkey and Cyprus to just south of the Eratosthenes Seamount (ESM). The aim of the project is to reveal the impact of the transition from subduction to continent-continent collision of Africa with Anatolia. Arrival picking, finite-differences ray-tracing and inversion of the offshore and on-offshore data produced a tomographic model across southern Cyprus, the accretionary wedge and the ESM. The main features of this model are (1) crustal P-velocities predominantly lower than 6.5 km/s beneath the ESM, (2) crustal thickness between 28 and 37 km, (3) an upper crustal reflection at 5 km depth beneath the ESM, (4) the likely presence of oceanic crust south of the ESM and a transform margin north of it and (5) a 12 km thick ophiolite sequence on Cyprus. Land shots on Turkey, also recorded on Cyprus, gravity data and geological and previous seismic investigations allow to derive a layered velocity model beneath Anatolia and the northern part of Cyprus. The main features of this model are (1) Moho depths of 38–45 km beneath the Anatolian plateau, (2) an upper and lower crust with large lateral changes in velocity and thickness, (3) a north-dipping subducting plate below Cyprus with a steepening of the dip-angle of the plate at about 45 km depth. Thus, the wide-angle seismic and gravity data provide detailed insights into the 2-D geometry and velocity structures associated with the Cyprus Arc collision zone. Finally, integrated analysis of the geophysics and geology allows a comprehensive interpretation of the crustal structure related to the collision process.
Spatial population mobility as social interaction : A Fijian Island's multi-local village community
(1999)
High-mountain regions provide valuable ecosystem services, including food, water, and energy production, to more than 900 million people worldwide. Projections hold, that this population number will rapidly increase in the next decades, accompanied by a continued urbanisation of cities located in mountain valleys. One of the manifestations of this ongoing socio-economic change of mountain societies is a rise in settlement areas and transportation infrastructure while an increased power need fuels the construction of hydropower plants along rivers in the high-mountain regions of the world. However, physical processes governing the cryosphere of these regions are highly sensitive to changes in climate and a global warming will likely alter the conditions in the headwaters of high-mountain rivers. One of the potential implications of this change is an increase in frequency and magnitude of outburst floods – highly dynamic flows capable of carrying large amounts of water and sediments. Sudden outbursts from lakes formed behind natural dams are complex geomorphological processes and are often part of a hazard cascade. In contrast to other types of natural hazards in high-alpine areas, for example landslides or avalanches, outburst floods are highly infrequent. Therefore, observations and data describing for example the mode of outburst or the hydraulic properties of the downstream propagating flow are very limited, which is a major challenge in contemporary (glacial) lake outburst flood research. Although glacial lake outburst floods (GLOFs) and landslide-dammed lake outburst floods (LLOFs) are rare, a number of documented events caused high fatality counts and damage. The highest documented losses due to outburst floods since the start of the 20th century were induced by only a few high-discharge events. Thus, outburst floods can be a significant hazard to downvalley communities and infrastructure in high-mountain regions worldwide.
This thesis focuses on the Greater Himalayan region, a vast mountain belt stretching across 0.89 million km2. Although potentially hundreds of outburst floods have occurred there since the beginning of the 20th century, data on these events is still scarce. Projections of cryospheric change, including glacier-mass wastage and permafrost degradation, will likely result in an overall increase of the water volume stored in meltwater lakes as well as the destabilisation of mountain slopes in the Greater Himalayan region. Thus, the potential for outburst floods to affect the increasingly more densely populated valleys of this mountain belt is also likely to increase in the future. A prime example of one of these valleys is the Pokhara valley in Nepal, which is drained by the Seti Khola, a river crossing one of the steepest topographic gradients in the Himalayas. This valley is also home to Nepal’s second largest, rapidly growing city, Pokhara, which currently has a population of more than half a million people – some of which live in informal settlements within the floodplain of the Seti Khola. Although there is ample evidence for past outburst floods along this river in recent and historic times, these events have hardly been quantified.
The main motivation of my thesis is to address the data scarcity on past and potential future outburst floods in the Greater Himalayan region, both at a regional and at a local scale. For the former, I compiled an inventory of >3,000 moraine-dammed lakes, of which about 1% had a documented sudden failure in the past four decades. I used this data to test whether a number of predictors that have been widely applied in previous GLOF assessments are statistically relevant when estimating past GLOF susceptibility. For this, I set up four Bayesian multi-level logistic regression models, in which I explored the credibility of the predictors lake area, lake-area dynamics, lake elevation, parent-glacier-mass balance, and monsoonality. By using a hierarchical approach consisting of two levels, this probabilistic framework also allowed for spatial variability on GLOF susceptibility across the vast study area, which until now had not been considered in studies of this scale. The model results suggest that in the Nyainqentanglha and Eastern Himalayas – regions with strong negative glacier-mass balances – lakes have been more prone to release GLOFs than in regions with less negative or even stable glacier-mass balances. Similarly, larger lakes in larger catchments had, on average, a higher probability to have had a GLOF in the past four decades. Yet, monsoonality, lake elevation, and lake-area dynamics were more ambiguous. This challenges the credibility of a lake’s rapid growth in surface area as an indicator of a pending outburst; a metric that has been applied to regional GLOF assessments worldwide.
At a local scale, my thesis aims to overcome data scarcity concerning the flow characteristics of the catastrophic May 2012 flood along the Seti Khola, which caused 72 fatalities, as well as potentially much larger predecessors, which deposited >1 km³ of sediment in the Pokhara valley between the 12th and 14th century CE. To reconstruct peak discharges, flow depths, and flow velocities of the 2012 flood, I mapped the extents of flood sediments from RapidEye satellite imagery and used these as a proxy for inundation limits. To constrain the latter for the Mediaeval events, I utilised outcrops of slackwater deposits in the fills of tributary valleys. Using steady-state hydrodynamic modelling for a wide range of plausible scenarios, from meteorological (1,000 m³ s-1) to cataclysmic outburst floods (600,000 m³ s-1), I assessed the likely initial discharges of the recent and the Mediaeval floods based on the lowest mismatch between sedimentary evidence and simulated flood limits. One-dimensional HEC-RAS simulations suggest, that the 2012 flood most likely had a peak discharge of 3,700 m³ s-1 in the upper Seti Khola and attenuated to 500 m³ s-1 when arriving in Pokhara’s suburbs some 15 km downstream.
Simulations of flow in two-dimensions with orders of magnitude higher peak discharges in ANUGA show extensive backwater effects in the main tributary valleys. These backwater effects match the locations of slackwater deposits and, hence, attest for the flood character of Mediaeval sediment pulses. This thesis provides first quantitative proof for the hypothesis, that the latter were linked to earthquake-triggered outbursts of large former lakes in the headwaters of the Seti Khola – producing floods with peak discharges of >50,000 m³ s-1.
Building on this improved understanding of past floods along the Seti Khola, my thesis continues with an analysis of the impacts of potential future outburst floods on land cover, including built-up areas and infrastructure mapped from high-resolution satellite and OpenStreetMap data. HEC-RAS simulations of ten flood scenarios, with peak discharges ranging from 1,000 to 10,000 m³ s-1, show that the relative inundation hazard is highest in Pokhara’s north-western suburbs. There, the potential effects of hydraulic ponding upstream of narrow gorges might locally sustain higher flow depths. Yet, along this reach, informal settlements and gravel mining activities are close to the active channel. By tracing the construction dynamics in two of these potentially affected informal settlements on multi-temporal RapidEye, PlanetScope, and Google Earth imagery, I found that exposure increased locally between three- to twentyfold in just over a decade (2008 to 2021).
In conclusion, this thesis provides new quantitative insights into the past controls on the susceptibility of glacial lakes to sudden outburst at a regional scale and the flow dynamics of propagating flood waves released by past events at a local scale, which can aid future hazard assessments on transient scales in the Greater Himalayan region. My subsequent exploration of the impacts of potential future outburst floods to exposed infrastructure and (informal) settlements might provide valuable inputs to anticipatory assessments of multiple risks in the Pokhara valley.
Pokhara (ca. 850 m a.s.l.), Nepal's second-largest city, lies at the foot of the Higher Himalayas and has more than tripled its population in the past 3 decades. Construction materials are in high demand in rapidly expanding built-up areas, and several informal settlements cater to unregulated sand and gravel mining in the Pokhara Valley's main river, the Seti Khola. This river is fed by the Sabche glacier below Annapurna III (7555 m a.s.l.), some 35 km upstream of the city, and traverses one of the steepest topographic gradients in the Himalayas. In May 2012 a sudden flood caused >70 fatalities and intense damage along this river and rekindled concerns about flood risk management. We estimate the flow dynamics and inundation depths of flood scenarios using the hydrodynamic model HEC-RAS (Hydrologic Engineering Center’s River Analysis System). We simulate the potential impacts of peak discharges from 1000 to 10 000 m3 s−1 on land cover based on high-resolution Maxar satellite imagery and OpenStreetMap data (buildings and road network). We also trace the dynamics of two informal settlements near Kaseri and Yamdi with high potential flood impact from RapidEye, PlanetScope, and Google Earth imagery of the past 2 decades. Our hydrodynamic simulations highlight several sites of potential hydraulic ponding that would largely affect these informal settlements and sites of sand and gravel mining. These built-up areas grew between 3- and 20-fold, thus likely raising local flood exposure well beyond changes in flood hazard. Besides these drastic local changes, about 1 % of Pokhara's built-up urban area and essential rural road network is in the highest-hazard zones highlighted by our flood simulations. Our results stress the need to adapt early-warning strategies for locally differing hydrological and geomorphic conditions in this rapidly growing urban watershed.
Pokhara (ca. 850 m a.s.l.), Nepal's second-largest city, lies at the foot of the Higher Himalayas and has more than tripled its population in the past 3 decades. Construction materials are in high demand in rapidly expanding built-up areas, and several informal settlements cater to unregulated sand and gravel mining in the Pokhara Valley's main river, the Seti Khola. This river is fed by the Sabche glacier below Annapurna III (7555 m a.s.l.), some 35 km upstream of the city, and traverses one of the steepest topographic gradients in the Himalayas. In May 2012 a sudden flood caused >70 fatalities and intense damage along this river and rekindled concerns about flood risk management. We estimate the flow dynamics and inundation depths of flood scenarios using the hydrodynamic model HEC-RAS (Hydrologic Engineering Center’s River Analysis System). We simulate the potential impacts of peak discharges from 1000 to 10 000 m3 s−1 on land cover based on high-resolution Maxar satellite imagery and OpenStreetMap data (buildings and road network). We also trace the dynamics of two informal settlements near Kaseri and Yamdi with high potential flood impact from RapidEye, PlanetScope, and Google Earth imagery of the past 2 decades. Our hydrodynamic simulations highlight several sites of potential hydraulic ponding that would largely affect these informal settlements and sites of sand and gravel mining. These built-up areas grew between 3- and 20-fold, thus likely raising local flood exposure well beyond changes in flood hazard. Besides these drastic local changes, about 1 % of Pokhara's built-up urban area and essential rural road network is in the highest-hazard zones highlighted by our flood simulations. Our results stress the need to adapt early-warning strategies for locally differing hydrological and geomorphic conditions in this rapidly growing urban watershed.
Glacial lakes in the Hindu Kush–Karakoram–Himalayas–Nyainqentanglha (HKKHN) region have grown rapidly in number and area in past decades, and some dozens have drained in catastrophic glacial lake outburst floods (GLOFs). Estimating regional susceptibility of glacial lakes has largely relied on qualitative assessments by experts, thus motivating a more systematic and quantitative appraisal. Before the backdrop of current climate-change projections and the potential of elevation-dependent warming, an objective and regionally consistent assessment is urgently needed. We use an inventory of 3390 moraine-dammed lakes and their documented outburst history in the past four decades to test whether elevation, lake area and its rate of change, glacier-mass balance, and monsoonality are useful inputs to a probabilistic classification model. We implement these candidate predictors in four Bayesian multi-level logistic regression models to estimate the posterior susceptibility to GLOFs. We find that mostly larger lakes have been more prone to GLOFs in the past four decades regardless of the elevation band in which they occurred. We also find that including the regional average glacier-mass balance improves the model classification. In contrast, changes in lake area and monsoonality play ambiguous roles. Our study provides first quantitative evidence that GLOF susceptibility in the HKKHN scales with lake area, though less so with its dynamics. Our probabilistic prognoses offer improvement compared to a random classification based on average GLOF frequency. Yet they also reveal some major uncertainties that have remained largely unquantified previously and that challenge the applicability of single models. Ensembles of multiple models could be a viable alternative for more accurately classifying the susceptibility of moraine-dammed lakes to GLOFs.
Speleothems, secondary cave carbonates, are important tools for climate reconstruction, especially as they often can be very precisely dated with the UTh method. If the U-Th method fails, dating becomes difficult, and often results in abandonment of samples and study sites. Radiocarbon dating is the only other radiometric dating technique applicable to the last similar to 50 ka, but presents complexities related to temporal variability of the reservoir effect in speleothems. Thus, radiocarbon dating of speleothems is not straightforward, and there are currently no publicly available tools to define proper age-depth relationships with this method. Here, we present an improved version of a previously published radiocarbon based age-depth modelling approach (star, Lechleitner et al., 2016b), which is now made freely available. The software is easy to use and provides the possibility to obtain reliable age-depth relationships, without prior knowledge of reservoir effects and their variability. In addition, star is able to detect and handle growth stops and phases with different growth rates. We test star on artificially constructed data sets and illustrate steps to improve the model performance. Furthermore, we apply the new approach to published radiocarbon data of U-Th dated stalagmites. This offers the possibility to investigate the strengths and weaknesses of the new approach with respect to potentially significant long term trends in the radiocarbon reservoir effect, which might otherwise remain undetected. In summary, we have produced a valuable software, which easily enables to construct age-depth relationships on the basis of reservoir effect disturbed radiocarbon measurements.
During the earliest Triassic microbial mats flourished in the photic zones of marginal seas, generating widespread microbialites. It has been suggested that anoxic conditions in shallow marine environments, linked to the end-Permian mass extinction, limited mat-inhibiting metazoans allowing for this microbialite expansion. The presence of a diverse suite of proxies indicating oxygenated shallow sea-water conditions (metazoan fossils, biomarkers and redox proxies) from microbialite successions have, however, challenged the inference of anoxic conditions. Here, the distribution and faunal composition of Griesbachian microbialites from China, Iran, Turkey, Armenia, Slovenia and Hungary are investigated to determine the factors that allowed microbialite-forming microbial mats to flourish following the end-Permian crisis. The results presented here show that Neotethyan microbial buildups record a unique faunal association due to the presence of keratose sponges, while the Palaeotethyan buildups have a higher proportion of molluscs and the foraminifera Earlandia. The distribution of the faunal components within the microbial fabrics suggests that, except for the keratose sponges and some microconchids, most of the metazoans were transported into the microbial framework via wave currents. The presence of both microbialites and metazoan associations were limited to oxygenated settings, suggesting that a factor other than anoxia resulted in a relaxation of ecological constraints following the mass extinction event. It is inferred that the end-Permian mass extinction event decreased the diversity and abundance of metazoans to the point of significantly reducing competition, allowing photosynthesis-based microbial mats to flourish in shallow water settings and resulting in the formation of widespread microbialites.
Management of agricultural soil quality requires fast and cost-efficient methods to identify multiple stressors that can affect soil organisms and associated ecological processes. Here, we propose to use soil protists which have a great yet poorly explored potential for bioindication. They are ubiquitous, highly diverse, and respond to various stresses to agricultural soils caused by frequent management or environmental changes. We test an approach that combines metabarcoding data and machine learning algorithms to identify potential stressors of soil protist community composition and diversity. We measured 17 key variables that reflect various potential stresses on soil protists across 132 plots in 28 Swiss vineyards over 2 years. We identified the taxa showing strong responses to the selected soil variables (potential bioindicator taxa) and tested for their predictive power. Changes in protist taxa occurrence and, to a lesser extent, diversity metrics exhibited great predictive power for the considered soil variables. Soil copper concentration, moisture, pH, and basal respiration were the best predicted soil variables, suggesting that protists are particularly responsive to stresses caused by these variables. The most responsive taxa were found within the clades Rhizaria and Alveolata. Our results also reveal that a majority of the potential bioindicators identified in this study can be used across years, in different regions and across different grape varieties. Altogether, soil protist metabarcoding data combined with machine learning can help identifying specific abiotic stresses on microbial communities caused by agricultural management. Such an approach provides complementary information to existing soil monitoring tools that can help manage the impact of agricultural practices on soil biodiversity and quality.
Water shortage is a serious threat for many societies worldwide. In drylands, water management measures like the construction of reservoirs are affected by eroded sediments transported in the rivers. Thus, the capability of assessing water and sediment fluxes at the river basin scale is of vital importance to support management decisions and policy making. This subject was addressed by the DFG-funded SESAM-project (Sediment Export from large Semi-Arid catchments: Measurements and Modelling). As a part of this project, this thesis focuses on (1) the development and implementation of an erosion module for a meso-scale catchment model, (2) the development of upscaling and generalization methods for the parameterization of such model, (3) the execution of measurements to obtain data required for the modelling and (4) the application of the model to different study areas and its evaluation. The research was carried out in two meso-scale dryland catchments in NE-Spain: Ribera Salada (200 km²) and Isábena (450 km²). Adressing objective 1, WASA-SED, a spatially semi-distributed model for water and sediment transport at the meso-scale was developed. The model simulates runoff and erosion processes at the hillslope scale, transport processes of suspended and bedload fluxes in the river reaches, and retention and remobilisation processes of sediments in reservoirs. This thesis introduces the model concept, presents current model applications and discusses its capabilities and limitations. Modelling at larger scales faces the dilemma of describing relevant processes while maintaining a manageable demand for input data and computation time. WASA-SED addresses this challenge by employing an innovative catena-based upscaling approach: the landscape is represented by characteristic toposequences. For deriving these toposequences with regard to multiple attributes (eg. topography, soils, vegetation) the LUMP-algorithm (Landscape Unit Mapping Program) was developed and related to objective 2. It incorporates an algorithm to retrieve representative catenas and their attributes, based on a Digital Elevation Model and supplemental spatial data. These catenas are classified to provide the discretization for the WASA-SED model. For objective 3, water and sediment fluxes were monitored at the catchment outlet of the Isábena and some of its sub-catchments. For sediment yield estimation, the intermittent measurements of suspended sediment concentration (SSC) had to be interpolated. This thesis presents a comparison of traditional sediment rating curves (SRCs), generalized linear models (GLMs) and non-parametric regression using Random Forests (RF) and Quantile Regression Forests (QRF). The observed SSCs are highly variable and range over six orders of magnitude. For these data, traditional SRCs performed poorly, as did GLMs, despite including other relevant process variables (e.g. rainfall intensities, discharge characteristics). RF and QRF proved to be very robust and performed favourably for reproducing sediment dynamics. QRF additionally excels in providing estimates on the accuracy of the predictions. Subsequent analysis showed that most of the sediment was exported during intense storms of late summer. Later floods yielded successively less sediment. Comparing sediment generation to yield at the outlet suggested considerable storage effects within the river channel. Addressing objective 4, the WASA-SED model was parameterized for the two study areas in NE Spain and applied with different foci. For Ribera Salada, the uncalibrated model yielded reasonable results for runoff and sediment. It provided quantitative measures of the change in runoff and sediment yield for different land-uses. Additional land management scenarios were presented and compared to impacts caused by climate change projections. In contrast, the application for the Isábena focussed on exploring the full potential of the model's predictive capabilities. The calibrated model achieved an acceptable performance for the validation period in terms of water and sediment fluxes. The inadequate representation of the lower sub-catchments inflicted considerable reductions on model performance, while results for the headwater catchments showed good agreement despite stark contrasts in sediment yield. In summary, the application of WASA-SED to three catchments proved the model framework to be a practicable multi-scale approach. It successfully links the hillslope to the catchment scale and integrates the three components hillslope, river and reservoir in one model. Thus, it provides a feasible approach for tackling issues of water and sediment yield at the meso-scale. The crucial role of processes like transmission losses and sediment storage in the river has been identified. Further advances can be expected when the representation of connectivity of water and sediment fluxes (intra-hillslope, hillslope-river, intra-river) is refined and input data improves.
Cosmic-ray neutron sensing (CRNS) is a non-invasive tool for measuring hydrogen pools such as soil moisture, snow or vegetation. The intrinsic integration over a radial hectare-scale footprint is a clear advantage for averaging out small-scale heterogeneity, but on the other hand the data may become hard to interpret in complex terrain with patchy land use.
This study presents a directional shielding approach to prevent neutrons from certain angles from being counted while counting neutrons entering the detector from other angles and explores its potential to gain a sharper horizontal view on the surrounding soil moisture distribution.
Using the Monte Carlo code URANOS (Ultra Rapid Neutron-Only Simulation), we modelled the effect of additional polyethylene shields on the horizontal field of view and assessed its impact on the epithermal count rate, propagated uncertainties and aggregation time.
The results demonstrate that directional CRNS measurements are strongly dominated by isotropic neutron transport, which dilutes the signal of the targeted direction especially from the far field. For typical count rates of customary CRNS stations, directional shielding of half-spaces could not lead to acceptable precision at a daily time resolution. However, the mere statistical distinction of two rates should be feasible.
Cosmic-ray neutron sensing (CRNS) is a non-invasive tool for measuring hydrogen pools such as soil moisture, snow or vegetation. The intrinsic integration over a radial hectare-scale footprint is a clear advantage for averaging out small-scale heterogeneity, but on the other hand the data may become hard to interpret in complex terrain with patchy land use.
This study presents a directional shielding approach to prevent neutrons from certain angles from being counted while counting neutrons entering the detector from other angles and explores its potential to gain a sharper horizontal view on the surrounding soil moisture distribution.
Using the Monte Carlo code URANOS (Ultra Rapid Neutron-Only Simulation), we modelled the effect of additional polyethylene shields on the horizontal field of view and assessed its impact on the epithermal count rate, propagated uncertainties and aggregation time.
The results demonstrate that directional CRNS measurements are strongly dominated by isotropic neutron transport, which dilutes the signal of the targeted direction especially from the far field. For typical count rates of customary CRNS stations, directional shielding of half-spaces could not lead to acceptable precision at a daily time resolution. However, the mere statistical distinction of two rates should be feasible.
Influence of the Main Border Faults on the 3D Hydraulic Field of the Central Upper Rhine Graben
(2019)
The Upper Rhine Graben (URG) is an active rift with a high geothermal potential. Despite being a well-studied area, the three-dimensional interaction of the main controlling factors of the thermal and hydraulic regime is still not fully understood. Therefore, we have used a data-based 3D structural model of the lithological configuration of the central URG for some conceptual numerical experiments of 3D coupled simulations of fluid and heat transport. To assess the influence of the main faults bordering the graben on the hydraulic and the deep thermal field, we carried out a sensitivity analysis on fault width and permeability. Depending on the assigned width and permeability of the main border faults, fluid velocity and temperatures are affected only in the direct proximity of the respective border faults. Hence, the hydraulic characteristics of these major faults do not significantly influence the graben-wide groundwater flow patterns. Instead, the different scenarios tested provide a consistent image of the main characteristics of fluid and heat transport as they have in common: (1) a topography-driven basin-wide fluid flow perpendicular to the rift axis from the graben shoulders to the rift center, (2) a N/NE-directed flow parallel to the rift axis in the center of the rift and, (3) a pronounced upflow of hot fluids along the rift central axis, where the streams from both sides of the rift merge. This upflow axis is predicted to occur predominantly in the center of the URG (northern and southern model area) and shifted towards the eastern boundary fault (central model area).
Silicon is a beneficial element for many plants and is deposited in plant tissue as amorphous bio-opal called phytoliths. The biochemical processes of silicon uptake and precipitation induce isotope fractionation: the mass-dependent shift in the relative abundances of the stable isotopes of silicon. At the bulk scale, delta Si-30 ratios span from -2 to +6 parts per thousand. To further constrain these variations in situ, at the scale of individual phytolith fragments, we used femtosecond laser ablation multi-collector inductively coupled plasma-mass spectrometry (fsLA-MC-ICP-MS). A variety of phytoliths from grasses, trees and ferns were prepared from plant tissue or extracted from soil. Good agreement between phytolith delta Si-30 ratios obtained by bulk solution MC-ICP-MS analysis and in situ isotope ratios from fsLA-MC-ICP-MS validates the method. Bulk solution analyses result in at least twofold better precision for delta Si-30 (2s on reference materials <= 0.11 parts per thousand) over that found for the means of in situ analyses (2s typically <= 0.24 parts per thousand). We find that bushgrass, common reed and horsetail show large internal variations up to 2 parts per thousand in delta Si-30, reflecting the various pathways of silicon from soil to deposition. Femtosecond laser ablation provides a means to identify the underlying processes involved in the formation of phytoliths using silicon isotope ratios.
Year-to-year variations in crop yields can have major impacts on the livelihoods of subsistence farmers and may trigger significant global price fluctuations, with severe consequences for people in developing countries. Fluctuations can be induced by weather conditions, management decisions, weeds, diseases, and pests. Although an explicit quantification and deeper understanding of weather-induced crop-yield variability is essential for adaptation strategies, so far it has only been addressed by empirical models. Here, we provide conservative estimates of the fraction of reported national yield variabilities that can be attributed to weather by state-of-the-art, process-based crop model simulations. We find that observed weather variations can explain more than 50% of the variability in wheat yields in Australia, Canada, Spain, Hungary, and Romania. For maize, weather sensitivities exceed 50% in seven countries, including the United States. The explained variance exceeds 50% for rice in Japan and South Korea and for soy in Argentina. Avoiding water stress by simulating yields assuming full irrigation shows that water limitation is a major driver of the observed variations in most of these countries. Identifying the mechanisms leading to crop-yield fluctuations is not only fundamental for dampening fluctuations, but is also important in the context of the debate on the attribution of loss and damage to climate change. Since process-based crop models not only account for weather influences on crop yields, but also provide options to represent human-management measures, they could become essential tools for differentiating these drivers, and for exploring options to reduce future yield fluctuations.
Thermokarst lakes cover nearly one fourth of ice-rich permafrost lowlands in the Arctic. Sediments from an athalassic subsaline thermokarst lake on Herschel Island (69°36′N; 139°04′W, Canadian Arctic) were used to understand regional changes in climate and in sediment transport, hydrology, nutrient availability and permafrost disturbance. The sediment record spans the last ~ 11,700 years and the basal date is in good agreement with the Holocene onset of thermokarst initiation in the region. Electrical conductivity in pore water continuously decreases, thus indicating desalinization and continuous increase of lake size and water level. The inc/coh ratio of XRF scans provides a high-resolution organic-carbon proxy which correlates with TOC measurements. XRF-derived Mn/Fe ratios indicate aerobic versus anaerobic conditions which moderate the preservation potential of organic matter in lake sediments. The coexistence of marine, brackish and freshwater ostracods and foraminifera is explained by (1) oligohaline to mesohaline water chemistry of the past lake and (2) redeposition of Pleistocene specimens found within upthrusted marine sediments around the lake. Episodes of catchment disturbance are identified when calcareous fossils and allochthonous material were transported into the lake by thermokarst processes such as active-layer detachments, slumping and erosion of ice-rich shores. The pollen record does not show major variations and the pollen-based climate record does not match well with other summer air temperature reconstructions from this region. Local vegetation patterns in small catchments are strongly linked to morphology and sub-surface permafrost conditions rather than to climate. Multidisciplinary studies can identify the onset and life cycle of thermokarst lakes as they play a crucial role in Arctic freshwater ecosystems and in the global carbon cycle of the past, present and future.
Phytoliths in particulate matter released by wind erosion on arable land in La Pampa, Argentina
(2022)
Silicon (Si) is considered a beneficial element in plant nutrition, but its importance on ecosystems goes far beyond that. Various forms of silicon are found in soils, of which the phytogenic pool plays a decisive role due to its good availability. This Si returns to the soil through the decomposition of plant residues, where they then participate in the further cycle as biogenic amorphous silica (bASi) or so-called phytoliths. These have a high affinity for water, so that the water holding capacity and water availability of soils can be increased even by small amounts of ASi. Agricultural land is a considerable global dust source, and dust samples from arable land have shown in cloud formation experiments a several times higher ice nucleation activity than pure mineral dust. Here, particle sizes in the particulate matter fractions (PM) are important, which can travel long distances and reach high altitudes in the atmosphere. Based on this, the research question was whether phytoliths could be detected in PM samples from wind erosion events, what are the main particle sizes of phytoliths and whether an initial quantification was possible.Measurements of PM concentrations were carried out at a wind erosion measuring field in the province La Pampa, Argentina. PM were sampled during five erosion events with Environmental Dust Monitors (EDM). After counting and classifying all particles with diameters between 0.3 and 32 mu m in the EDMs, they are collected on filters. The filters were analyzed by Scanning Electron Microscopy and Energy Dispersive X-Ray analysis (SEM-EDX) to investigate single or ensembles of particles regarding composition and possible origins.The analyses showed up to 8.3 per cent being phytoliths in the emitted dust and up to 25 per cent of organic origin. Particles of organic origin are mostly in the coarse dust fraction, whereas phytoliths are predominately transported in the finer dust fractions. Since phytoliths are both an important source of Si as a plant nutrient and are also involved in soil C fixation, their losses from arable land via dust emissions should be considered and its specific influence on atmospheric processes should be studied in detail in the future.
Flood polders are part of the flood risk management strategy for many lowland rivers. They are used for the controlled storage of flood water so as to lower peak discharges of large floods. Consequently, the flood hazard in adjacent and downstream river reaches is decreased in the case of flood polder utilisation. Flood polders are usually dry storage reservoirs that are typically characterised by agricultural activities or other land use of low economic and ecological vulnerability. The objective of this thesis is to analyse hydraulic, environmental and economic impacts of the utilisation of flood polders in order to draw conclusions for their management. For this purpose, hydrodynamic and water quality modelling as well as an economic vulnerability assessment are employed in two study areas on the Middle Elbe River in Germany. One study area is an existing flood polder system on the tributary Havel, which was put into operation during the Elbe flood in summer 2002. The second study area is a planned flood polder, which is currently in the early planning stages. Furthermore, numerical models of different spatial dimensionality, ranging from zero- to two-dimensional, are applied in order to evaluate their suitability for hydrodynamic and water quality simulations of flood polders in regard to performance and modelling effort. The thesis concludes with overall recommendations on the management of flood polders, including operational schemes and land use. In view of future changes in flood frequency and further increasing values of private and public assets in flood-prone areas, flood polders may be effective and flexible technical flood protection measures that contribute to a successful flood risk management for large lowland rivers.
Forests are a key resource serving a multitude of functions such as providing income to forest owners, supplying industries with timber, protecting water resources, and maintaining biodiversity. Recently much attention has been given to the role of forests in the global carbon cycle and their management for increased carbon sequestration as a possible mitigation option against climate change. Furthermore, the use of harvested wood can contribute to the reduction of atmospheric carbon through (i) carbon sequestration in wood products, (ii) the substitution of non-wood products with wood products, and (iii) through the use of wood as a biofuel to replace fossil fuels. Forest resource managers are challenged by the task to balance these multiple while simultaneously meeting economic requirements and taking into consideration the demands of stakeholder groups. Additionally, risks and uncertainties with regard to uncontrollable external variables such as climate have to be considered in the decision making process. In this study a scientific stakeholder dialogue with forest-related stakeholder groups in the Federal State of Brandenburg was accomplished. The main results of this dialogue were the definition of major forest functions (carbon sequestration, groundwater recharge, biodiversity, and timber production) and priority setting among them by the stakeholders using the pair-wise comparison technique. The impact of different forest management strategies and climate change scenarios on the main functions of forest ecosystems were evaluated at the Kleinsee management unit in south-east Brandenburg. Forest management strategies were simulated over 100 years using the forest growth model 4C and a wood product model (WPM). A current climate scenario and two climate change scenarios based on global circulation models (GCMs) HadCM2 and ECHAM4 were applied. The climate change scenario positively influenced stand productivity, carbon sequestration, and income. The impact on the other forest functions was small. Furthermore, the overall utility of forest management strategies were compared under the priority settings of stakeholders by a multi-criteria analysis (MCA) method. Significant differences in priority setting and the choice of an adequate management strategy were found for the environmentalists on one side and the more economy-oriented forest managers of public and private owned forests on the other side. From an ecological perspective, a conservation strategy would be preferable under all climate scenarios, but the business as usual management would also fit the expectations under the current climate. In contrast, a forest manager in public-owned forests or a private forest owner would prefer a management strategy with an intermediate thinning intensity and a high share of pine stands to enhance income from timber production while maintaining the other forest functions. The analysis served as an example for the combined application of simulation tools and a MCA method for the evaluation of management strategies under multi-purpose and multi-user settings with changing climatic conditions. Another focus was set on quantifying the overall effect of forest management on carbon sequestration in the forest sector and the wood industry sector plus substitution effects. To achieve this objective, the carbon emission reduction potential of material and energy substitution (Smat and Sen) was estimated based on a literature review. On average, for each tonne of dry wood used in a wood product substituting a non-wood product, 0.71 fewer tonnes of fossil carbon are emitted into to the atmosphere. Based on Smat and Sen, the calculation of the carbon emission reduction through substitution was implemented in the WPM. Carbon sequestration and substitution effects of management strategies were simulated at three local scales using the WPM and the forest growth models 4C (management unit level) or EFISCEN (federal state of Brandenburg and Germany). An investigation was conducted on the influence of uncertainties in the initialisation of the WPM, Smat, and basic conditions of the wood product sector on carbon sequestration. Results showed that carbon sequestration in the wood industry sector plus substitution effects exceeded sequestration in the forest sector. In contrast to the carbon pools in the forest sector, which acted as sink or source, the substitution effects continually reduced carbon emission as long as forests are managed and timber is harvested. The main climate protection function was investigated for energy substitution which accounted for about half of the total carbon sequestration, followed by carbon storage in landfills. In Germany, the absolute annual carbon sequestration in the forest and wood industry sector plus substitution effects was 19.9 Mt C. Over 50 years the wood industry sector contributed 70% of the total carbon sequestration plus substitution effects.
Physical and hydrological properties of peat as proxies for degradation of South African peatlands
(2018)
The physical and hydrological properties of peat from seven peatlands in northern Maputaland (South Africa) were investigated and related to the degradation processes of peatlands in different hydrogeomorphic settings. The selected peatlands are representative of typical hydrogeomorphic settings and different stages of human modification from natural to severely degraded. Nineteen transects (141 soil corings in total) were examined in order to describe peat properties typical of the distinct hydrogeomorphic settings. We studied degree of decomposition, organic matter content, bulk density, water retention, saturated hydraulic conductivity and hydrophobicity of the peats. From these properties we derived pore size distribution, unsaturated hydraulic conductivity and maximum capillary rise. We found that, after drainage, degradation advances faster in peatlands containing wood peat than in peatlands containing radicell peat. Eucalyptus plantations in catchment areas are especially threatening to peatlands in seeps, interdune depressions and unchannelled valley bottoms. All peatlands and their recharge areas require wise management, especially valley-bottom peatlands with swamp forest vegetation. Blocking drainage ditches is indispensable as a first step towards achieving the restoration of drained peatland areas, and further measures may be necessary to enhance the distribution of water. The sensitive swamp forest ecosystems should be given conservation priority.
The physical and hydrological properties of peat from seven peatlands in northern Maputaland (South Africa) were investigated and related to the degradation processes of peatlands in different hydrogeomorphic settings. The selected peatlands are representative of typical hydrogeomorphic settings and different stages of human modification from natural to severely degraded. Nineteen transects (141 soil corings in total) were examined in order to describe peat properties typical of the distinct hydrogeomorphic settings. We studied degree of decomposition, organic matter content, bulk density, water retention, saturated hydraulic conductivity and hydrophobicity of the peats. From these properties we derived pore size distribution, unsaturated hydraulic conductivity and maximum capillary rise. We found that, after drainage, degradation advances faster in peatlands containing wood peat than in peatlands containing radicell peat. Eucalyptus plantations in catchment areas are especially threatening to peatlands in seeps, interdune depressions and unchannelled valley bottoms. All peatlands and their recharge areas require wise management, especially valley-bottom peatlands with swamp forest vegetation. Blocking drainage ditches is indispensable as a first step towards achieving the restoration of drained peatland areas, and further measures may be necessary to enhance the distribution of water. The sensitive swamp forest ecosystems should be given conservation priority.
In this dissertation, I describe the mechanisms involved in magmatic plumbing system establishment and evolution. Magmatic plumbing systems play a key role in determining volcanic activity style and recognizing its complexities can help in forecasting eruptions, especially within hazardous volcanic systems such as calderas. I explore the mechanisms of dike emplacement and intrusion geometry that shape magmatic plumbing systems beneath caldera-like topographies and how their characteristics relate to precursory activity of a volcanic eruption. For this purpose, I use scaled laboratory models to study the effect of stress field reorientation on a propagating dike induced by caldera topography. I construct these models by using solid gelatin to mimic the elastic properties of the earth's crust with a caldera on the surface. I inject water as the magma analog and track the evolution of the experiments through qualitative (geometry and stress evolution) and quantitative (displacement and strain computation) descriptions. The results show that a vertical dike deviates towards and outside of the caldera-like margin due to stress field reorientation beneath the caldera-like topography. The propagating intrusion forms a circumferential-eruptive dike when the caldera-like size is small, whereas a cone sheet develops beneath the large caldera-like topography.
To corroborate the results obtained from the experimental models, this thesis also describes the results of a case study utilizing seismic monitoring data associated with the unrest period of the 2015 phreatic eruption of Lascar volcano. Lascar has a crater with a small-scale caldera-like topography and exhibited long-lasting anomalous evolution of the number of long-period (LP) events preceding the 2015 eruption. I apply seismic techniques to constrain the hypocentral locations of LP events and characterize their spatial distribution, obtaining an image of Lascar's plumbing system. I observe an agreement in shallow hypocentral locations obtained through four different seismic techniques; nevertheless, the cross-correlation technique provides the best results. These results depict a plumbing system with a narrow sub-vertical deep conduit and a shallow hydrothermal system, where most LP events are located. These two regions are connected through an intermediate region of path divergence, whose geometry and orientation likely is influenced by stress reorientation due to topographic effects of the caldera-like crater.
Finally, in order to further enhance the interpretations of the previous case study, the seismic data was analyzed in tandem with a complementary multiparametric monitoring dataset. This complementary study confirms that the anomalous LP activity occurred as a sign of unrest in the preparatory phase of the phreatic eruption. In addition, I show how changes observed in other monitored parameters enabled to detect further signs of unrest in the shallow hydrothermal system. Overall, this study demonstrates that detecting complex geometric regions within plumbing systems beneath volcanoes is fundamental to produce an effective forecast of eruptions that from a first view seem to occur without any precursory activity.
Furthermore, through the development of this research I show that combining methods that include both observations and models allows one to obtain a more precise interpretation of the volcanic processes.
Geothermal fields in subduction-related orogens are closely linked to areas characterized by young magmatic and tectonic activity, both in arc- and back-arc settings. The spatio-temporal interaction of Quaternary volcanic complexes with regional extensional and transtensional structures might favor a hydrothermal circuit between meteoric water and magmatic fluids. This study encompasses a kinematic analysis of fault structures from the high-enthalpy system located at the western flank of the Domuyo volcano in Argentina. An analysis of remote sensing data was applied to detect regional lineaments, lineament density, and to identify fracture patterns possibly associated with the different deformational stages documented in the area. These results were combined with detailed fracture analysis and kinematic study of mesoscale faults, as well as existing geological, structural, and geophysical data. The integration suggests that the fluid dynamics of the Domuyo geothermal field are directly conditioned by pre-existing basement structures that were reactivated as normal faults during Pliocene to-Quaternary times. Furthermore, the intensely fracture late Triassic - early Jurassic units are interpreted as the potential level for the reservoir. The fault reactivation was likely associated with the accommodation of regional extension along pre-existing fault structures, and locally enhanced by hydrothermal effects of the Domuyo geothermal field.
Southern Patagonia is a prime example of ongoing oceanic ridge collision and slab-window formation sustained over several million years. The impact of these phenomena on the thermal structure and exhumation of the crust have been mainly assessed with low-temperature thermochronology of bedrock samples. Here, we infer thermal histories from new and existing thermochronological data from the region of most recent ridge collision. In particular, we evaluate the potential far-reaching thermal effects of the evolving slab window, which have previously been considered responsible for patterns of late Miocene reheating associated with back-arc alkaline volcanism. Our model results define protracted cooling since similar to 15 Ma and stepwise exhumation since the late Miocene. The pattern of stepwise exhumation closely matches the onset of Patagonian glaciation at 7 Ma and the successive pulse of glacial incision coeval with neotectonic activity since 3-4 Ma that are also documented by independent geological and geomorphological evidence in the region. Importantly, our findings challenge the recently suggested lack of glacial erosion and incision since 5 Ma in this region. Furthermore, in contrast to previous modelling studies, we find that the available data do not evidence a previously proposed northward-propagating heating event associated with alkaline volcanism. We hypothesize that the anomalous alkaline volcanism in the Patagonian back-arc might be related to trench-orthogonal tears aligned with transform faults in the subducting plate. The substantial differences from the previous modelling procedure on some of the same samples is demonstrated to result from an important lack of convergence in model runs. (C) 2019 Elsevier B.V. All rights reserved.
Human transformation of the Earth’s land surface has far-reaching and important consequences for the functioning of hydrological and hydrochemical processes in watersheds. In nowadays land-use change from forest to pasture is a major issue in particular in the tropics. A sustainable management of deforested areas requires an in-depth understanding of the water and nutrient cycle. On this basis we compared the involved hydrological pathways for rainfall to reach streams and the nutrient budgets of a tropical rainforest and a pasture. In addition we studied the links of hydrochemical differences to differences of the relative importance of flowpaths. This study was conducted in the southwestern part of the Brazilian Amazon basin. An intensive hydrological and hydrochemical sampling and monitoring network was set up. The results indicate that the hydrology was modified in many ways due to land-use change. The most important alteration was the increased importance of the fast flowpath overland flow. Solute exports were in particular linked to the increased volume of overland flow that resulted from the land-use change. An additional reason for the increased nutrient exports from the pasture are the high concentrations of these nutrients in pasture overland flow probably as a due to cattle excrements. Tight nutrient cycles with minimal nutrient losses could not be maintained after the land-use change. This study provides the first attempt to quantify the respective nutrient losses.
Influence of land-use change on near-surface hydrological processes : undisturbed forest to pasture
(2010)
Soil compaction that follows the clearing of tropical forest for cattle pasture is associated with lower soil hydraulic conductivity and increased frequency and volume of overland flow. We investigated the frequency of perched water tables, overland flow and stormflow in an Amazon forest and in an adjacent 25-year-old pasture cleared from the same forest. We compared the results with the frequencies of these phenomena estimated from comparisons of rainfall intensity and soil hydraulic conductivity. The frequency of perched water tables based on rainfall intensity and soil hydraulic conductivity was expected to double in pasture compared with forest. This corresponded closely with an approximate doubling of the frequency of stormflow and overland flow in pasture. In contrast, the stormflow volume in pasture increased 17-fold. This disproportional increase of stormflow resulted from overland flow generation over large areas of pasture, while overland flow generation in the forest was spatially limited and was observed only very near the stream channel. In both catchments, stormflow was generated by saturation excess because of perched water tables and near-surface groundwater levels. Stormflow was occasionally generated in the forest by rapid return flow from macropores, while slow return flow from a continuous perched water table was more common in the pasture. These results suggest that deforestation for pasture alters fundamental mechanisms of stormflow generation and may increase runoff volumes over wide regions of Amazonia.
Prolonged dry periods, and increasingly the generation of smoke and dust in partially-deforested regions, can influence the chemistry of rainfall and throughfall in moist tropical forests. We investigated rainfall and throughfall chemistry in a palm-rich open tropical rainforest in the southwestern Brazilian Amazon state of Rondonia, where precipitation averages 2300 mm year(-1) with a marked seasonal pattern, and where the fragmentation of remaining forest is severe. Covering the transition from dry to wet season (TDWS) and the wet season (WS) of 2004-2005, we sampled 42 rainfall events on event basis as well as 35 events on a within-event basis, and measured concentrations of DOC, Na+, K+, Ca2+, Mg2+, NH4+ , Cl-, SO42- , NO3- and pH in rainfall and throughfall. We found strong evidence of both seasonal and within-event solute rainfall concentration dynamics. Seasonal volume-weighted mean (VWMS) concentrations in rainfall of DOC, K+, Ca2+, Mg2+, NH4+ , SO42- and NO3- were significantly higher in the TDWS than the WS, while VWMS concentrations in throughfall were significantly higher for all solutes except DOC. Patterns were generally similar within rain events, with solute concentrations declining sharply during the first few millimeters of rainfall. Rainfall and throughfall chemistry dynamics appeared to be strongly influenced by forest and pasture burning and a regional atmosphere rich in aerosols at the end of the dry season. These seasonal and within-event patterns of rainfall and throughfall chemistry were stronger than those recorded in central Amazonia, where the dry season is less pronounced and where regional deforestation is less severe. Fragmentation and fire in Rondonia now appear to be altering the patterns in which solutes are delivered to remaining moist tropical forests.
The replacement of undisturbed tropical forest with cattle pasture has the potential to greatly modify the hydrology of small watersheds and the fluxes of solutes. We examined the fluxes of water, Cl-, NO3--N: SO42---S, NH4+-N, Na+, K+, Mg2+ and Ca2+ in different flow paths in similar to 1 ha catchments of undisturbed open tropical rainforest and a 20 year-old pasture established from forest in the southwestern Brazilian Amazon state of Rondonia. Storm flow discharge was 18% of incident rainfall in pasture, but only 1% in forest. Quickflow predominated over baseflow in both catchments and in both wet and dry seasons. In the pasture, groundwater and quickflow were important flow paths for the export of all solutes. In the forest, quickflow was important for NO3--N export, but all other solutes were exported primarily by groundwater outflow. Both catchments were sinks for SO42--S and Ca2+, and sources of Na+. The pasture catchment also lost K+ and Mg2+ because of higher overland flow frequency and volume and to cattle excrement. These results show that forest clearing dramatically influences small watershed hydrology by increasing quickflow and water export to streams. They also indicate that tropical forest watersheds are highly conservative for most solutes but that pastures continue to lose important cations even decades after deforestation and pasture establishment.